Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	213,341	0	A	G	44.6%	111.0 / 33.2	85	G337G (GGA→GGG)	tilS	tRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (7/39);  new base G (37/0);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 7.43e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCCG  >  NC_000913/213196‑213491                                                                                                                                                  |                                                                                                                                                            ccTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAgg                                                                                                                                                          >  1:159012/1‑148 (MQ=255)      cccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAgtg                                                                                                                                                         >  1:1329339/1‑145 (MQ=255)      cccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACa                                                                                                                                                    >  2:1003092/1‑149 (MQ=255)       ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGG‑ACAGc                                                                                                                                                  >  1:618907/1‑150 (MQ=255)        cTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGCGGTGT‑aa                                                                                                                                                     >  2:1256308/1‑145 (MQ=255)         tGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGGAgt                                                                                                                                                          >  2:264328/1‑141 (MQ=255)          gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTa                                                                                                                                               <  2:233816/150‑1 (MQ=255)            ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑t                                                                                                                                            <  2:1115886/150‑1 (MQ=255)             tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGT‑AAAggct                                                                                                                                                >  1:312133/1‑143 (MQ=255)                  ttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGT‑ACAGCTTAA‑TGCggg                                                                                                                                       >  2:1253422/1‑149 (MQ=255)                   tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAg                                                                                                                                     <  1:929265/150‑1 (MQ=255)                     ggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGG‑ACAGCTTTA‑AGCgggggg                                                                                                                                    >  1:1135204/1‑149 (MQ=255)                     ggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCg                                                                                                                                  <  2:122492/151‑1 (MQ=255)                      gcgcGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGa                                                                                                                                 <  2:751535/151‑1 (MQ=255)                       cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTGCCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGat                                                                                                                                <  2:1008802/151‑1 (MQ=255)                       cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGG‑AAAGCTTAA‑TGAgggggggg                                                                                                                                  >  1:979101/1‑147 (MQ=255)                          gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGatat                                                                                                                              <  2:460321/150‑1 (MQ=255)                           tttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGG‑AAAACCTAATTGGGGGAGG‑GATAtt                                                                                                                             >  1:582894/1‑150 (MQ=255)                            ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTc                                                                                                                            <  2:647289/150‑1 (MQ=255)                            ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTc                                                                                                                            <  1:231485/150‑1 (MQ=255)                              gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGG‑ACAGCTTAA‑TGCGGGAGGCGGTATTTGcc                                                                                                                         >  2:1121823/1‑151 (MQ=255)                              gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGc                                                                                                                          <  1:245187/150‑1 (MQ=255)                               aaaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGT‑ACAGCTTAA‑AGCGGGGGGGGATATTCGccc                                                                                                                        >  1:429627/1‑151 (MQ=255)                               aaaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGccc                                                                                                                        <  1:879168/151‑1 (MQ=255)                                aaTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGT‑ACAGGTTAA‑TGCGGGAGGGGGTATTCGccc                                                                                                                        >  2:1320070/1‑150 (MQ=255)                                   ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGG‑ACAGCTTAA‑TACGGGGGGGGAGAATcccccccc                                                                                                                     >  2:585758/1‑150 (MQ=255)                                   ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTcc                                                                                                                     <  1:931942/150‑1 (MQ=255)                                   ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCcg                                                                                                                    <  2:769987/151‑1 (MQ=255)                                       cGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTg                                                                                                                <  1:710556/151‑1 (MQ=255)                                        gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCATCAACCGCTGGAATTACCGGCGGGGGTGGGTAGGGTGACA‑CGTTA‑TGCCGGAGGGGATTTTGGCCATccccctg                                                                                                                >  2:396019/1‑145 (MQ=255)                                          tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGG‑ACAGCTTAA‑TGCGGGGGGCGGTATTTGCCCCCcgc                                                                                                                   >  2:1013008/1‑145 (MQ=255)                                           aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACAGCTTAA‑AGCGGGAGGGGATATTTGcccccccc                                                                                                                   >  2:1316196/1‑142 (MQ=255)                                                    aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGT‑ACAGCTTAA‑TACGGGGGGGGGAATTTCCCCCCCCCGCGCGGACGAAGCgg                                                                                                    >  1:521958/1‑150 (MQ=255)                                                    aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGt                                                                                                   <  1:682420/151‑1 (MQ=255)                                                        gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgc                                                                                               <  2:605491/151‑1 (MQ=255)                                                        gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAg                                                                                                <  1:108839/150‑1 (MQ=255)                                                        gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAg                                                                                                <  2:939424/150‑1 (MQ=255)                                                             ggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgc‑gtgc‑g                                                                                        <  2:1065378/151‑1 (MQ=255)                                                              gATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgc‑gtgc‑gt                                                                                       <  2:878770/151‑1 (MQ=255)                                                                  aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGT‑TAAGCCTAA‑TTTGGGGGGGGATATTCGCCCTCCCCGGGCAGACGAAGCGGGCAgccggcggt                                                                                        >  1:549747/1‑141 (MQ=255)                                                                  aaaTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACAGCTTAA‑TGCGGGGGGGGATATTCGCCCCCCGCGTGCAGACGAAGCGGGCAGC‑GGGC‑Ggtt                                                                                     >  2:515015/1‑149 (MQ=255)                                                                    aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCaaa                                                                                 <  2:1194276/151‑1 (MQ=255)                                                                     tCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGT‑ACAACCTAA‑TGCGGGGGGGGATATTCGCCCTCCCCGGGCAGACGAAACGGGCAGC‑GGGC‑GTTTCaa                                                                                  >  2:172800/1‑149 (MQ=255)                                                                       cGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAgcg                                                                              <  2:350837/151‑1 (MQ=255)                                                                           aCCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACAGCTTAA‑TGCGGGGGGGGATATTTGCCCCCCCCGGGCAGACGAAGGGGGCAGC‑GGGG‑GTTTCAAAACGcc                                                                            >  2:184329/1‑149 (MQ=255)                                                                                 cAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGT‑ACAACTTAA‑TACGGGGGGGGGTATTTGCCCCCCCCGGGCAGACGAAGCGGGGCGC‑GGGG‑GGTTTAAAACGCCA‑GGAtt                                                                     >  1:930962/1‑149 (MQ=255)                                                                                 cAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGAGT‑AGAGCTTCA‑TTCCGGGGGGGGGATTTCACCCCCGCGGGCGGAAGAAGCGGGCAGG‑GTGCGGTTTCaaacggccgg                                                                          >  1:577219/1‑139 (MQ=255)                                                                                  aaaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGG‑ACAGCTTTA‑TTCGGGGGGGGGTATTTGCCCCCCCCGTGCAGGCGGAGGGGGCAGC‑GGGC‑GGTTCTAAGGGCCC‑GGGtt                                                                     >  1:178784/1‑148 (MQ=255)                                                                                   aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGAAGGGT‑ACAGCTTAA‑TGCGGGAGGGGATATTCGCCCTCCGCGTGCAGACGAAACGGGCAGC‑GGGC‑GTTTTAAAACGCCC‑gg                                                                        >  2:1066989/1‑144 (MQ=255)                                                                                     gCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACCGCTTTA‑TGCGGGGGGCGGTATTTGCCCCCCGCGTGCCGACGAAGGGGTCCGC‑GGGC‑GTTTTAAAACGCCC‑GGATtgc                                                                   >  2:1243964/1‑147 (MQ=255)                                                                                       gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGG‑ACAACTTTA‑TGCGGGGGGGGGTATTCGCCCTCCCCCGGCCGACGAAGCGGGCAGG‑GGGC‑GGTTCAAAACGCCA‑GGGTTGCTGCat                                                              >  1:922019/1‑150 (MQ=255)                                                                                        aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACAGCTTAA‑TGCGGGGGGCGGTATTTGCCCCCCCCGGGCAGAAGAAGCGGGCAGC‑GTGC‑GTTTCCAAGCGCCC‑GGATTGCTGCat                                                              >  1:186843/1‑149 (MQ=255)                                                                                        aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCata                                                             <  2:614291/150‑1 (MQ=255)                                                                                          aaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGT‑ACAGCTTAA‑TTCGGGGGGGGATATTTGCCCTCCCCCTGCAGACGAAGCGGGCAGC‑GGGC‑GGTTCAAAGCGCCA‑GGGATGCTGCctaat                                                           >  1:385354/1‑148 (MQ=255)                                                                                           aCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGGGGGT‑ACAGCCTTA‑TGCGGGGGGCGGTATTCTCCCCCCCCGTGCCGACGAAACGGGCCGC‑GTGCGGTTTCAAA‑CGCCC‑GGGTTTCTGCCTAtttt                                                         >  1:44819/1‑149 (MQ=255)                                                                                               tgtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGT‑ACAGCCTTA‑TTCGGGGGGGGGTATTCGCCCCCCCCCTGCAGACGAAACGGTCAGC‑GGGC‑GGTTCAAAAGGCCC‑GGGATGGTGCATATTTGCggg                                                     >  1:366123/1‑151 (MQ=255)                                                                                                 tgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGATGGGGAGTGT‑AAAGCTTAA‑TAAGGGGGGGGACATTTGCCCCCCCCGTGCAGAGGAACGGGGCAGC‑GCGC‑GGTTCAAAACGCCA‑GGGTTGGTGCATAATGTCggg                                                     >  2:98559/1‑149 (MQ=255)                                                                                                  gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCg                                                   <  1:87621/150‑1 (MQ=255)                                                                                                   ccGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGG‑ACAACTTAA‑TGCGGGGGGCGATATTTCCCCCCCCCGTGCAGACGAAGCGGGCAGCGGGGG‑GTTTCAAA‑CGCCAGGGTTTGCT‑CATATTTTCgggg                                                    >  1:1047517/1‑147 (MQ=255)                                                                                                    cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGG‑ACAGCTTTA‑TGCGGGGGGGGATATTTGCCCCCCCCGGGCAGACGAAGCGGTCAGC‑GGGC‑GGTTTAAAACGCCC‑GGGTTGCTTCATATTTTTggggg                                                   >  1:751398/1‑146 (MQ=255)                                                                                                    cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGGGGGGAGGGG‑ACAGCTTAA‑TGCGGGGGGCGATATTTGCCCCCCGCCTGCCGACGAAACGGTCAGC‑GGGC‑GGTTTCAAGCGCCC‑GGATTTGTGCAAATTTTCggggg                                                   >  1:253535/1‑146 (MQ=255)                                                                                                    cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTaa                                                <  1:1089204/151‑1 (MQ=255)                                                                                                     gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGG‑ACAGCTTTA‑TTCGGGGGGGGGAATTTGCCCCCCCCCTGCAGACGGAGCGGGCAGC‑GGGG‑GGTTCCAAACGCCA‑GGGTTGTTGCCTATTTTCggggggaa                                                >  1:884836/1‑145 (MQ=255)                                                                                                       ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAAcg                                              <  2:865418/150‑1 (MQ=255)                                                                                                       ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAAcg                                              <  1:47006/150‑1 (MQ=255)                                                                                                       ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAAcg                                              <  1:1116536/150‑1 (MQ=255)                                                                                                        gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACAGCTTAA‑AGCGGGAGGGGGTATTTGCCCCCCCCGTGCAGACGAAGCGGGCAGA‑GGGC‑GGTTCAAAGCGCCA‑CGGTTGCTGCATATTTTCCGGCGttaagg                                             >  2:907020/1‑145 (MQ=255)                                                                                                           gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACAACTTAA‑TGCGGGAGGCGATATTTGCCCCCCCCGGGCAGAAGAAACGGGCAGC‑GGGC‑GGTTTAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGGAAaggcgg                                          >  2:1059493/1‑150 (MQ=255)                                                                                                             cGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACg                                       <  2:665709/151‑1 (MQ=255)                                                                                                                ggCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTaa                                    <  2:85079/151‑1 (MQ=255)                                                                                                                      aaCCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaaa                              <  1:1062227/151‑1 (MQ=255)                                                                                                                       aCCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGGGCAGACGAAGCGGGCAGC‑GGGC‑GGTTCAAAGCGCCC‑GGATTGCTGCATATTGTCGGGCGGAACGGCGGACGGAAGCTaaa                              >  2:969108/1‑150 (MQ=255)                                                                                                                          gCTGGAATTACCGGCGGGGCTGGGGAGGGT‑ACAGCTTAA‑TGCGGGGGGCGGTATTTGCCCTCCCCGTGCAGACGAAGCGGTCAGC‑GGGCGGTTTCAAA‑CGCCC‑GGGTTGCTGCATCTTTTCGGGGGTAAACGCGGGCGGAAACCaaaaaa                           >  1:736006/1‑150 (MQ=255)                                                                                                                               aaTTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTg                     <  1:372163/151‑1 (MQ=255)                                                                                                                               aaTTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTg                     <  2:277961/151‑1 (MQ=255)                                                                                                                                   aCCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCa                  <  1:547018/150‑1 (MQ=255)                                                                                                                                   aCCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCa                  <  1:1003082/150‑1 (MQ=255)                                                                                                                                       gcggGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGc             <  1:711202/151‑1 (MQ=255)                                                                                                                                        cggGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCt            <  1:512062/151‑1 (MQ=255)                                                                                                                                         ggggCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg           >  1:141529/1‑151 (MQ=255)                                                                                                                                             cTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTTGCAAGAGGTggg         >  1:1224294/1‑149 (MQ=255)                                                                                                                                             cTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GCTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg       >  2:253080/1‑151 (MQ=255)                                                                                                                                              tGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTACATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGt      >  2:1140521/1‑151 (MQ=255)                                                                                                                                               gggAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGt      <  2:1173741/150‑1 (MQ=255)                                                                                                                                               gggAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTg     <  1:356102/151‑1 (MQ=255)                                                                                                                                                  aaGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgccg  >  2:596461/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                            CCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGT‑ACAGCTTAA‑TGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGC‑GTGC‑GTTTCAAAGCGCCA‑GGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCCG  >  NC_000913/213196‑213491

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