Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	2,577,283	+G	100%	intergenic (+205/‑224)	fdx → / → purU	ferredoxin VI/formyltetrahydrofolate deformylase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	2,577,283	1	.→G	100.0%	212.9	59	intergenic (+205/‑224)	fdx/purU	ferredoxin VI/formyltetrahydrofolate deformylase
Reads supporting (aligned to +/- strand):  new base (43/14):  ref base (0/2):  total (43/16)

CGGCGTATTTGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.GGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATCGCGCGCCCCGGAGCAA  >  NC_014034/2577183‑2577378                                                                                                           |                                                                                                  cGGCGTATTTGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.c.c                                                                                                    <  2:1959432/101‑1 (MQ=255)  ggCGTATTTGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.g                                                                                                  >  1:690591/1‑101 (MQ=255)  ggCGTATTTGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.g                                                                                                  <  1:4056912/101‑1 (MQ=255)     gTATTTGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.gggg                                                                                               <  2:2992636/101‑1 (MQ=255)      tATTTGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.ggggg                                                                                              <  1:213362/101‑2 (MQ=255)        tttGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.gggggcg                                                                                            >  1:4203170/1‑98 (MQ=255)        tttGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.gggggcg                                                                                            <  1:4144334/101‑4 (MQ=255)        tttGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.gggggcg                                                                                            >  1:2254173/1‑98 (MQ=255)        tttGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.gggggcg                                                                                            >  1:3570023/1‑98 (MQ=255)         ttGCACCAGGAAAAACCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.gggggcgg                                                                                           <  2:3131435/101‑5 (MQ=255)           gCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAg                                                                                          <  2:870035/101‑1 (MQ=255)           gCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAg                                                                                          <  1:1068080/101‑1 (MQ=255)           gCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAg                                                                                          >  2:821925/1‑101 (MQ=255)           gCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAg                                                                                          >  1:1573328/1‑101 (MQ=255)            cACCAGGAAAAGCCAACGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGc                                                                                         >  1:3896665/1‑101 (MQ=255)              ccAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGccc                                                                                       >  2:3336650/1‑101 (MQ=255)               cAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGcccc                                                                                      <  2:129150/101‑1 (MQ=255)               cAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGcccc                                                                                      <  1:3430696/101‑1 (MQ=255)               cAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGcccc                                                                                      <  2:3896411/101‑1 (MQ=255)                 ggAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCc..g                                                                                  <  1:122252/101‑1 (MQ=255)                 ggAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCc..g                                                                                  <  2:1295204/101‑1 (MQ=255)                   aaaaGCCAAAGCCGCCCCCGGCAGGACCGAAGG.CCCCGGGCCCCCGGGCCCC.TTTTATCTTGCCCAAAATACCCC.GGGGGT.C.CGGGGGCGGAGCCCCc..ggc                                                                                <  2:3989661/101‑1 (MQ=25)                    aaagccaaagccGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCc..ggcg                                                                               <  1:3336887/101‑1 (MQ=255)                         caaagccGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCg                                                                          >  2:3055833/1‑101 (MQ=255)                           agcccGCCCCCGGCAGG.GCGAAGGCCCCCAGCCCCCCGGCCCCT.TTTTATCTGCCCAAAAATACCCG.GGGGGC.CGG.GGGGCGGCCCCCCC..GGCGGGTCGcc                                                                        <  2:3010393/98‑1 (MQ=2)                               cGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGGCGCCgggg                                                                    >  1:4304899/1‑101 (MQ=255)                               cGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGGCGCCgggg                                                                    >  1:2729857/1‑101 (MQ=255)                                gTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGc                                                                   >  1:2058156/1‑101 (MQ=255)                                gTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGGCGCCGGGGc                                                                   >  1:3557623/1‑101 (MQ=255)                                  ccccGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGa                                                                 >  2:447409/1‑101 (MQ=255)                                     cgggcggcccGAGGC.CCTCAGGCCCCCGGCCCCC.TTTTATCTTGCCAAAAAGCCCCG.GGGGGG.CGG.GGGGCGGACCCCCC..GGCGGGTCGCCGGGGCGAAGc                                                              <  1:1541954/93‑1 (MQ=16)                                     cGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGGCGCCGGGGGGAAGc                                                              >  1:490123/1‑101 (MQ=25)                                        cAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCcgc                                                           >  2:2446149/1‑101 (MQ=255)                                        cAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGGCGCCGGGGCGAAGCcgc                                                           >  1:2910297/1‑101 (MQ=255)                                         aGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCcgcg                                                          >  2:63298/1‑101 (MQ=255)                                          ggACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGCAGCcgccg                                                         >  2:2016965/1‑99 (MQ=255)                                          ggACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCgcgg                                                         >  2:1541021/1‑101 (MQ=255)                                               gAAGG.CCTCAGGCCCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGaaa                                                    <  1:999735/101‑1 (MQ=255)                                               gAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGaaa                                                    >  1:198304/1‑101 (MQ=255)                                                aaGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGGCGCCGGGGGGAAGCCCGGGGGaaaa                                                   >  2:81202/1‑100 (MQ=21)                                                   g.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATccc                                                >  1:2385019/1‑101 (MQ=255)                                                     ccTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATcccg                                               <  2:2902608/101‑1 (MQ=255)                                                         aGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCCGCGAAATcccggccc                                           >  2:2126113/1‑101 (MQ=255)                                                            aCCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATcccggccccgg                                        >  2:3657705/1‑101 (MQ=255)                                                              cccGGGCCCCTTTTTACCTTGCCAAAAATATCCC.GGGGGC.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATcccggccccggc                                       <  2:4304630/101‑1 (MQ=18)                                                               ccGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGC.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCaa                                     <  2:2910076/101‑1 (MQ=255)                                                                cGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCCGCAAg                                    >  2:3792255/1‑101 (MQ=255)                                                                 gggCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCCAgc                                   >  2:3742298/1‑101 (MQ=255)                                                                 cggCCCC.TTTTATCCTTCCCAAAAGATCCC.GGGGGTGC.G.GGGGCGGGGCCCCC..GGGGGGTCCCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAgc                                   <  2:607427/100‑1 (MQ=25)                                                                     ccc.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGACGGCGGGGAGAAGGCCCGGCGGAAAcccggccccggggagacccc                               >  1:952614/1‑92 (MQ=17)                                                                      cc.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTg                              >  1:187956/1‑101 (MQ=255)                                                                      cc.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTg                              >  1:4302177/1‑101 (MQ=255)                                                                       c.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGGAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCttg                             >  1:3665279/1‑99 (MQ=255)                                                                       c.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTgg                             >  2:2175461/1‑101 (MQ=255)                                                                         ttttAT.TTCCCCAAAAAACCCG.GGGGGC.G.C.GGGGGGGGGCCCCCCGGGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTggg                            <  1:1017098/101‑1 (MQ=16)                                                                          tttATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGGGCCTGGGc                           >  2:875872/1‑101 (MQ=255)                                                                          tttATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGc                           >  1:678563/1‑101 (MQ=255)                                                                           ttATCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCg                          >  2:3661253/1‑101 (MQ=255)                                                                              tCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGa                       >  2:661937/1‑101 (MQ=255)                                                                              tCTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGa                       <  1:2855786/101‑1 (MQ=255)                                                                               cTTTCCCAAAAAACCCCGGGGGGG.C.G.GGGGGGGCGCCCCCCGGGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCggg                        <  1:3786514/101‑1 (MQ=18)                                                                               cTTGCCCAAAATGTCCT.GGGGGT.C.CGGGGGCGGGGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCTAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGaa                      <  1:2833426/101‑1 (MQ=25)                                                                               cTTGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGaa                      >  1:2506365/1‑101 (MQ=255)                                                                                ttGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGAAAGCGCCTGGGCGGGAAg                     >  2:1977441/1‑101 (MQ=255)                                                                                ttGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGGCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCCGGGCGGGGAg                     >  1:3009802/1‑101 (MQ=21)                                                                                 tGCCCAAAATATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGa                    <  2:1364871/101‑1 (MQ=255)                                                                                     cAAAATGTCCG.GGGGGT.CGG.GGGGCGGCGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATcgc                <  1:4106799/101‑1 (MQ=25)                                                                                       aaaTATCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATcgcgc              >  1:1173934/1‑101 (MQ=255)                                                                                            tCCC.GGGGGT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATCGCGCGcccc         >  2:1580515/1‑101 (MQ=255)                                                                                                 gggggT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGGGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATCGCGCGGCCCGGAg     >  2:2160053/1‑101 (MQ=37)                                                                                                 gggggT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGGAGATCGCGCGCCCCGGAg     >  2:4230508/1‑101 (MQ=255)                                                                                                 gggggT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGGAGATCGCGCGCCCCGGAg     >  1:177871/1‑101 (MQ=255)                                                                                                 gggggT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATCGCGCGCCCCgggg     >  1:68574/1‑99 (MQ=255)                                                                                                    ggT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATCGCGCGCCCCGGAGCaa  >  2:562666/1‑101 (MQ=255)                                                                                                    ggT.C.CGGGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATCGCGCGCCCCGGAGCaa  >  1:3915274/1‑101 (MQ=255)                                                                                                           |                                                                                                  CGGCGTATTTGCACCAGGAAAAGCCAAAGCCGTCCCCGGCAGGACCGAAGG.CCTCAGGACCCCGGGCCCC.TTTTATCTTGCCCAAAATATCCC.GGGGGT.C.C.GGGGCGGAGCCCCC..GGCGGGTCGCCGGGGCGAAGCCGCGGCGAAATCCCGGCCCCGGCAAGCGCCTGGGCGGGAAGATCGCGCGCCCCGGAGCAA  >  NC_014034/2577183‑2577378 Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG