Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,731,666	0	G	A	24.5%	113.6 / 11.9	53	G179E (GGA→GAA)	xylF	xylose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base A (0/13);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 6.32e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTAACGTTGTAAAAGAAGCCAAACAAGAAGGCATTAAAGTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAGGTATTAAACGTATTGCTGCCGGTA  >  NC_000913/3731417‑3731902                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   gTAACGTTGTAAAAGAAGCCAAACAAGAAGGCATTAAAGTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCgg                                                                                                                                                                                                                                                    >  2:539453/1‑250 (MQ=255)                   ccAAACAAGAAGGCATTAAAGTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTggg                                                                                                                                                                                                                                 <  2:532768/251‑2 (MQ=255)                       aCAAGAAGGCATTAAAGTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAGAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TTAcc                                                                                                                                                                                                                            <  1:356364/251‑1 (MQ=255)                       aCAAGAAGGCATTAAAGTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAGAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGAcc                                                                                                                                                                                                                            <  1:356365/251‑1 (MQ=255)                       aCAAGAAGGCATTAAAGTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGAcc                                                                                                                                                                                                                            <  2:92007/251‑1 (MQ=255)                               gCATTAAAGTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGtt                                                                                                                                                                                                                    <  1:478856/251‑1 (MQ=255)                                       gTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCACGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTg                                                                                                                                                                                                          >  2:343281/1‑251 (MQ=255)                                         aTTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGt                                                                                                                                                                                                        >  1:503997/1‑251 (MQ=255)                                            atcTTACGACCGTATGATTAAAGAGGCGGATATCGATTTTTATATTTCTGTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑Ta                                                                                                                                                                                                     <  2:80417/248‑1 (MQ=255)                                            aGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TAc                                                                                                                                                                                                    <  1:346063/251‑1 (MQ=255)                                             gCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TAcc                                                                                                                                                                                                   <  1:199257/251‑1 (MQ=255)                                               ttACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCgg                                                                                                                                                                                                 >  1:428685/1‑251 (MQ=255)                                                         aTGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTAGTTCCTGAGGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCAtt                                                                                                                                                                                      <  2:264313/251‑1 (MQ=255)                                                          tGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGt                                                                                                                                                                                                                                        >  2:14534/1‑205 (MQ=255)                                                          tGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGt                                                                                                                                                                                                                                        <  1:14534/205‑1 (MQ=255)                                                                  gATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TAt                                                                                                                                                                            >  1:417915/1‑251 (MQ=255)                                                                     gCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGa                                                                                                                                                                         <  1:247899/251‑1 (MQ=255)                                                                        gcgatCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGTCAAAAGCCCTGGTCGATATTGTTCCGCACGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCTCCGGACAAATGTAAGTGTTACAACCTTACGTTGATTCCGGCACAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACAGGAAAAC‑GCATTGAAAAT‑TATGGaac                                                                                                                                                                       >  2:268717/4‑249 (MQ=255)                                                                         atatCGATTTTTATATTTCCTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGaaaaa                                                                                                                                                                     <  2:159418/251‑2 (MQ=255)                                                                                  tttATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAAc                                                                                                                                                             <  1:367157/250‑1 (MQ=255)                                                                                       atTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCaa                                                                                                                                                       <  1:355520/251‑1 (MQ=255)                                                                                       atTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCa                                                                                                                                                        <  2:361085/250‑1 (MQ=255)                                                                                               cGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAAc                                                                                                                                                >  1:540116/1‑250 (MQ=255)                                                                                                           aGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATAGTTCCGCACGGTAA‑TTACTTCCGGATGGGCGGATCGACGGAAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTACGGAAAAATTAAAGTCGTTGG‑TGACCAAGGGGTG‑G‑ATGGCGGGT‑TACAGGAAAAC‑GCATTGAAAAT‑TATGGAAAAGGAGATAACGACAAATAATAACAAAATTGATGCTg                                                                                                                                   >  2:477649/1‑251 (MQ=255)                                                                                                              gggTGAACGGAAGGAAAAACCCCTGGTGGATATTGTTCCGCA‑GGTAAATTCCTCCCTGTGGGGCGGCTCCCCGGTAGTCACCAACGCAAGGCTGTCCCGCGCGGGACAAATGAAAGTGTTAACACCTTGCGTTGATTCGGGAAAAATTAAAGTCGTTGT‑AGACCAATGGTTT‑G‑ATGGCGTGT‑TCCCGGAAAAC‑GCATTAAAATT‑TATGGAAAACCCGCTAACCGCCAATAATAACAAAATTGATGCTGTAg                                                                                                                                <  1:256107/250‑1 (MQ=255)                                                                                                               gggaaaCTGCAGCCAAAAGCCCTTGTCGATATTGTCCCGCAATGTAT‑TTACTTCCTGATGGGCGACGCCCCTGTAGATAACAACGCCAAGCTTTCCCGCGCCGAACAAATGAAAGTTTTAAACCCTTACGTGGAGTCCGAAAAAATTAAAGTTTTTGT‑GGACCAATGGGTT‑G‑ATGGCTGTT‑TCCCGAAAAAC‑GCATTAAAATT‑TATGAAAAACGGGCTAGCCGCCAATAATAACAAAATTGATGCTGTAGt                                                                                                                               <  1:354891/247‑1 (MQ=255)                                                                                                                       gCAGGCAAAGGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTCCTTCCTGTTGGGCGGCTCGCGGGTAAATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAATTTTTAAAACCTTACGTTGATTCGGAAAAAATTAAATGCGTTGG‑TGACCAATGTGTT‑G‑ATGGCTGGT‑TCCCGAAAAAC‑GCATTTAAAAT‑TAGGGAAAACTCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTc                                                                                                                         <  1:157852/249‑1 (MQ=255)                                                                                                                        cAGGCAAAAGCCCTGGACGATTTTGTTCCGCAAGGTAA‑TTACTACCTGATGGGCGGCTCGGCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAACTGAAAGTGTTCAAACATTACGTTGATTCCGGACAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAG‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCaa                                                                                                                       >  2:432814/1‑250 (MQ=255)                                                                                                                         aGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAAc                                                                                                                     >  2:27318/1‑251 (MQ=255)                                                                                                                          gcaaaaaGCCCTGGTGAATATTGTTCGGAAAGGTAT‑TTACTTCGTAGGGGGCGGTTCGCCGTTAGATAACAACGGCAAGCTGTTCCGCGCGGAACAAAGAAAAGTGTTAAAGCCTTATGTTGATTCCGGAAAAATTAAAGTCGTTGG‑GAACCAAGGGGTT‑G‑ATGGCTGGT‑TACCGAAAAAC‑GCATATAAAAT‑TAGGGAAACCGCGATAGCCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACg                                                                                                                    <  1:490265/248‑1 (MQ=255)                                                                                                                              aacccccTGATCGATATTGTCCCGCAAGGTAA‑TTACTCCCTGGTGGGCGTCTCGCGGGTAAAAAAAAACGCCAATCTGTTCCGCGGCGGACAAATGAAAGTGTTAAAACCTTACGTGTATCCCGAAAAAATTAAAGTCGTTGT‑TGAACAGTGGGTG‑G‑ATGGCTGGT‑TACCGAAAAAC‑GCATTAAAAAT‑TATGGAAAACCCGCAAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGc                                                                                                                <  2:93805/247‑1 (MQ=255)                                                                                                                                aGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGcc                                                                                                               >  1:105125/1‑250 (MQ=255)                                                                                                                                               gttcccgaaaGGTAT‑TTCCTGCCTGATGGGGGGCGCGCCGGTAGAAAACAACGCCAAGCTGTCCCGCGAGGAAAAAATGAAAGTGTTAAATCCTTACGTTGTTTCGGAAAAAATTAAAGTCGTGGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGAAAAAC‑GCATTGAAAAT‑TATGAAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCaa                                                                                               <  2:13985/243‑1 (MQ=255)                                                                                                                                                gTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCACGCCGGTAGATAACACCGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAACGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTCGTTGCCTCAAACGATGCCACCGCAGGTGGGTCAAt                                                                                              >  2:384930/1‑251 (MQ=255)                                                                                                                                                  cccgaaaGGTAA‑TTACTTCCGGAGGGGGGGCCCGCCGTTAGATAACAACGCCAAGTTGTCCCGCGCCGGACAAATGAAAGTGTTAAACCTTTACGTTGATTCCGAAAAAATTAAAGTCGTGGG‑TGCCCAATGGGTT‑G‑AGGGCTGGT‑TACCGAAAAAC‑GCATGGAAAAT‑TATGGAAAACCCGATAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTc                                                                                            <  2:105125/246‑1 (MQ=255)                                                                                                                                                   tccgaaggtaaattgcatcctagtGGGCGGCTCGCCTGTAGATAACCACGCCAACCTGTCCGGCCCCAGACAAATGAATTTGTTAAATCCTTACGTTGGTTCGGAAAAAATTATAGTCTTTGT‑CGGCCAATGGGTT‑G‑CTGGCTTGT‑TACCGAAAAGC‑GCTGTAAAAAT‑TATGAAAGACCCGCAAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAAtt                                                                                             <  2:366108/227‑1 (MQ=255)                                                                                                                                                                       gggtGGGCGCCTCGCGGGTATAAAACACCGCCACGCTTTCCCGGGCCGAAAAAAGAAAATTGTAAAAACTTTACGTTGATTCCGAAAAAATAAAAGTCGTGGG‑TGACCAATGGGTT‑T‑AGGGGTGGT‑TCCCGAAAAAC‑GAATGAAAAAT‑TATGGAAACCGCGCTACCCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGt                                                                        <  2:368021/248‑1 (MQ=255)                                                                                                                                                                              cggcTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATc                                                                   >  1:408325/1‑249 (MQ=255)                                                                                                                                                                                  gctccGGTAAAAAACACCGCCAAGCTGTTCCGAGCCGGAAAAATAAAATTGTAAAAACCTTCCGTTGATTCCGAAAAAATAAAAGTCGTTGT‑TGCCCAA‑GGGTT‑GGAGGGTTGGT‑TACCG‑AAAACTGCATTAAAAAT‑TAGGAAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGAGGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGaa                                                             <  2:438428/248‑1 (MQ=255)                                                                                                                                                                                  acgcccgtaaaaaacaccGCAAAGCTGTTCCGCGCCAGAAAAAGGAATGTGGTAAACCCTCACTTAGTTTCCGAAAAAATGAATGGCGTTGG‑TGCCCATGGGGTG‑G‑ATGGCGGGT‑CCCCGAAAAAC‑GCAT‑AAAAATGTATGAAAACCGCGCAACCCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGTAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGaa                                                             <  2:9816/239‑1 (MQ=255)                                                                                                                                                                                           gATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAAt                                                    >  2:17666/1‑251 (MQ=255)                                                                                                                                                                                                     ccAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTGACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGACAGTAGCCATCTCCGGcc                                            >  2:339962/1‑249 (MQ=255)                                                                                                                                                                                                              cccGCGCGGGACAAATGAAAGTGATAAACCCTTGCGTGGATTACGAAAAAATTAAAGGCGTTGG‑TGACAGATGGGTG‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATAGAAAAT‑TATGGAAAACGGGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGAt                                 <  2:355508/250‑1 (MQ=255)                                                                                                                                                                                                              cccGCGCCGAACAAATAAAAGTGTTAAAACCTTCCGTTGATTCCGAAAAAATAAAATTCGTTGG‑TTACCAATGGTTG‑G‑ATGGCTGGT‑TACCGGAAAAC‑ACTTTGAAAAT‑TAGGGAAAACGCGCGATCCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGTAATTCAGGCATTAAGCGCGCAATGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGa                                  <  1:455340/249‑1 (MQ=255)                                                                                                                                                                                                                cgcgcgtgaaaaaggAAAGTGTTAAAACTTTACGTGGATTCAGGAAAAAGTAAGTTGGTTGG‑TG‑CCAAAGGGTTTG‑ATGGC‑GGTGTACCGGAAAAC‑GCATGAAAAAT‑TATGGAAAACGCGCTACCCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGAtct                               <  2:105390/237‑1 (MQ=255)                                                                                                                                                                                                                    ccGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGACGCGGATCTCGc                            >  1:77447/1‑250 (MQ=255)                                                                                                                                                                                                                      ggACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAgg                         >  1:357506/1‑251 (MQ=255)                                                                                                                                                                                                                      ggACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGACACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAgg                         >  1:523639/1‑251 (MQ=255)                                                                                                                                                                                                                       gACAAATAAAAGTGTTAAACCCTTACGTTGATTCCGGAAAAATAAAAGTCGTGGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGAAAAAC‑GCATTAAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAGGt                        <  1:27318/251‑1 (MQ=255)                                                                                                                                                                                                                        aaaaatgaaatgGTTAAACCCTTACGTTGATCCCGAAAAAATAAAAG‑CTTTGGGTGACAAATGGGTT‑G‑ATGGCGGTT‑TACCGAAAACC‑GCATAAAAAAT‑TATGAAAAGCGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAGGTa                       <  2:428685/249‑1 (MQ=255)                                                                                                                                                                                                                            aTGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCACTGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAGGTATTa                    >  2:21613/1‑250 (MQ=255)                                                                                                                                                                                                                               aaaGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGCAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGAGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAGGTATTAAACg                >  2:1146/1‑251 (MQ=255)                                                                                                                                                                                                                                       aaaaCCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTAGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAGGTATTAAACGTATtgctg        >  2:54522/1‑251 (MQ=255)                                                                                                                                                                                                                                              tacttggATTCCGAAAAAATTAAAGTCGTGGG‑TGCCCAAGGGTTT‑G‑AGGGCTGTT‑AACCGAAAAAC‑GAATAGAAAAT‑TAGGGAAAACGCGCAACCCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAGGTATTAAACGTATTGCTGCCGGTa  <  2:149588/244‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                   GTAACGTTGTAAAAGAAGCCAAACAAGAAGGCATTAAAGTATTAGCTTACGACCGTATGATTAACGATGCGGATATCGATTTTTATATTTCTTTCGATAACGAAAAAGTCGGTGAACTGCAGGCAAAAGCCCTGGTCGATATTGTTCCGCAAGGTAA‑TTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGG‑TGACCAATGGGTT‑G‑ATGGCTGGT‑TACCGGAAAAC‑GCATTGAAAAT‑TATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCGCAAGGTTTATCAGGGAAAGTAGCAATCTCCGGCCAGGATGCGGATCTCGCAGGTATTAAACGTATTGCTGCCGGTA  >  NC_000913/3731417‑3731902

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