Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,878,536	G→A	Y528Y (TAC→TAT)	gyrB ←	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,878,536	0	G	A	100.0%	193.0 / NA	52	Y528Y (TAC→TAT)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/37);  total (15/37)

ACGGTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAG  >  NC_000913/3878243‑3878828                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                       aCGGTCTGCTCATCATAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCGCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGAGAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTGTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTGGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                                  <  2:413094/300‑1 (MQ=255) aCGGTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                                  <  1:195578/300‑1 (MQ=255) aCGGTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                                  <  2:532344/300‑1 (MQ=255) aCGGTCTGCTCAGCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACg                                                                                                                                                                                                                                                                                                  <  1:414404/300‑1 (MQ=255)   ggTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTg                                                                                                                                                                                                                                                                                                <  1:563179/300‑1 (MQ=255)        gCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTATAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGAGGCTGAGCGATATAGACGTGAAcgc                                                                                                                                                                                                                                                                                           >  2:567642/1‑300 (MQ=255)             tcaGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCa                                                                                                                                                                                                                                                                                      >  1:580565/1‑300 (MQ=255)                     gTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGAACGCGTTCAACGAttt                                                                                                                                                                                                                                                                               >  2:291055/1‑299 (MQ=255)                                                 ataAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAAc                                                                                                                                                                                                                                                <  1:587353/300‑1 (MQ=255)                                                    aGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAAcagc                                                                                                                                                                                                                                             <  1:589822/300‑1 (MQ=255)                                                        ttcagggCATTGCTTTCNTGGATAACGACGCTCCATNTACGATTGTCATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACCTTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGANGCGATAGAAGAAGGTc                                                                                                                                                                                                                                                   <  1:506634/294‑1 (MQ=255)                                                              gCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGaa                                                                                                                                                                                                                                   <  2:393330/300‑1 (MQ=255)                                                              gCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGaa                                                                                                                                                                                                                                   <  1:239936/300‑1 (MQ=255)                                                                aTTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAAtg                                                                                                                                                                                                                                 <  1:581281/300‑1 (MQ=255)                                                                               aaCGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCTTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGa                                                                                                                                                                                                                    <  1:105461/300‑1 (MQ=255)                                                                                acgacgCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGa                                                                                                                                                                                                                    >  2:569219/1‑299 (MQ=255)                                                                                        ccA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCAt                                                                                                                                                                                                           >  1:545980/1‑300 (MQ=255)                                                                                        ccA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATAAAGACGTGAACGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCa                                                                                                                                                                                                            >  2:594734/1‑299 (MQ=255)                                                                                                              tCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCATCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGGAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGc                                                                                                                                                                                         <  2:104768/300‑1 (MQ=255)                                                                                                                tGCGTCTCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCNGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCtg                                                                                                                                                                                       <  2:477191/300‑1 (MQ=255)                                                                                                                tGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCNGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCtg                                                                                                                                                                                       <  2:477137/300‑1 (MQ=255)                                                                                                                      gcgTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGGCGTCCGCATCGGTCATGATGATGATGCTGTGAtc                                                                                                                                                                                  >  2:207632/1‑298 (MQ=255)                                                                                                                      gcgTGGTACTCAGATACCAGTTTCTCTAACACTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCTGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATaa                                                                                                                                                                                 <  2:303809/300‑1 (MQ=255)                                                                                                                         ttgggcttagagaCCAGTTTCTCGAACACTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACTTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGc                                                                                                                                                                              <  2:81934/288‑1 (MQ=255)                                                                                                                              cTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGttt                                                                                                                                                                         >  1:571322/1‑300 (MQ=255)                                                                                                                                   aTACCAGTTTCTATAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGTCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCGTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCg                                                                                                                                                                    <  1:238382/300‑1 (MQ=255)                                                                                                                                   aTACCAGTTCCTCTAACACTCCGCCATCCAATGCAGTTGCAC‑‑TTGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATATAGATCTGGTACTGATCCATCTCTTCTTCGGCTTTAATGTACTATTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGACTTCCTGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGGACGAATGTGCGAACCGTTGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCg                                                                                                                                                                    <  1:238275/300‑1 (MQ=255)                                                                                                                                        aGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGTCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGtt                                                                                                                                                                <  2:386336/299‑1 (MQ=255)                                                                                                                                            tctcTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACt                                                                                                                                                           <  1:341921/300‑1 (MQ=255)                                                                                                                                                tAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAATGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTc                                                                                                                                                          <  2:605793/297‑1 (MQ=255)                                                                                                                                                 aaCGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTTTGATAACGCAGTTTGTCCGGGTTGTACTCGTc                                                                                                                                                       <  2:14003/299‑1 (MQ=255)                                                                                                                                                      ttCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTNGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGAc                                                                                                                                                  <  2:472748/299‑1 (MQ=255)                                                                                                                                                            agccaATGCCGGTGCAC‑‑TTGCGTTGGTGTGCATCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCTGGTTGTACTCGTCACGACCGATAcc                                                                                                                                           <  2:331441/300‑1 (MQ=255)                                                                                                                                                            agccaATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATAcc                                                                                                                                           <  1:41005/300‑1 (MQ=255)                                                                                                                                                                 aTGCCGTTGCAC‑‑TGTCATTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGc                                                                                                                                      <  2:173212/300‑1 (MQ=255)                                                                                                                                                                        tGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAgcgc                                                                                                                               <  1:383107/300‑1 (MQ=255)                                                                                                                                                                          cAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAgcgc                                                                                                                               >  1:33124/1‑298 (MQ=255)                                                                                                                                                                                ggCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGAt                                                                                                                         >  1:462714/1‑300 (MQ=255)                                                                                                                                                                                    ttGGTGTGCAGCGTTGCGCCGTCCAGCTCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCTCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGc                                                                                                                     <  2:289514/300‑1 (MQ=255)                                                                                                                                                                                    ttGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGc                                                                                                                     <  2:489878/300‑1 (MQ=255)                                                                                                                                                                                      ggTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGt                                                                                                                   >  1:350748/1‑300 (MQ=255)                                                                                                                                                                                              gCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGNCTGAGCGATATAGACGTGACCGCGTTCAACGANTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACtt                                                                                                           >  2:469888/1‑300 (MQ=255)                                                                                                                                                                                                                     gagaTCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCg                                                                                    <  1:580858/300‑1 (MQ=255)                                                                                                                                                                                                                                  gATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAAAGAGCATCTTATCGAAGCGCGCTTtct                                                                       >  2:253367/1‑300 (MQ=255)                                                                                                                                                                                                                                                      ttAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTTAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTAccc                                                   <  2:109364/300‑1 (MQ=255)                                                                                                                                                                                                                                                       tAATGTACTATTCCTGCTTACCTTTCTTCACTTTGTACAGCGGCGGCGGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGTTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGATAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTTAGGATTTTAccc                                                   <  2:79289/299‑1 (MQ=255)                                                                                                                                                                                                                                                                           ccTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCCGAATCGCCt                                >  2:299228/1‑298 (MQ=255)                                                                                                                                                                                                                                                                              ttcttcGCTTTTTACAGCTGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGTTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTTAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTc                           <  2:48256/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                       ttGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTa                  >  1:490684/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                ggcggcTGAGCGATAAAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGACTGGTTCTTGCGGTTACGACCCTGc         <  1:594734/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                   ggcTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCtt       <  1:569472/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       gatcGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAACCTTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAATAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCag  <  2:374484/297‑1 (MQ=255)                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                       ACGGTCTGCTCATCAGAAAGGTCAGCTTCCGTCAACGTCGGCTGATAGATAAGCTCTTTCAGCATTGCTTTC‑‑GGATAACGACGCTCCA‑‑TACGATTG‑‑ATCATTTTCTGCGTCGCGTTGTACTCAGATACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCAC‑‑TGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATCGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGA‑‑CGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGTCGACGTCCGCATCGGTCATGATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAG  >  NC_000913/3878243‑3878828

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