Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,847	C→T	G497G (GGG→GGA)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,847	0	C	T	98.8%	268.2 / NA	86	G497G (GGG→GGA)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (55/30);  total (55/31)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑AGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGA  >  NC_000913/66559‑67135                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                 cTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTNTCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATacac                                                                                                                                                                                                                                                                                       <  1:619118/300‑1 (MQ=255) cTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCATACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATacac                                                                                                                                                                                                                                                                                       >  1:58227/1‑300 (MQ=25)              gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgc                                                                                                                                                                                                                                                                          <  1:182033/300‑1 (MQ=255)                   ttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGAATCAGGTGTTACACACCTGATGCGACTAATTAGCGACGAAATCCGTAATACACTTCGTTCCagcgcagcgc                                                                                                                                                                                                                                                                     >  2:802103/1‑300 (MQ=25)                        tATCCGGCCTACAAAATCGCACCGTGGAGTCCTGATAAAACGCGCCAACGGCACATCAGGCGTTGAATACCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGAGAAGACGCGCCAGCGGCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGTGTGTAGGCC‑TGATAAGACGCGTCAGCTTCGCATCAGGCGCTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCtt                                                                                                                                                                                                                                                                <  2:462969/300‑1 (MQ=18)                              gCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGAAGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGAAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGACGCGACTACTTAGCGACGAAGTCCGTAATACAATTCGTTCCAGCGCAGCGCGGCTTTAAACg                                                                                                                                                                                                                                                          >  2:771856/1‑300 (MQ=21)                                               gtgtAGGCCTGATAAGACGCGCCAGCGTCGCCTCAGGCGTTGACTGCCGGATGCGCTTTGCTTATCCGGCATACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTCCACACCTGATGCGACTACTTAGCGACGAAGTCCGTAATACACTTCGTTCAAGCGCAGCGTGTCTTTAAACGCTggcgggcgggtatcg                                                                                                                                                                                                                                         >  2:339047/1‑293 (MQ=18)                                               gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAACATCGGAACGCGTAGGCCTGATACGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTAGTTAGCGACGAAATCCGTAATACACTTAGTTCCAGCGCGGCGCGGCTTTAAACGCTGGCAGGCgggtatcg                                                                                                                                                                                                                                         >  2:85826/1‑293 (MQ=21)                                               gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGCNGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGGCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTATTTAACCGCTGGTAGGcaggtatcg                                                                                                                                                                                                                                         >  2:204876/1‑292 (MQ=21)                                               gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAACCGCTGGCAGGCgggtatcg                                                                                                                                                                                                                                         >  1:667592/1‑293 (MQ=21)                                               gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgggtatcg                                                                                                                                                                                                                                         >  1:133003/1‑293 (MQ=25)                                               gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTAAACATCTGATGCGACTACTTCGCGACGAAGTCCGTACTACACTTCGTTCCAGCGCAGCGCGTATTTAAACGCTGGCAggagggtagcg                                                                                                                                                                                                                                         >  2:514770/1‑291 (MQ=18)                                               gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGCAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTCGCGACGAAATCCGGAATACACTTCATTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCgggtatcg                                                                                                                                                                                                                                         >  1:782183/1‑293 (MQ=21)                                               gtgtAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCAGGATGCGCTTTGCTTATCCGGCCTACAACATCGCAGCGTGTAGGAC‑TGATAAGACGCGTCAGCTTCGCATCAGGCGTTCCCCACGTGATGCGACTACTCAGCGACGAAATCCGTCATACACTTCGTTCACGCGCAGCGAGTCTTTAAACGATGGAAGGCGGGTGTCg                                                                                                                                                                                                                                         >  2:513529/1‑300 (MQ=18)                                                    ggCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatc                                                                                                                                                                                                                                    >  1:402925/1‑300 (MQ=25)                                                      ccTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatcaa                                                                                                                                                                                                                                  <  1:193619/300‑1 (MQ=25)                                                                 gcgcCAGCGTCGCATCAGGCTTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTNCAAAATCTGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGAtt                                                                                                                                                                                                                       <  2:640177/300‑1 (MQ=21)                                                                       gggTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTCGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTCTGCTTATCCGGCCTACAAAATCGAAGCGCGTACGCC‑TGATAAGACGCGTCAGAGTCGCATCAGGCGTTACACCCCTGATGCGAATACTTAGCGCAGCAATCCGTAATACACTTCGTTCCAACGCAGNGCGGCTTTCAACACTGGCAGGTGGGTATCGTTCGCACTCACCGTGGTNTCAATg                                                                                                                                                                                                                 >  2:708827/3‑300 (MQ=12)                                                                       gCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATg                                                                                                                                                                                                                 <  1:419422/300‑1 (MQ=25)                                                                            gCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACNCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTNTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTTAATGTCGTg                                                                                                                                                                                                            >  2:790891/1‑300 (MQ=21)                                                                               tCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCTGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCAt                                                                                                                                                                                                         <  1:112552/300‑1 (MQ=25)                                                                                     gTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGc                                                                                                                                                                                                   >  1:22979/1‑300 (MQ=25)                                                                                     gTTGAATGCCGGATGCACTTTGCTTATCCGGNCTACCAAATCAGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGACGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAANTCGCAGCGTGTAGGCC‑CGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACATGATGCGACTAATTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTCTCAATCAACGTGATTTCAATGTCGTGCCGCTCGGc                                                                                                                                                                                                   >  2:615721/1‑300 (MQ=21)                                                                                                   gcgcGTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGGCGCATCAGGCGTTGAATGCCGTATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCatat                                                                                                                                                                                     <  2:78100/300‑1 (MQ=25)                                                                                                    cgcTTTGCTTATCCGGCCTACCAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCATCGCATCAGGCATTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGAGTTACACGCCTGATGCGACTACTTAGCGACGACATCCGTAATACCCTTCGTTCCAGCGCAGCGCGTCTTTAAATGTTGGCAGGCGGGTACTGTTATCAATCACAATGATTCCAATGTAGTGCATCTAGGCGAACCGGCGCGTATc                                                                                                                                                                                    >  2:20987/1‑300 (MQ=17)                                                                                                    cgcTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATc                                                                                                                                                                                    >  1:780532/1‑300 (MQ=25)                                                                                                        ttGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTg                                                                                                                                                                                <  2:415036/300‑1 (MQ=25)                                                                                                           cTTATCCGGCCTACAAAATCTGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACGAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCCGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAgg                                                                                                                                                                             <  2:46412/300‑1 (MQ=25)                                                                                                                  ggCCTACAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcg                                                                                                                                                                      <  1:61371/300‑3 (MQ=25)                                                                                                                        caaaatcggaacatgTAGGCCCGATAAGACGCGCCAGCGTCACAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TAATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGc                                                                                                                                                                <  2:377662/286‑1 (MQ=25)                                                                                                                        cAAAATCGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGc                                                                                                                                                                <  1:538346/300‑1 (MQ=25)                                                                                                                             tcggaacacgTAGGCCTGACAAGACGCGCCGGTGTCGCATCAGGCGTTGAATGCCGGATACGCTTTGCTTATCCGGCCTACGAAATCGCAGCGTGTAGGTC‑TGATAAGACGCGGCAGCGTCGCATTAGGTGTTACACACCTGATGCGACTGCGTAGCGACGAGATTCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAATGCTGGTATGCGGGTATCGTTATCAATCACCGTGATTTCAATTTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGTTTCAGCGCATGGCTGAAg                                                                                                                                                           <  2:315585/292‑1 (MQ=255)                                                                                                                               ggaacgcgTAAGCCTGATAAGGCGCGCCAGCGTCACATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGNGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGAc                                                                                                                                                         <  1:623248/296‑1 (MQ=37)                                                                                                                                    gtgTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGGAACGCGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTAGGCGAACTGGCGAATATCGTTGAGGTTCAGCGCATGGTTGGAGACGGTAt                                                                                                                                                    >  2:351400/3‑300 (MQ=21)                                                                                                                                      atagtccTGATAAGACTCGCCAGCGTCGCATTAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATTTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATgg                                                                                                                                                  <  2:523717/295‑1 (MQ=37)                                                                                                                                        cgaccTGATGAGACGCGCCAGCGTCGCATCGGCCATCGAATGCCGGAGGCGCTGTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGTCC‑TGGTAAGACGCGTCGGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCTTAATACATTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGAGTAGCGGTATCAATCACCTTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTg                                                                                                                                                <  2:428756/297‑1 (MQ=255)                                                                                                                                         ggCCTGATAAGACGCGCCAGCGTCGCATCAGNCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACANACCTGATGCGACTACTTAACGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTGTCGTTATCAATCACCATGATTTAAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGATGAAGACGGTATGGTgc                                                                                                                                               >  2:603204/1‑300 (MQ=37)                                                                                                                                         ggCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCTGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTgc                                                                                                                                               <  1:365213/300‑1 (MQ=37)                                                                                                                                           ccTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTgcgc                                                                                                                                             >  1:712630/1‑300 (MQ=37)                                                                                                                                            cTGATAAGCCGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGTGTTACACACCTGCTGCGACCACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGCGATTTCAATGCCGTGCATCTCGGCGAACTGCCGCATATCGTTGCGGTTCAGCGCATGACTGAAGACGGTATGGTgcgca                                                                                                                                            >  2:142029/1‑299 (MQ=25)                                                                                                                                             tGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGAATAATTAGCGACGAAATCCGTCATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGTGCATATCGTTGAGGTTTAGCGCATGGCTGAAGACGGTATGGTgcgcgc                                                                                                                                           >  2:17786/1‑300 (MQ=25)                                                                                                                                                     cgcgCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCg                                                                                                                                   <  1:549513/300‑1 (MQ=37)                                                                                                                                                           agcttCGCATCAGGCGCGGGATGCTGGATACGCTTTGCTTGTCCGGCCTACCAAATCGCAGCGTGTAGGAC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGCTACACACCTGATGCGACTACTTAGCGACGAGACCCGTAATACACTCCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGGTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATc                                                                                                                             <  2:419662/296‑1 (MQ=255)                                                                                                                                                                cGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGc                                                                                                                        <  1:398494/300‑1 (MQ=38)                                                                                                                                                                   gttgactccttcgtGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATAAACTTAGTTCCAGCGCAGCGCGTCTTTACACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCCGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTc                                                                                                                     >  1:559127/14‑300 (MQ=37)                                                                                                                                                                        gCGTTGAATCCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGTGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAg                                                                                                                <  2:141280/300‑1 (MQ=38)                                                                                                                                                                                         cgcTTTGCTTATCCTGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGt                                                                                               <  2:412739/300‑1 (MQ=255)                                                                                                                                                                                         cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGTGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGCAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCCGGCGAGTATCGTTATCAATCACCGTGGTTTCACTGTCGTGCATCTCGGCGAACTGGCGAATATCGTTGAGGTTCAGCGCATGGCTGACGATGGTATGGTGCGCGCCACCAGCGAGGATTCACGCTTCGGACGCAGTTGGTAGATCCGGt                                                                                               >  2:463092/1‑300 (MQ=255)                                                                                                                                                                                         cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATTAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCAGAAGCAGTTGGCAGATCCggc                                                                                               >  2:155979/1‑299 (MQ=38)                                                                                                                                                                                         cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGCTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCAACAGCGAGGATCCACGCTTCGGAAGCAGTTGTCAAATCCGGt                                                                                               >  2:468487/1‑300 (MQ=38)                                                                                                                                                                                                  tATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCAACGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTc                                                                                      >  1:707115/1‑300 (MQ=255)                                                                                                                                                                                                     ccGGCCTACAAAATCGCAGCGTGTAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCcac                                                                                   <  1:300646/300‑1 (MQ=255)                                                                                                                                                                                                        gtaggcagaaTCGCAGCTTGTGGGCT‑TGGTAAGACGAGTCAGCGGCGCAGCAGGCGTTACACACCTGAGGCGACTGCTTGGTGATGAATTCTGTGATACACTTCGTTGCAGCGCAGCGCGTCTTGAAACGCGGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCTGCGAACGGTCACATATCGGTGAGGCTCAGCGCGTGGCTGAGGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTGCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAgc                                                                                <  2:436505/292‑1 (MQ=255)                                                                                                                                                                                                        gCCTACAAAATCGCAGCGTGTAGTCCTAG‑TAAGACGCGTCAGCGTCGCAGCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAgc                                                                                <  2:163044/300‑1 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGCCGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGNTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTNAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCCCAGCGCATTCGCCACCGGCAg                                                               >  2:684735/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCNTCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTCATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGNTATCGTTATCAATCACCGTGATTTCAATGTCGTGCCTCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCAACAGCGAGGATCCACGATTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCAACGGCAg                                                               >  2:591652/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCCTCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAACCGCTGGCAGGCGGGTATCGTTATCACTCACCGTGCTTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGTGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATTCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                                                               >  2:460327/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATGTCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGTGCGCCAACAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGGTCCGGTTGCGCTTTCCACCGCGCATTCGCCACcagcag                                                               >  2:379405/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCCCAGCGCATTCGCCACCGGCAg                                                               >  1:638473/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                                                               >  1:225062/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                                                               >  2:504611/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTAGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                                                               >  1:638435/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTAGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGCATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCAACCGGCAg                                                               >  2:108240/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGCATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATAGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTAGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                                                               >  2:415678/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCATAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGNTATCGTTATCAATCACCGTGATTTCAGTGTCGTGCATCTCGGCGAACTGGCGCATATCGTTNAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGGAGTTGGCAGATCCGGTTGTGCTTTCCACAGCGCATTCGCCACCGGCAg                                                               >  2:690746/1‑300 (MQ=255)                                                                                                                                                                                                                         gtgtAGGCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGACACGANTACCTAGCGACAAAATCCGTAATCCACTTCGTTCCAGCACAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGAGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGCAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTNTGAAGCAGTTGGCAGATTCGGTTGCGCTTTCCACAGCGCATTCGCCaacgatag                                                               >  2:791166/1‑293 (MQ=255)                                                                                                                                                                                                                              ggCC‑TGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGATGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAATTTCg                                                          >  2:400011/1‑300 (MQ=255)                                                                                                                                                                                                                                gcctgATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCgg                                                         >  1:703873/5‑300 (MQ=255)                                                                                                                                                                                                                                 cctgATAAGACGCGTCAGCTTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAAAGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGc                                                        <  1:784396/297‑1 (MQ=255)                                                                                                                                                                                                                                    gATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACTCGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGGTCCGGTTGTGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAgg                                                     >  2:307936/1‑300 (MQ=255)                                                                                                                                                                                                                                            gcgTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCgg                                             <  1:402713/300‑1 (MQ=255)                                                                                                                                                                                                                                             cgTCAGCGTCACATCAGGCGTTACACACCTGATGCGACTACTTATCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGgt                                            <  2:273608/300‑1 (MQ=255)                                                                                                                                                                                                                                                  gCGTCGCATCAGGCGTTCCACCCCTTATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTCAACGCCGGCAGGCGGGTATCGTTATCAACTACCGTGACTTCAATGTCGTGCATCTCAGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATTGTGCGCGCCACCAGCGAGGATCCACACTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGACAACGGTAGTTCCGGAAGGGGGTGCGGTGtttt                                       >  2:80763/1‑300 (MQ=255)                                                                                                                                                                                                                                                        cATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCAACGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGCCGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTAGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCCTTCGCCCACGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCgt                                 >  2:138183/1‑300 (MQ=255)                                                                                                                                                                                                                                                                gTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGAtac                         >  1:211090/1‑298 (MQ=255)                                                                                                                                                                                                                                                                gTTACACACCTGATGCGACTACTTAACGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTAGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGCAGCAGTTGGCAGATCCGGTTGCGCTCTCCCCCGCGCATTCGCCCCCGGCAGTTTCGGCAGGGAGTGCGGTGCTTTCACCGTGTCGAtac                         >  2:457213/1‑298 (MQ=255)                                                                                                                                                                                                                                                                 ttACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCCGCGCAGCGCGTCTTTAAACGCTAGCAGGCGGGTATCGTTATCAATCACCGTGATTTCACTGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCTGGTTGCGCTTTCCACAGCGCATTCGCCACCGGTAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGAtaca                        >  2:242560/1‑297 (MQ=255)                                                                                                                                                                                                                                                                 ttACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCCGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTAGAtaca                        >  2:315871/1‑297 (MQ=255)                                                                                                                                                                                                                                                                 ttACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGAtaca                        >  1:308356/1‑297 (MQ=255)                                                                                                                                                                                                                                                                 ttACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTAGTGCATCTAGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTAGGAAGCAGTTGGCAGATCAGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTATCGAtaca                        >  2:285587/1‑297 (MQ=255)                                                                                                                                                                                                                                                                  tACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAg                       <  1:256912/300‑1 (MQ=255)                                                                                                                                                                                                                                                                      cacctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTaa                   >  1:809023/6‑300 (MQ=255)                                                                                                                                                                                                                                                                      cacctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTAGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGNATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTaa                   >  1:604874/6‑300 (MQ=255)                                                                                                                                                                                                                                                                             tGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAgcag            >  2:402417/1‑297 (MQ=255)                                                                                                                                                                                                                                                                                    aCTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGCGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCNCCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCAACGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGAcaataa     >  2:200343/1‑295 (MQ=255)                                                                                                                                                                                                                                                                                       tAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGa  <  2:430830/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                 CTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCC‑AGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAAGCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGA  >  NC_000913/66559‑67135

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