Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,540,351	(A)6→5	coding (189/228 nt)	feoA →	ferrous iron transporter, protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,540,346	0	A	.	94.9%	329.0 / 12.7	79	coding (184/228 nt)	feoA	ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base . (35/39);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATCC  >  NC_000913/3540247‑3540441                                                                                                    |                                                                                                ctTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCa                                                                                                 <  2:509880/100‑1 (MQ=255) ctTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCa                                                                                                 <  2:528699/100‑1 (MQ=255)        tGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCaaaaagat                                                                                          >  2:316198/1‑97 (MQ=255)         gTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCaaaaagatc                                                                                         <  1:964916/100‑5 (MQ=255)          ttACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCt                                                                                       <  1:55260/100‑1 (MQ=255)            aCCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTgg                                                                                     <  1:2153743/100‑1 (MQ=255)             ccTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGc                                                                                    >  1:846248/1‑100 (MQ=255)              cTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGcc                                                                                   <  2:496687/100‑1 (MQ=255)                ggCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCtt                                                                                 >  1:1228293/1‑100 (MQ=255)                ggCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCtt                                                                                 >  2:541163/1‑100 (MQ=255)                  ctcctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGCGCCTGGTATTACGC‑AAAAAGATCTGGCCttat                                                                               >  1:8111/1‑100 (MQ=255)                  ctcctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCttat                                                                               >  2:1335534/1‑100 (MQ=255)                    cctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCttatta                                                                             >  2:581229/1‑100 (MQ=255)                     ctctTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCt                                                                                  >  1:423686/1‑94 (MQ=255)                     ctctTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAg                                                                            >  1:2035118/1‑100 (MQ=255)                      tctTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGa                                                                           <  2:2089688/100‑1 (MQ=255)                      tctTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCT‑AAAATAAAAAATAAAGATCTGGCCTTATTAGa                                                                           <  2:1428782/100‑1 (MQ=255)                        ttttAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAACCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAg                                                                         <  1:742034/100‑1 (MQ=255)                         cgatctgtgGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGt                                                                        <  2:639436/95‑1 (MQ=255)                           tAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTgg                                                                      <  2:328236/100‑1 (MQ=255)                           tAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTgg                                                                      <  2:869801/100‑1 (MQ=255)                            aaTGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTgg                                                                      >  2:845243/1‑100 (MQ=255)                            aaTGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGa                                                                     <  1:581230/100‑1 (MQ=255)                              tgtgGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAg                                                                   <  1:1761460/100‑1 (MQ=255)                               gtggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGtagatcg                                                                  >  1:639437/1‑94 (MQ=255)                               gtggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGc                                                                  >  1:450392/1‑100 (MQ=255)                                tggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCg                                                                 <  1:1247959/100‑1 (MQ=255)                                 ggtgCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCgg                                                                >  2:1256650/1‑100 (MQ=255)                                    gcgcGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGttt                                                             >  1:1282348/1‑100 (MQ=255)                                     cgcgTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTc                                                            >  1:1232335/1‑100 (MQ=255)                                      gcgTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTcc                                                           <  2:2081209/100‑1 (MQ=255)                                      gcgTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTcc                                                           >  2:1513137/1‑100 (MQ=255)                                      gcgTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTcc                                                           <  1:148544/100‑1 (MQ=255)                                      gatctgCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTcc                                                           <  2:766360/95‑1 (MQ=255)                                        gTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTg                                                         <  2:1833500/100‑1 (MQ=255)                                         tCGCTCCGCTCGGCGACGCCATTAATATCGAAAGCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTCTTAGGAGTGGAAGCGGTTTCCTg                                                         <  2:464252/100‑1 (MQ=255)                                         tCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGt                                                        >  2:2335735/1‑100 (MQ=255)                                         tCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACG‑TAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGt                                                        <  2:1279700/100‑1 (MQ=255)                                          cGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGtt                                                       <  2:2271579/100‑1 (MQ=255)                                           gCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTccagatc                                                      >  1:766361/1‑95 (MQ=255)                                             tCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAAt                                                    >  1:196094/1‑100 (MQ=255)                                             tCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAAt                                                    <  2:486994/100‑1 (MQ=255)                                              ccACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATa                                                   >  2:1876254/1‑100 (MQ=255)                                               cACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATAc                                                  >  1:2025895/1‑100 (MQ=255)                                                   cGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTg                                              <  2:1483825/100‑1 (MQ=255)                                                        aCCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATaac                                         <  1:1618039/100‑1 (MQ=255)                                                           ccATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATaacaac                                      >  1:1567365/1‑100 (MQ=255)                                                            cATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAacaaca                                     <  1:1006967/100‑1 (MQ=255)                                                            cATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAacaaca                                     >  1:579114/1‑100 (MQ=255)                                                             aTTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAAcaacaa                                    <  2:1228984/100‑1 (MQ=255)                                                             aTTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAAcaacaa                                    <  2:8111/100‑1 (MQ=255)                                                               tCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATg                                  <  2:2146478/100‑1 (MQ=255)                                                                  tatCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaa                               >  2:1071463/1‑100 (MQ=255)                                                                  tatCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaa                               <  2:423684/100‑1 (MQ=255)                                                                   atCGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaaa                              >  1:291229/1‑100 (MQ=255)                                                                     cGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaaaaa                            >  2:1127986/1‑100 (MQ=255)                                                                     cGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaaaaa                            >  2:848071/1‑100 (MQ=255)                                                                     cGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaaaaa                            >  2:2279347/1‑100 (MQ=255)                                                                     cGAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGaaaaaa                            >  2:952560/1‑100 (MQ=255)                                                                      gAAACCCGTCGTGTTAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAAt                           <  1:573582/100‑1 (MQ=255)                                                                      gAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAAt                           >  2:490318/1‑100 (MQ=255)                                                                      gAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAAt                           >  2:2282228/1‑100 (MQ=255)                                                                      gAAACCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAAt                           >  1:1330551/1‑100 (MQ=255)                                                                       aaaCCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAAtt                          <  1:490320/100‑1 (MQ=255)                                                                        aaCCCGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTa                         >  1:1222241/1‑100 (MQ=255)                                                                          cccGTCGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAAc                       <  1:58390/100‑1 (MQ=255)                                                                           ccGTTGTGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAAcc                      <  2:751066/100‑1 (MQ=255)                                                                              tcgtGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCAtt                   <  2:450390/100‑1 (MQ=255)                                                                               cgtGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCAtt                   <  1:637245/100‑1 (MQ=255)                                                                               cgtGTGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTg                  >  1:170957/1‑100 (MQ=255)                                                                               cgtGTGAGCCTAATAATAAGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTg                  <  2:1282345/100‑1 (MQ=255)                                                                                gtgtGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTgg                 <  1:1707702/100‑1 (MQ=255)                                                                                gtgtGAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAATGTAAGCGGTTTTCTGTTTATACGGTGATAAAAACAATGAAGAAATTAACCATTg                  <  2:1449660/99‑1 (MQ=255)                                                                                   tgAGCCTGGTATTACGC‑AAAATGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTCAATACGGTGATAAAAAAAAAAAAAAAATTAAcc                      <  2:1405874/92‑1 (MQ=255)                                                                                   tgAGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCtt              <  2:1567361/100‑1 (MQ=255)                                                                                     aGCCTGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTaa            >  1:248698/1‑100 (MQ=255)                                                                                         tGGTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTgg        <  2:1228290/100‑1 (MQ=255)                                                                                          ggTATTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGt       >  1:1819757/1‑100 (MQ=255)                                                                                             aTTACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAAt    <  1:1256653/100‑1 (MQ=255)                                                                                              ttACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATc   <  2:119493/100‑1 (MQ=255)                                                                                               tACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATcc  >  2:641097/1‑100 (MQ=255)                                                                                               tACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATcc  <  2:1244947/100‑1 (MQ=255)                                                                                               tACGC‑AAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATcc  <  2:1554839/100‑1 (MQ=255)                                                                                                    |                                                                                                CTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTAATTGGTAATCC  >  NC_000913/3540247‑3540441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑