Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	21,520	0	C	A	73.1%	117.5 / 37.9	78	G38G (GGC→GGA)	ribF	bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (12/9);  new base A (33/24);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAA  >  NC_000913/21274‑21763                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    ggTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGc                                                                                                                                                                                                                                                 <  1:262722/251‑1 (MQ=255)    cAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAgg                                                                                                                                                                                                                                              <  2:524030/251‑1 (MQ=255)       tttCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTATGATTTCGGTACTCAATCGCAGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGga                                                                                                                                                                                                                                             >  1:558764/1‑249 (MQ=255)                 tGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGAGCGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCCGGAAGAAGGgcgca                                                                                                                                                                                                                                  >  2:128504/1‑250 (MQ=255)                           cATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGTTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACtt                                                                                                                                                                                                                        >  1:453160/1‑250 (MQ=255)                           cATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACtt                                                                                                                                                                                                                        >  1:14633/1‑250 (MQ=255)                            aTGGCCATGAACAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGATGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATACTCTCAGCCAGGCCCCGAA‑AGAAGGGTGTGTGCTGACTCTTGGTAACTTCGACGGCGGGCATAGCGGTCATCGCGCGATTTTTCAGGGATTGCAGGAAGAAGGGCGCAAGCGCGACTTAc                                                                                                                                                                                                                      >  2:430784/1‑251 (MQ=255)                             tGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACACGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTAcc                                                                                                                                                                                                                     >  1:629756/1‑251 (MQ=255)                              ggCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCg                                                                                                                                                                                                                    >  1:137119/1‑251 (MQ=255)                               gCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCC‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCgg                                                                                                                                                                                                                   >  1:3696/1‑251 (MQ=255)                                     aaaaGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCggtgatgg                                                                                                                                                                                                             >  2:434669/1‑251 (MQ=255)                                      aaaGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCggtgatggt                                                                                                                                                                                                            <  1:185542/251‑1 (MQ=255)                                      aaaGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCggtgatggg                                                                                                                                                                                                            >  2:158424/1‑250 (MQ=255)                                       aaGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCggtgatggtg                                                                                                                                                                                                           >  1:295247/1‑251 (MQ=255)                                           cACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGGAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAGGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGc                                                                                                                                                                                                       >  1:215740/1‑251 (MQ=255)                                             cttcttTCTGCTTTCGGTCCTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACACGAAGCTGCTACGCGGCCTACATAATCTCAGCCAGGCCCCGCA‑AGACGGGTGAGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCATGCCGGAAGAAGGGCGCAATCGCAACTTACCGGTGCTGGTGTTGCtt                                                                                                                                                                                                     >  2:400775/1‑251 (MQ=255)                                             cttcttTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCtt                                                                                                                                                                                                     >  1:15034/1‑251 (MQ=255)                                               tctttctGATTTCGGTACTCAATCGCCGGTTAACCTTGGCCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTGACAGGGATTGCAGGAAGAAGGGAGCAAGCGCAACTTACCGgtgatggtgatgttt                                                                                                                                                                                                     <  2:59180/249‑4 (MQ=255)                                               tctttctGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCtttt                                                                                                                                                                                                   >  2:275416/1‑251 (MQ=255)                                                     tGATTTCGGTCCTCAATCGCCGGTTAACCCTGACCGCTGTACAAAGCATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑TAAAGGGCGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCACGCGCAACTTACCGGTGATGGTGATGCTTTTTGAAcc                                                                                                                                                                                             >  2:295480/1‑251 (MQ=255)                                                       aTTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTc                                                                                                                                                                                           <  1:509658/251‑1 (MQ=255)                                                           cGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAAc                                                                                                                                                                                        <  1:508054/250‑1 (MQ=255)                                                               aCTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACt                                                                                                                                                                                    <  1:156470/250‑1 (MQ=255)                                                                  cAATCGCCGGTTAACCTTGACCGATGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGCAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCACCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGa                                                                                                                                                                                 >  2:530511/1‑250 (MQ=255)                                                                    aTCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGCAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACt                                                                                                                                                                              >  2:413829/1‑251 (MQ=255)                                                                           gTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACCTGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGACATGTTTGCt                                                                                                                                                                       >  1:191074/1‑251 (MQ=255)                                                                              aaCCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTAcc                                                                                                                                                                    <  2:254393/251‑1 (MQ=255)                                                                               aCCTTGACCGCTGTACAAGGTATCCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTAATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACAGAAACTGTTTGCTACCg                                                                                                                                                                   >  2:369376/1‑251 (MQ=255)                                                                                 cTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGAt                                                                                                                                                                 >  1:493910/1‑251 (MQ=255)                                                                                    gACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGCCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATaa                                                                                                                                                               >  2:61033/1‑250 (MQ=255)                                                                                     aCCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAg                                                                                                                                                             <  1:524717/251‑1 (MQ=255)                                                                                                  gTATACTCGGACGCTTTTCACTGTTTTGAGCCAGACATGAAGCTGGTACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGCAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCACGCGCAACTTACAGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCacggct                                                                                                                                                >  1:33863/1‑246 (MQ=255)                                                                                                    ataCTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGTTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCacggctg                                                                                                                                               <  2:453160/250‑7 (MQ=255)                                                                                                              gATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCg                                                                                                                                    >  1:529638/1‑251 (MQ=255)                                                                                                               aTTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCgg                                                                                                                                   >  2:228554/1‑251 (MQ=255)                                                                                                                 tttCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCggg                                                                                                                                  >  2:550006/1‑250 (MQ=255)                                                                                                                     aCTGTTATGATCCAGACATGAAGCTGATACGCGGCCTACATACTCTCAGCCAGGCCCCGCA‑CGAAGGGTGTGTGCCGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGCAGACGGGCGCAATCGCAACTTACCGGTGCTGGTGATGCTTTTTGAACCTCACCCACTGGAACTCTTTGCTACCGATAAAGCCACGGCCCGGCTGACCCGGcaccgggacaa                                                                                                                              >  2:254172/1‑240 (MQ=255)                                                                                                                          tttGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGCCGCCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCg                                                                                                                        >  1:365422/1‑251 (MQ=255)                                                                                                                              aGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTAc                                                                                                                    <  2:505693/251‑1 (MQ=255)                                                                                                                                         aaGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAgtg                                                                                                         >  2:523873/1‑251 (MQ=255)                                                                                                                                          aGCTGATACGCGGCATACATACTCTCAGCCACGCCCCGCA‑AGACGGGCGCGTGCTGACTATTGGTAATTTCGACGGCGTGCATAGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCACCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCACCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTCCCTTGCCGAgtgg                                                                                                        >  1:525842/1‑250 (MQ=255)                                                                                                                                             tGATACGCGGCATACATATGCTCAGCCAGTCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGGAATGTCGACGGCGTGCATCGCGGTCCTCGCGCGCTGTTACAGGTCTTGCATGAAGAAGGGCGCAAGCGCAACTTTCCGGTGATGGGGATGCGTGGTGAACCGCAACCGCTGGAACTGTTTGCAACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGGGACCTGGCAGAGtgtg                                                                                                       <  2:625991/249‑1 (MQ=255)                                                                                                                                             tGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGCTGGTGATGCTTTTTGCATCTCAACCACTGGAACTGTTTGCTACCGCTAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGtgtg                                                                                                       >  1:625991/1‑249 (MQ=255)                                                                                                                                               aTACGCGGCATACCTAATCTCAGCCAGGCCTCGCA‑AGAAGGGTGTGTGCTGGCTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTCCCATGCAGAGTGTGGCGt                                                                                                   <  2:350668/251‑1 (MQ=255)                                                                                                                                                  cgcgGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGa                                                                                                <  2:627883/251‑1 (MQ=255)                                                                                                                                                        atacatcATCTCAGCCAGGCCCCGCA‑AGACGGGTGTGTGCTGACTATTGGTCATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACgg                                                                                          >  2:416804/1‑250 (MQ=255)                                                                                                                                                          acataATCGCAGCCAGGCCACGAACGGAAGGG‑GGGTGCTGACTATTGATAATTTCGACGGCGTGCAGAGCGGTCATCGCGCGCTGTTACAGGGATAGCAGGAAGAAGGGCGCAAGCGCAACTTCCCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACACGGCTGCGGGAAACACTGCGTTACCTTGCAGAGGGGGGCGTTGATTACGTg                                                                                         >  2:540151/1‑250 (MQ=255)                                                                                                                                                           cataATCTCAGCCTGGCACCGCA‑AGATGGGTGTGTGCTGACTAGTGGTAATTGCGACGGCTTGCATCGCGGCCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTAACGGTGATGTTGATGCTTTTTGAACCTCAACCATTGGAACTTTTTGCTACCAATAAAGCCACGGCACGGCTGACCCGGCTGCGTAAAACACTGCGTTACCTTGCGGAGTGTGGCGTTGATTACGTGCt                                                                                       <  1:214665/251‑1 (MQ=255)                                                                                                                                                           cataATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCt                                                                                       <  2:70939/251‑1 (MQ=255)                                                                                                                                                                  ccAGCCATGCCCCGCA‑AGAAGGGTGTTGGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCAGCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCATGCTCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTTCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCTTTGATTACGTGCTgtgcgtt                                                                                <  2:378376/250‑2 (MQ=255)                                                                                                                                                                   cAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTgtgcgtgc                                                                               <  1:431343/251‑1 (MQ=255)                                                                                                                                                                   cAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTgtgcgt                                                                                 <  2:243981/249‑1 (MQ=255)                                                                                                                                                                   cAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTgtgcgtgc                                                                               <  2:3413/251‑1 (MQ=255)                                                                                                                                                                    aGCCAGGCCCCGAA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTCCGTGCTGtgcgtgcg                                                                              >  2:490240/1‑251 (MQ=255)                                                                                                                                                                      ccAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGGGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCATAGTGTGGCGTTGATTACGTGCTGTGCGTGCGtt                                                                            <  2:518426/251‑1 (MQ=255)                                                                                                                                                                         ggCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCg                                                                         >  1:24440/1‑251 (MQ=255)                                                                                                                                                                              cGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAgg                                                                    <  1:623833/251‑1 (MQ=255)                                                                                                                                                                                       gggTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGGTTCGACAGGCGTTTcgc                                                            <  1:430784/251‑1 (MQ=255)                                                                                                                                                                                          agtgtgCTGCCTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAGGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGAGCGTGCGTTTAGACAGGCGTTTCgcggc                                                         >  2:582260/2‑251 (MQ=255)                                                                                                                                                                                           gtgtgCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCgcggcg                                                        >  1:133570/1‑251 (MQ=255)                                                                                                                                                                                             gtgCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGtt                                                      <  1:48508/251‑1 (MQ=255)                                                                                                                                                                                             gtgCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGt                                                       <  1:360800/250‑1 (MQ=255)                                                                                                                                                                                              tgCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAGGCCACGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGAGCGTGCGTTTCGACAGGCGTTTCGCGGCGtt                                                      >  1:89633/1‑250 (MQ=255)                                                                                                                                                                                                      aTTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCa                                             <  1:434669/251‑1 (MQ=255)                                                                                                                                                                                                        tGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCaa                                            <  2:14633/250‑1 (MQ=255)                                                                                                                                                                                                             aTTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTtc                                      <  1:156420/251‑1 (MQ=255)                                                                                                                                                                                                              tttCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCtgt                                                                                     >  2:109986/1‑203 (MQ=255)                                                                                                                                                                                                              tttCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCtgt                                                                                     <  1:109986/203‑1 (MQ=255)                                                                                                                                                                                                               ttCGACGGCGGGCATCGCCGTCCTCGCGCGCTGTTACAGTGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACTTCAACCACTGGAACTGTTTGCTACCTATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAATACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCTTGCGTTTCGACATGCTTTTCGCGGCGTTAACCGCGCAAAATTtca                                     <  1:128133/250‑1 (MQ=255)                                                                                                                                                                                                                               gcgGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTg                    >  1:480377/1‑251 (MQ=255)                                                                                                                                                                                                                                     aTCGCGCGCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAACATTTCATCAGCGATCTTCTGGTGAAGctt              >  2:472167/1‑249 (MQ=255)                                                                                                                                                                                                                                       cgcgcgCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGcttt             <  2:24440/250‑4 (MQ=255)                                                                                                                                                                                                                                       cgcgcgCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCtttt            <  2:529638/251‑1 (MQ=255)                                                                                                                                                                                                                                       cgcgcgCTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCtttt            >  1:481855/1‑251 (MQ=255)                                                                                                                                                                                                                                       agcgcgcTGATACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGCGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGGGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTACATCAGCGATCTTCTggggaggctttt            >  2:483275/2‑241 (MQ=255)                                                                                                                                                                                                                                         cgcgcTGTTACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTgc          <  2:137119/251‑1 (MQ=255)                                                                                                                                                                                                                                               gTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCACCCACTGGCACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCCTCAGCGATCTACTGGTGAAGCTTTTGCGCGTa     >  2:250817/1‑250 (MQ=255)                                                                                                                                                                                                                                                 tACAGGGATTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGGGCGTTTCGCCAGGCGTTTCGCGGCGTTAACCGCGCAAACTTTCATCAGCGATCTTCTGGTGAGGCTTTTGCGCGTaaat  >  2:603366/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GGTCAGTTTCTTCCCGTGATTTGCGCCATGGCAATGAAAAGCCACTTCTTTCTGATTTCGGTACTCAATCGCCGGTTAACCTTGACCGCTGTACAAGGTATACTCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCA‑AGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAA  >  NC_000913/21274‑21763

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