Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	213,341	0	A	G	48.5%	71.5 / 35.0	68	G337G (GGA→GGG)	tilS	tRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (5/30);  new base G (33/0);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.85e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTG  >  NC_000913/213191‑213488                                                                                                                                                       |                                                                                                                                                       aGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGa                                                                                                                                                        <  1:273372/151‑1 (MQ=255)       cTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGTGCTGGGGAggg                                                                                                                                                    >  2:88140/1‑147 (MQ=255)          aCCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAg                                                                                                                                               <  1:212644/151‑1 (MQ=255)               gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCtt                                                                                                                                            <  2:25154/149‑1 (MQ=255)                  tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATg                                                                                                                                        <  2:478839/150‑1 (MQ=255)                   aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGc                                                                                                                                       <  1:548748/150‑1 (MQ=255)                      tttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTAAAGCTTAATTCggg                                                                                                                                    >  1:803982/1‑150 (MQ=255)                        tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAg                                                                                                                                  <  1:840448/150‑1 (MQ=255)                         gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                                 <  2:9756/150‑1 (MQ=255)                          ggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGAGGTTAGCTTaaatagggggggg                                                                                                                               >  2:177955/1‑140 (MQ=255)                           gcgcGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGAGAGGTTTATTCgggggg                                                                                                                                 >  2:300318/1‑148 (MQ=255)                           gcgcGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCgggggg                                                                                                                                 >  1:31660/1‑148 (MQ=255)                           gcgcGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGa                                                                                                                              <  1:19594/151‑1 (MQ=255)                            cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTAACGCTTTATGCgggggg                                                                                                                                 >  2:1037415/1‑147 (MQ=255)                            cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGat                                                                                                                             <  1:436671/151‑1 (MQ=255)                             gcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAAGGTGCAGCCTAAAGCgggggggg                                                                                                                               >  2:892914/1‑146 (MQ=255)                                 ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATGCGGGAGGCGATA‑tt                                                                                                                         >  1:931740/1‑149 (MQ=255)                                 ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACCGATTTATGCGGGGGGGGGTA‑tt                                                                                                                         >  2:288756/1‑149 (MQ=255)                                  tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGGGGGGGGAGGACATCTTAATgagggagggggaattt                                                                                                                         >  1:657084/1‑141 (MQ=255)                                  tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGATA‑tt                                                                                                                         >  2:126816/1‑148 (MQ=255)                                  tGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGGGATA‑tttg                                                                                                                       >  2:558271/1‑148 (MQ=255)                                   gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGc                                                                                                                      <  2:712508/150‑1 (MQ=255)                                      aTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTc                                                                                                                  <  2:415452/151‑1 (MQ=255)                                               tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCCTTAAGCGGGGGGCGGAA‑TTCGCCCCCCCCGTg                                                                                                            >  2:899622/1‑148 (MQ=255)                                                 tCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATTCGGGGGGGGATA‑TTTGCCCTCcgc                                                                                                               >  1:1038332/1‑143 (MQ=255)                                                  cAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCGGTA‑TTCGCCCTCCGCGTGCAg                                                                                                         >  1:987589/1‑148 (MQ=255)                                                        cAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCgg                                                                                                <  1:638136/151‑1 (MQ=255)                                                             gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAg                                                                                            <  1:249106/150‑1 (MQ=255)                                                              tggtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACCGCTTCAAGCGGGGGGGGGTA‑TTTGCCCTCCCCCTGCAGGCGAAGCgg                                                                                                >  2:661874/1‑145 (MQ=255)                                                                gtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtg                                                                                        <  1:759898/151‑1 (MQ=255)                                                                  ggATTGAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtgc                                                                                       <  2:948947/150‑1 (MQ=255)                                                                    aTTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑Gtt                                                                                   <  2:358351/151‑1 (MQ=255)                                                                          tCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGGGGGGACAGGCTTATTCGGGGGGGGGTA‑TTCGCCCTCCGCGTGCAGACGAAGGGGGCCgcgtggggt                                                                                    >  2:386236/1‑141 (MQ=255)                                                                           ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGGGGTA‑TTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGG‑GGTTTaaa                                                                              >  1:535285/1‑149 (MQ=255)                                                                               cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGcc                                                                         <  2:436961/150‑1 (MQ=255)                                                                                  cGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCa                                                                        <  1:911260/148‑1 (MQ=255)                                                                                   gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGGGAGCTTAATGCGCGGGGAGAGATTTCCCCCTCC‑CCTGCTGACGAGGCGGGCCGGGTGCGGTTTCaaacgccgcg                                                                       >  2:151558/1‑142 (MQ=255)                                                                                        aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATtgctg                                                              <  2:788780/151‑1 (MQ=255)                                                                                         aGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATGCGGGGGGGGATA‑TTTGCCCCCCGCGGGCCGACGAAGCGGTCAGGGGGG‑GGTTTAAAACGCCA‑GGGtt                                                                  >  1:315058/1‑146 (MQ=255)                                                                                           cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCat                                                           <  2:490415/151‑1 (MQ=255)                                                                                           cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAGGGGTCAGCTTAAAGTCGGGGGCGGTA‑TTTGCCCTCCGCGTGCCGAAGAAGCGGGCGGCGGGG‑GTTTCCAAACGGCA‑ggggtggttctt                                                           >  1:487576/1‑141 (MQ=255)                                                                                             aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTAAAGCTTAAAGAGGGGGGGGAAA‑TTTGCCCCCCCCGGGCAGACGAAAAGGGCAGGGGGC‑GGTTTTAAACGCCC‑GGGGTGCTGCaaa                                                          >  2:1066153/1‑148 (MQ=255)                                                                                             aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCatat                                                         <  2:898618/151‑1 (MQ=255)                                                                                                 cATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTc                                                     <  1:65039/151‑1 (MQ=255)                                                                                                    tgtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTTAAGCGGGGGGGGCTA‑TTTGCCCCCCGCGTGCAGACGAAGCGGGCAGCGGGC‑GGTTTAAAACGCCA‑GGATTGCTGCATATTGTCggg                                                  >  1:936707/1‑151 (MQ=255)                                                                                                       gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGt                                               <  2:1047637/151‑1 (MQ=255)                                                                                                         cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGAGGCGGTA‑TTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGC‑GGTTTAAAACGCCA‑GGGTTGCTGCATATTGTCGGGCGTa                                              >  1:372051/1‑150 (MQ=255)                                                                                                           tGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTTATTCGGGGGGGGGTA‑TTTTCCCCCCCCGGGCAGAAGAAGCGGGCAGGGCGG‑GGTTTTCAACGCCC‑GGGTTGCCGCATATTTTCggggg                                                >  2:92122/1‑144 (MQ=255)                                                                                                            ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAAcg                                           <  1:288827/150‑1 (MQ=255)                                                                                                             gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAGGGGACAGCTTAATGCGGGAGGCGGTA‑TTTGCCCCCCCCCTGCAGACGAAACGGGCAGCGGGC‑GGTTCAAAACGCCA‑GGATTGCTGCATATTGTCGGGGGTAAcgg                                          >  1:989195/1‑150 (MQ=255)                                                                                                               aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCGTTATGCGGGAGGGGAGA‑TTCGCCCTCCGCGGGCAGACGAAGCGGGCAGCGGGC‑GGTTTAAAAAGCCC‑GGGTTGTTGCATATTTTTGGGGGTGAcgggg                                        >  2:656614/1‑148 (MQ=255)                                                                                                               aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGGGGTA‑TTCGCCCTCCCCGCGCAGACGAAGCGGGCAGCGGGC‑GGTTCAAAACGCCA‑GGATTGCTGCATATTGTCGGGCGGAAcggcgg                                       >  2:776799/1‑151 (MQ=255)                                                                                                                   gTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTAATGCGGGGGGCGGTA‑TTCGCCCTCCGCGGGCAGACGAAGCGGGCAGCGGGG‑GGTTCAAAACGCCC‑GGGTTGGTGCATATTTTCggggggaaaggggg                                       >  1:93671/1‑136 (MQ=255)                                                                                                                   gTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGGGATA‑TTCGCCCTCCGCGGGCAGACGAAGAGGGGAGCGTGG‑GGTTTTAAGCGCCC‑GGGTTGCTGCCTATTTTCCGGGGGAACGGGGGGCgg                                   >  2:532980/1‑150 (MQ=255)                                                                                                                   gTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAGGGGACAGCTTTATGCGGGGGGGGGTA‑TTCCCCCCCCGCGTGCCGAAGGAACGGGCAGCGGGC‑GGTTTAAAACGCCACGGATTG‑TGCATATTGTTGGGCGTGAcggcggg                                      >  2:105757/1‑147 (MQ=255)                                                                                                                    tGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAAAGGACAATCTAATGAGGGAGGGGGCA‑TTCGCCCCCCCCGTGCAGAAGAAAGGGGCACCGGGG‑GGTTCAAAACGCCA‑GGGATGCTGCCAAATGTCGGGGGGTAcggcggg                                      >  2:932739/1‑146 (MQ=255)                                                                                                                     ggCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTa                                  <  2:613641/150‑1 (MQ=255)                                                                                                                         tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGc                               <  1:558399/149‑1 (MQ=255)                                                                                                                          cAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGGAGAGCTTAATTGGGGGGGGGATA‑TTTTCCCCCCGCGGGGAGAAGAAGCGGGCGGCGTGC‑GGTTCAAAGCACCA‑GGATTTGTGCATATTGTCGGGCGGAAACGCGGGCGGAAggt                              >  1:424945/1‑147 (MQ=255)                                                                                                                           aaCCGCTGGAATTACCGGCGGGGCTGGGAAGAGTACAGCTTAATGCGGGAGGGGGAA‑TTCGCCCTCCCCGTGCAGACGGAGCGGGCAGCGGGC‑GGTTCTAAGCGCCA‑GGATTGCTGGATATTGTCGGGCGTAAAGGCGGGCGGAAGGTaa                            >  1:169715/1‑150 (MQ=255)                                                                                                                             ccGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGa                         <  1:817134/151‑1 (MQ=255)                                                                                                                               gCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaa                        >  1:478089/1‑150 (MQ=255)                                                                                                                                          cggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAaga            <  2:1050893/151‑1 (MQ=255)                                                                                                                                           ggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAag             >  1:539064/1‑149 (MQ=255)                                                                                                                                               gggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGAAAGCTAAAGAAAATCTgg                 >  2:162135/1‑141 (MQ=255)                                                                                                                                                ggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTggg      >  2:493519/1‑151 (MQ=255)                                                                                                                                                 gCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGc     <  1:191924/151‑1 (MQ=255)                                                                                                                                                    gggAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg    <  1:934481/149‑1 (MQ=255)                                                                                                                                                    gggAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTg  <  1:59244/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATA‑TTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTG  >  NC_000913/213191‑213488

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