Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	77,776	C→A	G52G (GGC→GGA)	setA →	sugar exporter SetA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	77,776	0	C	A	98.9%	242.9 / NA	89	G52G (GGC→GGA)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (42/45);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCATGTGTGACTGAGTATTGGTG‑TAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGAT  >  NC_000913/77527‑78019                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    gCATCATGTGTGACTGAGTATTGGTG‑TAAAATCACCCGCCCGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTTTGAGTCGTGAGGGTGGCGCGCAACCTTTCTGGATCGGcc                                                                                                                                                                                                                                                    >  2:586892/1‑251 (MQ=255)     aatGTGGGACTGAGT‑TTGGTGGAAAAACCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCtt                                                                                                                                                                                                                                                <  2:497189/250‑1 (MQ=255)      atGTGTGACTGAGTATTGGTG‑TAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCttt                                                                                                                                                                                                                                               >  1:474185/1‑251 (MQ=255)        gtgtgACTGAGTATTGGTG‑TAAAATCACCCGACAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACTATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCttttt                                                                                                                                                                                                                                             >  1:841854/1‑250 (MQ=255)             actggttATTGGTG‑AAAAATCCCCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACgg                                                                                                                                                                                                                                        <  1:399170/245‑1 (MQ=255)             aCTGAGTATTGGTG‑TAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGAGCGCAACCTTTCTGGATCGGAATCTTTTCTACgg                                                                                                                                                                                                                                        >  2:591337/1‑251 (MQ=255)               agAGTATGGGTG‑AAAAACCACCCGCCAGCAGATTTTACCTGCTGGTTTTTTTTAATCTCGCCTCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAGCCGTTCGGGATCGGACTCTTTTATCCGGtt                                                                                                                                                                                                                                      <  2:150133/250‑2 (MQ=255)                   tATTGGTG‑TAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATg                                                                                                                                                                                                                                  >  1:114882/1‑251 (MQ=255)                    tttGGTG‑AAAAACCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATg                                                                                                                                                                                                                                  <  2:900009/249‑1 (MQ=255)                       ggTG‑TAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACCTTGAGCTTATTTCTGAGTCGTGCGGTTGGCGCGCCACCTTTCTGGACCGGACTCTTTTATACGGTGAATGCTAt                                                                                                                                                                                                                              >  2:673678/1‑251 (MQ=255)                        ggttaaaaaTCCCCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTAGGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTa                                                                                                                                                                                                                               <  1:499142/245‑1 (MQ=255)                             aaaaTCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTAGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGCTGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGGGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTg                                                                                                                                                                                                                         >  2:466461/1‑251 (MQ=255)                                  catccGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATc                                                                                                                                                                                                                    <  1:703907/248‑1 (MQ=255)                                     ccGCCAGCAGATTATTCCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTTAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGc                                                                                                                                                                                                                 <  2:48720/251‑1 (MQ=255)                                      cGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCg                                                                                                                                                                                                                >  1:441629/1‑251 (MQ=255)                                          agcagATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAACAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGCACGGAGTTTACGCGGCATTTATGCAGGTCGCTTTTATGATGTGGGTGGCCTGGGAGCTCCAGGCTCC‑TACATTGAGATTATTTCTGAGTCTTGAGGTTGCCGCGAAACCTTTCTGGTTCGGACTCTTTTATACGGTGATTGATCTTGCTGGGATCGTCGTagg                                                                                                                                                                                                            >  2:381061/1‑249 (MQ=255)                                               aTTATACCTGCTGGTTTTTTTTATTCTCGACGCGCTAAAACGGGACCGTATGATCTGGATAATGCCGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGATTTTATGTTGGGGGTGGCCGGGGAGCTACAGGCTCC‑TCCATTGCGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGACTGCTATTGCTGGGATCGGCGTAAGCctct                                                                                                                                                                                                       >  2:488250/1‑251 (MQ=255)                                                ttATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGAATGATCTGGATAATGACGTTGGCTCGCCGAATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGCTCGGACTCTTTTATCAGGTGATTGCTATTGCGGGGATCGGCGTAAGCCTCTg                                                                                                                                                                                                      >  2:760213/1‑251 (MQ=255)                                                 tataCCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTg                                                                                                                                                                                                      <  1:673678/250‑1 (MQ=255)                                                  ataCCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGGGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGCGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATGCGGTGAATGCTATTTCTGGGATCGGCGTAAGCCTCTGGt                                                                                                                                                                                                    >  2:286252/1‑251 (MQ=255)                                                   taCCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGtt                                                                                                                                                                                                   >  1:177514/1‑251 (MQ=255)                                                    aCCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGtt                                                                                                                                                                                                   <  1:688522/250‑1 (MQ=255)                                                        gctgGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGGGCATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCa                                                                                                                                                                                               >  2:695164/1‑250 (MQ=255)                                                        gctgGTTTTTTTTATTCTCGCCGCGATAAAAAGGGAACGTATGATCTGGATAATGAAGATGGCTCGCCGTATGAACGGGGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTgaca                                                                                                                                                                                               >  1:91844/1‑247 (MQ=255)                                                                   ttATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGCACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGACGTGAGGTTGGCGCGCAACCTTTCTGGATCGGCATCTTTTATACGGAGAATGCTGTTGATGGGGTCGGCGTAAGCCTCTGGTTGGCAACAAGTTCTGAc                                                                                                                                                                                   >  2:627499/1‑251 (MQ=255)                                                                               ggcTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGTATCCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCGCTGGTTGGCAAAGCGTTCTGACAGTCAGGGCGAt                                                                                                                                                                       <  2:135965/250‑1 (MQ=255)                                                                               cgcTAAAAAGGGAACGTATGATCTGGATAATGACGATGTCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCGCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGAt                                                                                                                                                                       <  1:274181/251‑1 (MQ=255)                                                                                gcTAAAAAGGGTACGTACGATCTGGATAAAGACGATGGCTCGCCGTAGGAACGGTGTTTACGCGGGATTTATGCTGGTCGCTGTTAGGATGGGGGGGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGGTGGCGCGCAACCTTTCTGGATCGGACTCTTTTAGACGGTGAATGCGATTGCTGGGAACGGCGTAAGCCTCTGGTTGGCAAATCGTTCTGACAGTCAGGGCGATc                                                                                                                                                                      <  2:83710/251‑1 (MQ=255)                                                                                gcTAAAAAGGGAACGTAGGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATc                                                                                                                                                                      <  2:829764/251‑1 (MQ=255)                                                                                   aaaaaGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCgg                                                                                                                                                                    <  1:143504/250‑1 (MQ=255)                                                                                     aaaGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGa                                                                                                                                                                 <  2:272753/251‑1 (MQ=255)                                                                                       aGGGAACGTATGATCTGGATAATTAAGATGCCTCGCCGTATAAACGGTGTTTACGCGGCATTTATGTTGGTCGCTTTTATGATGGGGGCGGCCGGGGCGCTACAGGCTCA‑TACATTGAGCTTTTTTCTGAGGCGCGAGGTTGACGCGCAACCTTTCTGAATCGTACTATTTTATACGGGAAATACTATAGCTGTGCTGGGCGTACGCCTATGGCTGGCAAAACGTTCTGACAGTCAGGGCGATCgacggac                                                                                                                                                               >  2:797881/1‑245 (MQ=255)                                                                                       aGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGCGATTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAGCCTTTCTGGTTCGGACTCTTTTATCCGGTGAATGCTATGGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGCCAGTCAGGGCGATCGGCGaa                                                                                                                                                                >  2:406168/1‑250 (MQ=255)                                                                                       aGGGAACGTATGATCTGGATAAGGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGTTCGGCGaaa                                                                                                                                                               <  1:846437/251‑1 (MQ=255)                                                                                            aCGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCGTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGTCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGa                                                                                                                                                          <  1:555848/251‑1 (MQ=255)                                                                                                       ggATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGCTCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAATCAGGGAGAGCGGCGAAAACTGTTTATATTTTGc                                                                                                                                               >  2:6048/1‑251 (MQ=255)                                                                                                       ggATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATAtttt                                                                                                                                                 <  1:681129/249‑1 (MQ=255)                                                                                                        gATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCt                                                                                                                                              <  1:695164/251‑1 (MQ=255)                                                                                                          tAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGt                                                                                                                                            >  1:612495/1‑251 (MQ=255)                                                                                                           aaTGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGtt                                                                                                                                           <  2:177514/251‑1 (MQ=255)                                                                                                           aaTGACGATGGCACGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGGCGCTTTTATGATGGGGGTGGCCGGGTCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTGCTGACAGGCAGGGCGATCGGCGAAAACTGATGATATTTTGCTGt                                                                                                                                            <  2:490883/250‑1 (MQ=255)                                                                                                            aTGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATTTTTTGCTgctt                                                                                                                                          >  2:368802/1‑248 (MQ=255)                                                                                                            aTGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGttt                                                                                                                                          <  1:846555/251‑1 (MQ=255)                                                                                                               aCGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGAt                                                                                                                                       >  1:561012/1‑251 (MQ=255)                                                                                                                 gATGGTTCGCCGTATGAACGGCGTTTACACGGCATTTATGCTGGTCGCTTTTATGATGGTGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTAAGGTTGGCGCCCCACCTTTCTGAATCGGAATCTTTTATACGGTGAATGCTATTGCTGGTGTCGGCGTACGCATCTGGTTGGCAAAACGTTCAGACAGTCCGGGCGATCGGCGAAAAATGATTCTATTTTGTTTTTTGATgg                                                                                                                                     >  2:537430/1‑251 (MQ=255)                                                                                                                 gATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTTC‑TTCATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATgg                                                                                                                                     <  2:452126/251‑1 (MQ=255)                                                                                                                           agTATGAACGGAGTTTACGAGGCCTTTATGCTGGCCGCATTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGAATCGGCGTAAGCCTCTGGTGGGCAAAACGGTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCggac                                                                                                                            >  2:385448/2‑248 (MQ=255)                                                                                                                            ggaTGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAAt                                                                                                                          <  1:677557/249‑1 (MQ=255)                                                                                                                            gTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAAt                                                                                                                          >  1:66086/1‑251 (MQ=255)                                                                                                                              aggAACGGTGGTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTAGACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATgt                                                                                                                        <  2:13745/249‑2 (MQ=255)                                                                                                                              aTGAACGGTGTTTACGAGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCTGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGATTCGTGAGGTTGGCGCGCAACCTTTCTGGTTCGGACTATTTTATGCGTTGATTGCTATTGCTGTGTTCGACGTAAGCCTCTGGTTGGAAAAACGTTCTGTTCGTCAGGGCGATCGGCGAAAGCTGATTATATTTTGCTGATTTATGGCTGTCGGCAATgg                                                                                                                        >  2:698718/1‑250 (MQ=255)                                                                                                                               tGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGGGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATgcg                                                                                                                       >  1:429566/1‑251 (MQ=255)                                                                                                                                       gtTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGCGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGTCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGtt                                                                                                               >  2:279515/1‑251 (MQ=255)                                                                                                                                       gtTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGtt                                                                                                               <  1:877070/251‑1 (MQ=255)                                                                                                                                       ggttACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGtt                                                                                                               <  1:152672/249‑1 (MQ=255)                                                                                                                                                     tATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGt                                                                                                  <  1:470407/250‑1 (MQ=255)                                                                                                                                                     tATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTc                                                                                                 >  1:726002/1‑251 (MQ=255)                                                                                                                                                          tGGTCGATTTTATGATGGGGTTTGACGGGGCGCTGCAGGCTCC‑ACCATCAAGCTTATTTCTGAGTCGTGAGGGTGGCGCGCAACATTTCTGGATCGGACACTTTTATACGGTGACTGCTATTGCTGGGATCGGCGTAAGCCTCTGGTGGGCAAAACGTTCTGACAGTCAGGGCGTTCGGCGAAATCTGTTTATATTTTGCTGTTTGATTGCTATCGGCAATGCGCTATTGTTTGCATTTAAGTGTCATTAt                                                                                            >  2:445877/1‑251 (MQ=255)                                                                                                                                                            gTCGCGTTTATGATGGGGGTGGCCGTGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGGTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTGTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCt                                                                                          <  1:408834/251‑1 (MQ=255)                                                                                                                                                              agattttATGAGGGTGGTGGTCGGGGCGCTACCTG‑TCCATACATTGAGTTTGTTTCTTAGTCGTGAGGGTGGCGCGCAACCTTTCGGGATCTGACTCTTGTATACGGTGAATGCGATTGCTGGGATCGGCGGAAGCCTCTGGTTGGCAAACTGGTCTGACATTTAGAGCGATCGGTGAAAACTGATTATATTTTGCTGTTTGATGGTTACCAGCAATGCGCTATTGTTTTCATTTTATAGTCATTAtatg                                                                                         <  2:605484/247‑4 (MQ=255)                                                                                                                                                                    tATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGTTCGGACTCTTTTATACGGTGAATGATATTGCTGGGACCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGGCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGGTTGCATTTAATCGTTATTACCTGAAGCtt                                                                                   <  2:570551/250‑1 (MQ=255)                                                                                                                                                                    tATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCATCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCtt                                                                                   >  2:46155/1‑250 (MQ=255)                                                                                                                                                                     atgatgGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTAt                                                                                 <  1:767325/251‑1 (MQ=255)                                                                                                                                                                      tgatgGGGGTAGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCCACATTTCTGGATCGGACTCTTTTATCCGGTGACTGCTATTGCTGGGATCGGCGTAAGCCACTGGTTGGCAACACGTTCTGACAGTCAGGGCGATCAGCGAACACTGATTATATTTTGCTGTTTGATTGCTATCGGCCATGCGCTATTGTTTACATTTCATAGTCATTATCTGTCGCttctc                                                                                >  2:972549/1‑248 (MQ=255)                                                                                                                                                                          ggggggagcaggggCGCTATAGGCTCC‑GGCATTGAGCTTATGTCTGAGTCGTTAGGTTGGCGCGCATGCTTGCTTGATCGGACTCTTTTATACGGTTAATGCTATTGCTTGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCTATCGGCGAAAACTTATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTGAATCGTCATTATCTGTCGTTTATCACCt                                                                            <  2:661724/241‑1 (MQ=255)                                                                                                                                                                          gggggTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCAcc                                                                             <  1:475975/250‑1 (MQ=255)                                                                                                                                                                          gggggTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCAcc                                                                             <  1:711168/250‑1 (MQ=255)                                                                                                                                                                               tGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTgtggg                                                                       >  1:968320/1‑250 (MQ=255)                                                                                                                                                                                   cGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGc                                                                   >  1:192968/1‑251 (MQ=255)                                                                                                                                                                                     gggCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCtt                                                                 <  1:884860/251‑1 (MQ=255)                                                                                                                                                                                           tACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAGAACGTTCTGACAGTCAGGGCGGTCGGCGGGAACTGATTATATTTTGCTGTTTGATGGATATGGGCAATGCGCTATTGTTTGCATTTAATCGACATTATCTGACGCTTATCACCTGTGGTGTGCGTCTgac                                                            >  1:39368/1‑248 (MQ=255)                                                                                                                                                                                            aCAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGc                                                            <  1:976977/249‑1 (MQ=255)                                                                                                                                                                                            aCAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCGTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGc                                                            <  2:474185/249‑1 (MQ=255)                                                                                                                                                                                            aCAGGCTCC‑TACAGTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGAGCTGACTCTTTTAGACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGTCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCTCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCAt                                                          <  1:980304/251‑1 (MQ=255)                                                                                                                                                                                              agactac‑TACATTGAGCATATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATCCGGTGACTGCTATTGCTGGGCTCGGCGTAAGACTCTGGTTGGCAAAACGTTCTGACAGTCAGGGAGATCGGCGCAAACTGATTATATTTTGCTGTTTGATGGCGATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGCCGCTTATCACAGGTGGTGTGCTTCTGGCctct                                                        >  2:751393/7‑251 (MQ=255)                                                                                                                                                                                               ggCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCAtctc                                                       >  1:11226/1‑251 (MQ=255)                                                                                                                                                                                                        cATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAAt                                               >  1:364508/1‑251 (MQ=255)                                                                                                                                                                                                         aTTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCTTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTTCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATa                                              <  1:753145/251‑1 (MQ=255)                                                                                                                                                                                                           tGAGCTTATTTCTGAGTCGTGAGGTTGGCGCTCAACCTTTCTGGATCGGACTCTTTTATACGGTGAAGGCTATTGCTGGGAGCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCTGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATTGCTTTCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTTACGCTTATCACCTGTGGTGTGCTGCTGGCATCTCTGGCCAATAc                                             <  2:595182/250‑1 (MQ=255)                                                                                                                                                                                                            gAGCTTATTTCGGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTGTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTAGCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACg                                            <  1:825321/250‑1 (MQ=255)                                                                                                                                                                                                                   tttCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGAGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGcc                                    >  2:734956/1‑251 (MQ=255)                                                                                                                                                                                                                     tCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCCATGCcac                                  >  1:575355/1‑251 (MQ=255)                                                                                                                                                                                                                         aGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGtt                              >  1:886602/1‑251 (MQ=255)                                                                                                                                                                                                                            cGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTAt                            <  1:751393/250‑1 (MQ=255)                                                                                                                                                                                                                                       cgcgcAACCTTTCTGGATCGTACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGgcgc                >  1:469454/1‑251 (MQ=255)                                                                                                                                                                                                                                         cgcAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCgg              >  2:760853/1‑251 (MQ=255)                                                                                                                                                                                                                                            aaCCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGCTCGGAGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCACTGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACATGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTCTTTGCTCTGGCGCGGGa            >  2:443596/1‑250 (MQ=255)                                                                                                                                                                                                                                             aCCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTTTTGCTGTGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGGCATTATCTGACGCTTATCACCTGTGGTGTGCTTCGGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGAGCGGGAat          <  2:349162/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ttCTGGATCGTACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCACCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCGGGCGCGGGAatat        <  2:114882/249‑1 (MQ=255)                                                                                                                                                                                                                                                      gATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGAt  >  1:205287/1‑250 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    GCATCATGTGTGACTGAGTATTGGTG‑TAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCC‑TACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGAT  >  NC_000913/77527‑78019

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