Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	161,653	0	T	G	26.3%	159.2 / 11.0	76	T127P (ACC→CCC)	thpR	RNA 2',3'‑cyclic phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (20/36);  new base G (20/0);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.33e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATCGGGTGTGATCACATCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAG  >  NC_000913/161405‑161898                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      cATCGGGTGTGATCACATCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGGGGCGTAAAGGGGGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑GATGTCACCGCCTCGCCGGCGTCGCGCCATAAGGTa                                                                                                                                                                                                                                                      >  1:261852/1‑250 (MQ=255)           gATCACATCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCCATAAGGGAATAtgggggtg                                                                                                                                                                                                                                           >  1:151312/1‑245 (MQ=255)           gATCACATCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGAt                                                                                                                                                                                                                                            <  1:204966/250‑1 (MQ=255)               acaTCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCCAACGAGGAGGCGGAAAGGGGGAACTCCGTCACCGCATAC‑GAAC‑AGTTAAAAACTGGCGGCGGG‑GTTGTCCCCCCCCCCCTTGCCTCCCGCAAAAAGGGAATTTGTGGGTGaaa                                                                                                                                                                                                                                        >  2:215027/1‑250 (MQ=255)                caTCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑GATGTCACCGCCTCGCTGGCGTCGCGCCATAAGGGTATATGTGGATGAAACg                                                                                                                                                                                                                                      >  1:382452/1‑251 (MQ=255)                caTCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGGGTGCGTCCACGGGCAAACGGGGAGGCGTAAAAGGTGAACTCCCTCACCGCATAC‑GACC‑AGTTAAAACCTTGCGGCGGG‑ATTTTCACCCCCTCCCTTGCGTCGCGCCATAAAGGAATATGTGGATGAAACg                                                                                                                                                                                                                                      >  1:100944/1‑251 (MQ=255)                 aTCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑GTTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGGAATATGTGGGTGaaaagg                                                                                                                                                                                                                                     >  1:370399/1‑248 (MQ=255)                         aaaaaCGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGCGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑GTTGTCACCGCCTCGCTGGCGTCGCGCACTAAAGTAATATGTGGGTGAAACGGACGCTTGc                                                                                                                                                                                                                             >  1:177085/1‑251 (MQ=255)                                   gTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCCATAAGGGAATATGTGGATGAAACGGACGATTGCTTTTAAaac                                                                                                                                                                                                                    >  1:468060/1‑250 (MQ=255)                                      ccGCTGAATTAGCGTCGCGCTGTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTTGCGGCGTGT‑GGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCGCCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAAcc                                                                                                                                                                                                                <  2:528098/250‑1 (MQ=255)                                        gCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCCCGGGCAAAAGAGGAGGCGTCAAGGGGGAACTCCGTCACCCCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑GTTGTCCCCGCCCCGCTGGCGCCGCGCCAAAAGGGAACATGTGGCTGGAACCGACGctttctttttaaactaacgc                                                                                                                                                                                                              >  2:273188/1‑239 (MQ=255)                                        gCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGGCCACGGGCAAACGAGGAGGCGTAAAGGGGGAACTCCGTCACCGCATAC‑GACCCGGTTAAAAC‑TGGCGGCGGG‑ATTGTTACCGCCTCGCTGGCGTCGCGCAATCAGGGAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGc                                                                                                                                                                                                              >  1:341199/1‑251 (MQ=255)                                          tGAATTAGCGCCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGGGTGCGTCCACGGGCAAACGAGGAGGCGGAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑CGTTAAAACCTGGCGGCGGG‑GTTGTCACCGCCTCGCTGGCGGCGCGCAATAAGGGAATATGGGGGTGAAACCGACGGTTGCTTTGTAAACAACCGCTg                                                                                                                                                                                                            >  1:353307/1‑251 (MQ=255)                                           gAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATGGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATCC‑GTCC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCt                                                                                                                                                                                                             <  1:7695/249‑1 (MQ=255)                                               tAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGGGTGCGTCCACGGGCAAACGAGGAGGCGGAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTGATATGTGGATGAAACGGACGAATTCTTTTTAAAAAACCGCTGCgggg                                                                                                                                                                                                       >  2:267654/1‑250 (MQ=255)                                                aGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑GTTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGGTTGCTTTTAAAAAAACCGCTGCGGGcc                                                                                                                                                                                                      >  1:261737/1‑250 (MQ=255)                                                      gcgcgcTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCtgt                                                                                                                                                                                                <  1:387620/251‑1 (MQ=255)                                                          gcTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGG‑AGCGTGTGCGTCCTCGGGAAAAAGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTTTCCACGCCTCGCTGGGGTCGCGCCATAAGGGAATATGTGGGTGAAACGGACGATTTCTTTGTAAAAAACCGCTGCGGGCAGCCcgtgcac                                                                                                                                                                                            >  2:227608/1‑248 (MQ=255)                                                               agagGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCGGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAg                                                                                                                                                                                        <  1:95239/250‑1 (MQ=255)                                                                 agGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTCGCGGCGTAT‑AGACTGTGAGACCACGGGGAAAAGAGGAGGTGTAAAGGGTGAAATCCGTCACCGCATAC‑CACC‑CGTTAAAACCTGGCGGCGGG‑GTTGTCACCGCCTCGCCGGCGTCGCGCGATAAGGGAATATGTGGATGAAAAGGACCCTTGCTTTTAAAACAACAGCTGCGGGCAGCCTGTGCACAGAGCat                                                                                                                                                                                     >  1:449306/1‑251 (MQ=255)                                                                     ggAGAAAATTCCATTGCGACAATCCTTGTTATTGCTTTAGCGCACAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTc                                                                                                                                                                                 <  2:516228/251‑1 (MQ=255)                                                                     ggAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCCCGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGGTTAAACCCGGCGGGGGG‑ATTTTCACCGGCTTGCTGGCCTGCCGCAATAAGGGAATAATCGGGTGAAAAGGACGATTGCTTTGTAAACAACCGCAGCAGGCAGCCCGTGAAAGGAGCacactt                                                                                                                                                                                 >  2:96303/1‑246 (MQ=255)                                                                          aaaTTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGGGGCGTAAAGGGGGGACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑GTTGTCAACCCCTCGCTGGCGGCGCGCAATAAGGGAATATGTGGGAGAAAAGGGCGGATGCTTTGAAAACAACCCCTGCGGGCGACCTGTTAACggggccaatttcgca                                                                                                                                                                             >  2:305771/1‑237 (MQ=255)                                                                            aTTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTCCCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAAAGAGGAGGCGTCAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGATAAAACCTGGCGGCGGGGATT‑TCCCCGCCCCGCTGGCGTCCCGCCATAAGGTAATATGTGGATGAAACGGGCAGTTGCTTTGAAAAAAACCGCCGCGGGGAGCCCGTGAACGGAGCATATTCGcc                                                                                                                                                                              >  1:153817/1‑247 (MQ=255)                                                                              tCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAA‑GACC‑CATTAAAACCTGGCGGCGGG‑AGTGTCCCCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAAACGCCGCGGGCCGCCTGTGCACGGAGCATATTAGCCAGCTgg                                                                                                                                                                        >  1:228259/1‑251 (MQ=255)                                                                                     cGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGc                                                                                                                                                                 <  1:267654/251‑1 (MQ=255)                                                                                      gACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAAAGAGGAGGCGTAAAGGGTGGACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTAACCGCCCCGCCGGCGTCGGGCAATAAGGGAATATGGGGGTGAAACGGACGATTGATTTTAAAAAAACCGCTGCGGGCCGCCTGTGAACGGGGCATATTCGCCAGATGGATTAAGcc                                                                                                                                                                >  1:543734/1‑251 (MQ=255)                                                                                        cAATCCTTGGTTTTTCGGTAGCCCCCCGCGTTTTAGCGGCG‑GTGAGCGTGTGCGTCCACGGGCAAACGGGGAGGCGTAAAGGGGGAACCCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCcgc                                                                                                                                                              <  2:150951/251‑1 (MQ=255)                                                                                         aaTCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAAGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCCGGCGTCGCGCAATAAGGGAATATGTGGATGAAACGGACGGTTGCTTTGAAAACAACCCCTGCGGGCAGCCTGTGGACGGAGGATAATCGCCAGCTGGATTAAGCcgcg                                                                                                                                                             >  1:540199/1‑251 (MQ=255)                                                                                          aTCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCCCGGGCCAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑GTTGTCACCCCCCCGCCGGCGCCGCGCACTAAGGGAATATGTGGGTGAAAAGGACGATTGCTTTGAAAACAAACGCCGCGGGCAGCCCGTGAACGGAGCATATTTCCCCGCTGGATTaaaccccgg                                                                                                                                                            >  2:496560/1‑244 (MQ=255)                                                                                                gTTATTGCGTTAGCTCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTAAACTCCTTCACCGCATAC‑GACC‑AGTTAAACCCTGGAGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCTCGTggcgg                                                                                                                                                       <  2:310466/250‑1 (MQ=255)                                                                                                      gCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACg                                                                                                                                                 <  1:481155/250‑1 (MQ=255)                                                                                                           aGCGCCCGGGGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACGCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATccc                                                                                                                                           <  2:198346/251‑1 (MQ=255)                                                                                                                ccAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAacc                                                                                                                                      <  1:338102/251‑1 (MQ=255)                                                                                                                 cAGCGTTTTAGCGGCGTGT‑AGAGAGTGCGTCCACGGGCAAACGAGGAGGCGGAAAAGGTGGACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑AATGTTACCGCCTCGCTGGCGTCGGGCAATAAGGGAATATGGGGCGGAAACGGAAGATTGCTTTGGAAACAACCCCCGCGGGGAGCCTGGGAACGGGGCATATTCGCCCGCTGGATTAAGCCGCGGggcggggcacgcctcccccaccc                                                                                                                                     >  2:518459/1‑233 (MQ=255)                                                                                                                  aGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGGAATATGTGGGTGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGGGGCGGCTGACGCATCCCTAaccac                                                                                                                                    >  1:570632/1‑251 (MQ=255)                                                                                                                        tttaggggcggctagcttctgCGTCCCCGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCCCCTCATAC‑GACC‑AGTTAAAACCTGGGGGGGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGTTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAaccacacca                                                                                                                                <  2:324604/231‑1 (MQ=255)                                                                                                                            cggcgTGT‑AGCGGGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGGGAACTCCGTCAACGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGACTCGCTGGCGGCGCGCAATAAGGTAATATGTGGGTGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGAATATTCGCCAGCTGGATTAAGCCGCGTGGCCGCGTACGCATCCCTAAACACCCCACAcgcg                                                                                                                          >  2:568963/1‑251 (MQ=255)                                                                                                                                gtgt‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGGGGCGGCTGACGCATCCCTAACCACACCACACGCGAACg                                                                                                                      >  1:548877/1‑251 (MQ=255)                                                                                                                                    gaGCGTGGGCGGCCACGGGCAAACGTGGAGGCGTAAAGGGTGAACCCCGTCCCCGCATAC‑GCCC‑AGTTAAACCCTGGCGGGGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCTCACTCGAACGCAg                                                                                                                   <  2:443554/250‑1 (MQ=255)                                                                                                                                        gtgtGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACCCAGCCAt                                                                                                               >  2:421521/1‑251 (MQ=255)                                                                                                                                          gggggtcgactggcaatagaggagGCGTAACGGGTGAGCCCCGCCACCGCATAA‑GCCC‑GGTTAAAGCCTGGCGGGGGG‑TTTCACACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCCTATTCGCCAGCTGGATTAAGCCGCGTGGCGGATGACGCATCCCTAACCACACCACACGCGAACGCAGCTAtt                                                                                                              <  2:261659/232‑1 (MQ=255)                                                                                                                                            gCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGt                                                                                                            <  1:10921/250‑1 (MQ=255)                                                                                                                                              gTCCCCGGTCAAACGAGGAGGCGTAAAGGTTTACCGCCGTCACCGCATAC‑GACC‑AGTTAAAACCAGGCGGCGGG‑ATGTTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGATCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACTCCACACGCGAACGCAGCCATTGTCCg                                                                                                         <  1:453645/251‑1 (MQ=255)                                                                                                                                                ccACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGc                                                                                                       <  1:275360/251‑1 (MQ=255)                                                                                                                                                          aCGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGc                                                                                             >  1:457803/1‑251 (MQ=255)                                                                                                                                                                          ggTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCGGGCGGGGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAg                                                                             <  1:20400/251‑1 (MQ=255)                                                                                                                                                                           gTGAACTCCGGCACCGCATAC‑GACC‑AGTTAAATCCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAgg                                                                            <  1:184739/251‑1 (MQ=255)                                                                                                                                                                                ccccGTCACCGAATACCGAC‑‑AGTTAAAACCTGGCGGAGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGAc                                                                       <  1:34083/249‑1 (MQ=255)                                                                                                                                                                                    gTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAAt                                                                   >  1:565225/1‑251 (MQ=255)                                                                                                                                                                                              aC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCt                                                         <  1:237648/251‑1 (MQ=255)                                                                                                                                                                                               c‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTa                                                        <  2:261737/251‑1 (MQ=255)                                                                                                                                                                                                   cc‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTaaaa                                                     <  1:68860/251‑1 (MQ=255)                                                                                                                                                                                                       gTTAAACCCTGGGGGGGGG‑ATTGTCACCGCCTAGCTGGCGTCGAGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAagag                                                  <  2:100944/251‑1 (MQ=255)                                                                                                                                                                                                          aaaaCCGGGCGGGGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGaaa                                               <  2:13283/251‑1 (MQ=255)                                                                                                                                                                                                               cgggcgggGGG‑TTTGTCGCCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGcc                                          <  2:526828/249‑1 (MQ=255)                                                                                                                                                                                                                    ggcggG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGGAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTtctc                                     >  1:445761/1‑251 (MQ=255)                                                                                                                                                                                                                     gatgg‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATAGGTGGATGAAACGGACTATTTCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGAGTAAGCCGCGTTGCGGCTGCCGCATCCCTAACCACACCACGCGCGAACGCAGCCTTTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTctct                                    <  2:432577/248‑1 (MQ=255)                                                                                                                                                                                                                       ggg‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCttt                                  >  1:318275/1‑251 (MQ=255)                                                                                                                                                                                                                       ggg‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCttt                                  <  1:363689/251‑1 (MQ=255)                                                                                                                                                                                                                           aTTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTtctc                               >  1:269868/1‑251 (MQ=255)                                                                                                                                                                                                                               ccACCGCCGCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGgcg                           <  2:445761/250‑1 (MQ=255)                                                                                                                                                                                                                                    gCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCAc                      <  2:468060/251‑1 (MQ=255)                                                                                                                                                                                                                                      cTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCATTCTCTTTCTCGGCGCTCACtt                    >  2:150550/1‑251 (MQ=255)                                                                                                                                                                                                                                        cGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTc                   <  1:421521/250‑1 (MQ=255)                                                                                                                                                                                                                                        cGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTc                   <  1:46539/250‑1 (MQ=255)                                                                                                                                                                                                                                         gCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACCCGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGGGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAA‑GCGCCTTCTCTTTCTCGGCGCTCACTTCGcc                <  1:96303/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGc                 >  2:460598/1‑251 (MQ=255)                                                                                                                                                                                                                                          cTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGcc                >  2:539339/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cgcgCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCTGCGCTCACTTCGCCTaaaaa          >  2:187364/1‑250 (MQ=255)                                                                                                                                                                                                                                                 cgcgCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTaaaaa          <  1:346448/250‑1 (MQ=255)                                                                                                                                                                                                                                                 cgcgCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAAt         >  2:428400/1‑251 (MQ=255)                                                                                                                                                                                                                                                    gcAATAAGGGAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGAGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGCTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGcc      >  1:534635/1‑251 (MQ=255)                                                                                                                                                                                                                                                      aaTAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGcc      >  1:482426/1‑249 (MQ=255)                                                                                                                                                                                                                                                       aTAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGcc      >  1:107011/1‑248 (MQ=255)                                                                                                                                                                                                                                                         aaGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGACCGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCTCTTCGCCTAAAAATGCCagag  <  2:353307/250‑1 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      CATCGGGTGTGATCACATCGGCTAAAAAACGTTTGTAACGCTGAATTAGCGTCGCGCGCTGTAGAGGGGGAGAAAATTCCATTGCGACAATCCTTGTTATTGCGTTAGCGCCCAGCGTTTTAGCGGCGTGT‑AGCGTGTGCGTCCACGGGCAAACGAGGAGGCGTAAAGGGTGAACTCCGTCACCGCATAC‑GACC‑AGTTAAAACCTGGCGGCGGG‑ATTGTCACCGCCTCGCTGGCGTCGCGCAATAAGGTAATATGTGGATGAAACGGACGATTGCTTTGAAAACAACCGCTGCGGGCAGCCTGTGAACGGAGCATATTCGCCAGCTGGATTAAGCCGCGTGGCGGCTGACGCATCCCTAACCACACCACACGCGAACGCAGCCATTGTCCGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAG  >  NC_000913/161405‑161898

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