Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,705,985	0	T	C	21.5%	173.7 / 43.4	80	intergenic (‑195/+113)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base C (16/1);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 5.60e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATC  >  NC_000913/3705887‑3706078                                                                                                   |                                                                                                   gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                   >  1:411160/1‑100 (MQ=255)    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTCTcc                                                                                                >  2:987141/1‑100 (MQ=255)       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                             <  1:1010741/100‑1 (MQ=255)       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                             <  1:593242/100‑1 (MQ=37)        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  1:106762/100‑1 (MQ=255)           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                         >  2:1439212/1‑100 (MQ=255)             cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                      <  1:1211545/100‑1 (MQ=255)               cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                                    <  1:267284/100‑1 (MQ=255)                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                   <  1:478482/100‑1 (MQ=255)                  cccGCAGAGGCAAAGCGCCACCGGGTATAAAGATTGTTGCCACAGTGAGATTTGTTTATGCCGGATGCGGCGTTAACGACTTATCAGGCCTAC‑AACATCg                                                                                 <  1:1101272/98‑1 (MQ=18)                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                 <  1:321228/100‑1 (MQ=255)                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                 <  2:371093/100‑1 (MQ=255)                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                               <  2:887462/100‑1 (MQ=255)                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                               <  1:324218/100‑1 (MQ=255)                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                              <  2:604494/100‑1 (MQ=255)                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                             >  2:232380/1‑100 (MQ=255)                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                             >  2:1527325/1‑100 (MQ=255)                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                        <  1:202366/100‑1 (MQ=255)                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                       >  2:690466/1‑100 (MQ=255)                             acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGATGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                      <  1:550814/100‑1 (MQ=255)                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGCAAATTCAAt                                                                     >  1:726380/1‑100 (MQ=255)                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                     <  1:1279620/100‑1 (MQ=255)                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                     >  2:976319/1‑100 (MQ=255)                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                     <  2:838362/100‑1 (MQ=255)                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                >  2:1148658/1‑100 (MQ=255)                                    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                               <  1:226481/100‑1 (MQ=255)                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                           >  2:1485127/1‑100 (MQ=255)                                            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑At‑tt                                                   >  1:1099855/1‑100 (MQ=255)                                             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑At‑ttt                                                  >  1:164209/1‑100 (MQ=255)                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCg                                                <  2:849645/100‑1 (MQ=255)                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCg                                                >  1:789106/1‑100 (MQ=255)                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCg                                                <  1:976327/100‑1 (MQ=255)                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGt                                               >  1:1486182/1‑100 (MQ=255)                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTa                                              >  1:251155/1‑100 (MQ=255)                                                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAgg                                            <  2:1187436/100‑1 (MQ=255)                                                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAgg                                            <  2:569334/100‑1 (MQ=255)                                                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAgg                                            >  1:847838/1‑100 (MQ=255)                                                      ttaaCAAAGTGCGCTTTGTTTATGCAGGATTCGTCGTTAACGCCTGATCCGGCCTAC‑AAAATCGTTCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGcc                                          <  2:153546/96‑1 (MQ=17)                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTg                                        >  2:1325421/1‑100 (MQ=255)                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGa                                       >  2:1526449/1‑100 (MQ=255)                                                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGAt                                      >  2:1348016/1‑100 (MQ=255)                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATa                                     <  1:294522/100‑1 (MQ=255)                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATa                                     >  2:131025/1‑100 (MQ=255)                                                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATaa                                    >  2:987632/1‑100 (MQ=255)                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGc                                  <  1:273675/100‑1 (MQ=255)                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCg                                 >  2:139166/1‑100 (MQ=255)                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCgg                                <  2:563755/100‑1 (MQ=255)                                                                   cTTTGTTTATGGCGGATGCGGCGTTAAAGCCTTATCCGGCCTAG‑AAAAGCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑GTTCGTAGGCCTGATAAGCGGAgcg                            <  2:1417513/100‑1 (MQ=17)                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAgcg                            >  1:44386/1‑100 (MQ=255)                                                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCa                          <  1:103631/100‑1 (MQ=255)                                                                       ttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                        >  1:1572164/7‑100 (MQ=18)                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                        >  2:485241/1‑100 (MQ=255)                                                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                       <  1:232381/100‑1 (MQ=255)                                                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAg                      >  1:260919/6‑100 (MQ=18)                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                      >  1:358572/1‑100 (MQ=255)                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAgg                     >  1:1299694/5‑100 (MQ=21)                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑‑AC‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGc                    >  2:649396/3‑100 (MQ=18)                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCA‑ATTCAA‑CTT‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                    >  1:1146109/2‑100 (MQ=18)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑‑AG‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGc                    >  1:249526/2‑99 (MQ=17)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑‑AT‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                   >  2:1130087/2‑100 (MQ=17)                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑‑GA‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGCa                   >  2:461992/2‑100 (MQ=18)                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                  >  1:1317673/1‑100 (MQ=255)                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑‑AT‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                  >  1:993691/1‑99 (MQ=21)                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                >  1:1230649/9‑100 (MQ=18)                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑‑AG‑CCACGTAGGCCTGATAAGCGTAGCGCATCATGCagt                 >  1:947137/9‑98 (MQ=17)                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt               >  1:585039/7‑100 (MQ=18)                                                                                ggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt               >  1:399837/7‑100 (MQ=18)                                                                                ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt               <  1:134947/100‑1 (MQ=255)                                                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt              >  1:831241/6‑100 (MQ=18)                                                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAATCCTTGCGAATTCAATATATTGTATGG‑‑AA‑ATTC‑‑AGGCCTGATAAGCGTAGCGCATCAGGCttttcac             >  1:320365/6‑97 (MQ=37)                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg             <  1:345541/100‑1 (MQ=255)                                                                                   taaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑‑AC‑GAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaactacgt            >  1:1515347/4‑93 (MQ=25)                                                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc            >  2:1003362/1‑100 (MQ=255)                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt         <  2:505645/100‑1 (MQ=255)                                                                                       gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATCAGCGGAGCGCATCAGGCCTTttttccttt        >  2:131786/1‑100 (MQ=255)                                                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg       <  1:690469/100‑1 (MQ=255)                                                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg       <  1:987149/100‑1 (MQ=255)                                                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc     >  2:391425/1‑100 (MQ=255)                                                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa    <  2:1300810/100‑1 (MQ=255)                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc  <  1:592866/100‑1 (MQ=255)                                                                                                   |                                                                                                   GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATC  >  NC_000913/3705887‑3706078

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