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	92.9%	249.8 / 2.0	87	Y528Y (TAC→TAT)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base G (3/3);  new base A (39/40);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTC  >  NC_000913/3878289‑3878782                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          agataATCTCTTTCAGCAGTGCGTTCGGATAACGACGCTC‑CTTACGGTTGATCATTTTCTGCGTCGCGGTGTACTCAGA‑TTCCAGTTTCGCTAACGCCTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTGCTTCACTTTGTACAGCGGCGGCTGAGCGATATAg                                                                                                                                                                                                                                                        <  2:407902/251‑1 (MQ=255) agataAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAg                                                                                                                                                                                                                                                        <  1:465627/251‑1 (MQ=255)     aaGCTCTTTCAGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCCGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCATGCTTGCCTTTCTTCACTTTTTACAGCGTCGGCTGATCGATATAGACGt                                                                                                                                                                                                                                                    >  1:663503/1‑251 (MQ=255)     aaGCTCTTTCAGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGAGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTAATCCATC‑GCTTCGCAGCCTTTCATGCACTGTTCCTGCATGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGt                                                                                                                                                                                                                                                    >  1:448540/1‑251 (MQ=255)       gCTCTTTCAGCATTGCTTTCGTATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGa                                                                                                                                                                                                                                                  <  2:523838/251‑1 (MQ=255)          ctTTCAGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACcg                                                                                                                                                                                                                                               >  1:282877/1‑251 (MQ=255)            ttCAGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGAcc                                                                                                                                                                                                                                                <  1:342675/248‑1 (MQ=255)            ttCAGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACcgcg                                                                                                                                                                                                                                             >  1:677414/1‑251 (MQ=255)               aGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTc                                                                                                                                                                                                                                          >  1:519035/1‑251 (MQ=255)                 cATTGCTTGCGGTTAACGACG‑TCTCATCCTATTGATCATTTTCTGCGTCGCGTTGTCCTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCaa                                                                                                                                                                                                                                        <  2:260469/251‑1 (MQ=255)                   ttGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACg                                                                                                                                                                                                                                      >  1:222450/1‑251 (MQ=255)                        ttCGGATAACGACGCTC‑CATACGATTGAACATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCACTGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCCTC‑GCTTCGTCGTCTTTACTGTACTGTTCCTGCTTGCCTTACTTCACTTTGTACCGCGGCGGCTGAGCGATCTAGCCGTGACTGCGTTCAACGATTTc                                                                                                                                                                                                                                 >  2:341459/1‑251 (MQ=255)                        ttCGGATAACCTCGTTC‑CATTTGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TCCCAGTTTCTCTAACGGTTCGCCAGCCAATGCCGGTGCACTGGCGGTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTGTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCTTTCAACGAttt                                                                                                                                                                                                                                  <  2:643869/250‑1 (MQ=255)                                   acgCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGCCCGCGTTCAACGATTTCCGGCATCTGAc                                                                                                                                                                                                                      >  2:135199/1‑251 (MQ=255)                                        ccctaatattttatccttttCTGCGTCGCGTGGTACT‑AGAATACCAGTTGCTCCACCGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATa                                                                                                                                                                                                                  <  2:645706/235‑1 (MQ=255)                                             tcGATTCATCAGTTTCTGCGCCGCGTTGTACTCAGA‑TACCAGTTTCTCTAGCGCTTCGCCAGCCTATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGTTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGAGTTCAACGATTTCCGGCATCTGACGAaagaag                                                                                                                                                                                                              <  2:140366/249‑1 (MQ=255)                                              cGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGCTAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGCGCGATCTAGACGTGCCCGCGTTCAACTCTTTCCGGCATCTGCCGTTagaag                                                                                                                                                                                                              >  2:110063/1‑249 (MQ=255)                                                aTTGATCATTTTCTGCGTCGCGTTGTACTCCGA‑TACAAGTTTCTCTAACGCTTCGCCAGCCACTGCCGGTGCACTGGCGTTGGTGTGCAGCTTTGCGCCGTCAAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTATTTACTGTACTGTTCCTGCTTGCCTTTCTTCCCTTTGTCCAGCGGCGGCTGTGCGCTATAAGCGCGACCGCGTTCCACGATTTCCGGCATCTGTCGATAGAAGCAgg                                                                                                                                                                                                          >  2:433680/1‑251 (MQ=255)                                                  tGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGt                                                                                                                                                                                                         <  1:471251/250‑1 (MQ=255)                                                   gATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCa                                                                                                                                                                                                       <  1:508376/251‑1 (MQ=255)                                                                    cgTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGGCAACAGCAGCGTACGAAtg                                                                                                                                                                                      >  2:214128/1‑251 (MQ=255)                                                                     gtggtatccAGA‑TACTAGTTTCCTTACCGCTCGGCCAGACAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTGGCGCCGTCCTGCGCGATAGAGATCTGGTACTGATCCAGC‑TCTGCGGCGTCTGTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTGGTACAGCTGCGGCTGAGCGATATAGACGTGACCGCGTGCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAAtgt                                                                                                                                                                                     <  2:67968/243‑1 (MQ=255)                                                                       tGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGc                                                                                                                                                                                   >  1:249438/1‑251 (MQ=255)                                                                       tGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGc                                                                                                                                                                                   >  1:612257/1‑251 (MQ=255)                                                                       tGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGAAAGAAGAAGGGCAACCGCAGCGTACGAATGTGc                                                                                                                                                                                   >  2:189497/1‑251 (MQ=255)                                                                            tCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGcc                                                                                                                                                                              >  1:447714/1‑251 (MQ=255)                                                                             gaGA‑TCCCAGTTTGTCTAACTCTCCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCTGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGGTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGGTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGcc                                                                                                                                                                              <  2:674264/249‑1 (MQ=255)                                                                             cAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCg                                                                                                                                                                             >  1:602510/1‑251 (MQ=255)                                                                                     cAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGt                                                                                                                                                                     >  1:460873/1‑251 (MQ=255)                                                                                     cAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGt                                                                                                                                                                     <  1:450450/251‑1 (MQ=255)                                                                                       gTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTcc                                                                                                                                                                   <  1:107011/251‑1 (MQ=255)                                                                                       gTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTcc                                                                                                                                                                   <  1:164842/251‑1 (MQ=255)                                                                                       gTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTcc                                                                                                                                                                   <  1:214128/251‑1 (MQ=255)                                                                                         ttCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTAAAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCCTCTGACGATAGAAGAGGGGCAACAGAAGCGTACGCATGGGCGAGCCGT‑CGACGTCCGc                                                                                                                                                                 >  1:22290/1‑251 (MQ=255)                                                                                              tAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGAGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCgg                                                                                                                                                            >  1:468029/1‑251 (MQ=255)                                                                                                  gCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑Cat                                                                                                                                                       >  1:558937/1‑251 (MQ=255)                                                                                                       gccagccaATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑Catg‑atg                                                                                                                                                  >  2:70403/1‑250 (MQ=255)                                                                                                       gccagccaATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑Catg‑atg                                                                                                                                                  >  1:582434/1‑250 (MQ=255)                                                                                                            ccaACGCCGGGGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTCCTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑Catg‑atgatgatg                                                                                                                                            <  2:447714/251‑1 (MQ=255)                                                                                                             caATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGc                                                                                                                                           <  1:24641/251‑1 (MQ=255)                                                                                                                   cGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGATTGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGAt                                                                                                                                     >  1:600918/1‑251 (MQ=255)                                                                                                                            ggCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCGTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTTACCGCGTTCAACGATTTCCGGCAGCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGAGGCTGTGATAACGCAGt                                                                                                                             <  1:617482/250‑1 (MQ=255)                                                                                                                                     gtgcagatttGCGCCGACCCGCGCAATAGCGCTCCGGTGCTGATCTATC‑GCGTCGCCGTCTTTCATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGGCGTTACCGCGCTCAACGCTTTCTGGCATCTGACGCTTGCAGAAGGTCAACAACAGCGGGCTGATGTGCGAGCCGTTCGA‑GTCCGAATCGGT‑CATT‑ATGATGATGATGTGATTACGCCGGTTGTACgggg                                                                                                                   >  2:354771/9‑250 (MQ=255)                                                                                                                                         aGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTa                                                                                                               <  1:360999/251‑1 (MQ=255)                                                                                                                                            gTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTc                                                                                                            >  2:680030/1‑251 (MQ=255)                                                                                                                                               gcgcCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTc                                                                                                         >  1:362595/1‑251 (MQ=255)                                                                                                                                                  ccGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTTTCCGGGTTGTACTCGTCACg                                                                                                      >  1:36741/1‑251 (MQ=255)                                                                                                                                                       cAGCGCGATAGAGATATGGTACTGATCAATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAACGATATAAACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGCGCCGT‑CGACATCCGCATCGGG‑CATG‑ATGATGCTGCTGTGATCACGCAGTTTGTCCGGGTTGTACTCGTCACGGCCgg                                                                                                 >  2:350273/1‑250 (MQ=255)                                                                                                                                                        aGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGAt                                                                                                >  1:301036/1‑251 (MQ=255)                                                                                                                                                        aGCGCGATAGAGAGCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGAt                                                                                                <  1:394119/251‑1 (MQ=255)                                                                                                                                                            ggatatgtgtCTGTTACTGATCC‑TCTGCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCGTCACTTTGTTCAGCGGCGGCTGAGCGATATAGACTGGACCTCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTTCTAATGTGCTAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGGTTGTACGGGTTGTACTCGTCACGACCGATACca                                                                                            <  2:18121/242‑1 (MQ=255)                                                                                                                                                              atagatttCTGGTAATGATCCATA‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTCCAGCGCCGGCTGTGCGATATAGACATGACCGCGTTAAACGATTTCCGGCGTCTGACGATAGAAGATGGTAACCAGCAGCGTACGAATGTGCGAGCCGT‑CGGCGTCAGCATCGGTACA‑G‑ATGATGATGCTGTGCTAACACAGTTTGTCCGGGTTGTACTCGTCAAGACCGAAACcaac                                                                                          >  2:462903/8‑249 (MQ=255)                                                                                                                                                              aTAGAGATCTGGTACTGATCCATC‑GGTTCGTCGGCGTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGGGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑AGGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACcaca                                                                                          <  1:710275/251‑1 (MQ=255)                                                                                                                                                                agagacttGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGc                                                                                        <  2:519035/244‑1 (MQ=255)                                                                                                                                                                               aTCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATa                                                                          >  1:439138/1‑250 (MQ=255)                                                                                                                                                                                    cc‑GCTTCGTGGTCTTTAATGTCCTGTGCCTGCTGGCCTTTCTCCGCTTTGTACAGCGGCGGCCGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAAGGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTc                                                                    <  2:222450/250‑1 (MQ=255)                                                                                                                                                                                     c‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGGACAGCGGCGGCTGAGCGATATAGACGTGCCCGCGTTCAACGCTTTCCGGCATCTGACGATAGCAGAAGGTCCACAGCAGCGTACGAATGTGCGCGCCGT‑CGCCGTCCGCATCGGA‑CATGCAT‑ATGATGCTGTGCTAACGCAGTTTGTCCGGGTTGTACTCGTCACGGCCGATACCCCAGCCAAGAGCGGTGGTAAGCGTcc                                                                   >  2:352628/1‑250 (MQ=255)                                                                                                                                                                                           cggcgtcTTTAATGTACTGTTCCTGCTTGCGTTTCTTCACTGTGTACAGCGGCGGCTGAGCGATGGAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCAcg                                                              <  2:9496/248‑2 (MQ=255)                                                                                                                                                                                            gtcgtcTTTGATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACtt                                                             >  1:251101/1‑251 (MQ=255)                                                                                                                                                                                            gtcgtcTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGGCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCCTCGGT‑CCTG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCCCAGCCAAGCGCGTGGGTAAGCGTCGCCCCtt                                                             >  2:38420/1‑251 (MQ=255)                                                                                                                                                                                                   ttAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTga                                                        <  1:135199/249‑1 (MQ=255)                                                                                                                                                                                                     caTGTACTGGTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCCGAGCGCTGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGGCGATAGAAGAAGGGCAACCGCAGCGTACGAATGTGCGAGCCGT‑CGGCGTCCGCATCGGC‑CATG‑ATGCTGATGCTGTGATAACGCAGTTTGTCCGGGTTATACTc                                                                                                            >  2:337389/2‑195 (MQ=255)                                                                                                                                                                                                     aaTGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCCGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGC‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTATACTc                                                                                                            <  1:337389/195‑1 (MQ=255)                                                                                                                                                                                                       ggTACTGTACCTGCTTTCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGAGATAGACGTGACCGCGTCCAACGATTTCCGGCATCTGACGATAGAATAAGGTCAACAGCAGCGTACGAATGTGCGAGGCGT‑CGACGTCCGCATCGGT‑CATT‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCTGTGATAAGCGTCGCCACTTCCTGAGAagag                                                  <  1:286958/250‑1 (MQ=255)                                                                                                                                                                                                        gTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGagaag                                                    >  1:147460/1‑248 (MQ=255)                                                                                                                                                                                                        gTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAagaga                                                 <  2:243111/251‑1 (MQ=255)                                                                                                                                                                                                        gTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAagag                                                  <  1:70403/250‑1 (MQ=255)                                                                                                                                                                                                            tGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCAt                                             >  1:452314/1‑251 (MQ=255)                                                                                                                                                                                                             gtcccggCTTGCCTTTCTTCGCTTTGTACAGCGGAGGGTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCa                                              <  2:159607/243‑1 (MQ=255)                                                                                                                                                                                                              ttCCTGCTTTCCTTTCTTCCCTTTGTACAGCGGCGGCTGAGCGATATGGACGTGACCGCGTTCAACGATTGCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGCCGCCACTTCCTGAGAAGAGAGCATc                                            <  2:293645/250‑1 (MQ=255)                                                                                                                                                                                                               tCCTGCTTGCCTTTTTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTAAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCt                                           <  2:142676/250‑1 (MQ=255)                                                                                                                                                                                                                ccTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCtt                                          <  1:303647/250‑1 (MQ=255)                                                                                                                                                                                                                ccTGCGTGCCTTTCTTCCCTTTGTACAGCGGCGGCTTAGCGATATAGACGGGACCGCGTTCAACGATTTCCGGCATCGGACGAGAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTa                                         <  2:602510/251‑1 (MQ=255)                                                                                                                                                                                                                   gCTTGCCTTTCTTCCCTTTGTACAGCGTCGGCTTAGTGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTAt                                        <  1:201716/249‑1 (MQ=255)                                                                                                                                                                                                                      tGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTAGCGaa                                    <  1:665766/250‑1 (MQ=255)                                                                                                                                                                                                                        ccTTTCTTCCCTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAgcg                                 <  1:188622/251‑1 (MQ=255)                                                                                                                                                                                                                         cTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGGGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAgcgc                                >  1:391537/1‑251 (MQ=255)                                                                                                                                                                                                                           ttcttcACTTTGTACAGCGGCGGCTGGGCGATATAGACGTGACCGCGTTCAACGCTTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGAGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAgcgcgc                              >  2:591991/1‑251 (MQ=255)                                                                                                                                                                                                                           ttcttcACTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGGGATAAGCGTCGCCCATTCCTGAGAAGAGAGCATCTTATCGAAgcgcgc                              >  2:35691/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccTTTGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTtctc                        <  1:53632/250‑1 (MQ=255)                                                                                                                                                                                                                                   tttGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTCTCGAAGCGCGCTTTCTCg                       >  1:82893/1‑250 (MQ=255)                                                                                                                                                                                                                                    ttGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGcc                                                                                       >  1:193157/1‑185 (MQ=255)                                                                                                                                                                                                                                    ttGTACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGcc                                                                                       <  2:193157/185‑1 (MQ=255)                                                                                                                                                                                                                                       tACAGCGGCGGCTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGtt                  >  1:65872/1‑251 (MQ=255)                                                                                                                                                                                                                                           gctgcggcgtaGCGATGGAGTCGTGACCGCGTGTACCGATTGCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCTTACGAATTTGCTAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATACCGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAg               <  2:22290/240‑1 (MQ=255)                                                                                                                                                                                                                                                ggcTGAGCGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGAtttt         >  1:67173/1‑251 (MQ=255)                                                                                                                                                                                                                                                       cGATATAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTc  <  1:697820/251‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          AGATAAGCTCTTTCAGCATTGCTTTCGGATAACGACGCTC‑CATACGATTGATCATTTTCTGCGTCGCGTTGTACTCAGA‑TACCAGTTTCTCTAACGCTTCGCCAGCCAATGCCGGTGCACTGGCGTTGGTGTGCAGCGTTGCGCCGTCCAGCGCGATAGAGATCTGGTACTGATCCATC‑GCTTCGTCGTCTTTAATGTACTGTTCCTGCTTGCCTTTCTTCACTTTGTACAGCGGCGGCTGAGCGATGTAGACGTGACCGCGTTCAACGATTTCCGGCATCTGACGATAGAAGAAGGTCAACAGCAGCGTACGAATGTGCGAGCCGT‑CGACGTCCGCATCGGT‑CATG‑ATGATGATGCTGTGATAACGCAGTTTGTCCGGGTTGTACTCGTCACGACCGATACCACAGCCAAGCGCGGTGATAAGCGTCGCCACTTCCTGAGAAGAGAGCATCTTATCGAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTC  >  NC_000913/3878289‑3878782

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