Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,101,673	0	T	A	34.0%	95.3 / 14.8	53	L4F (TTA→TTT)	csgG	curli secretion channel
Reads supporting (aligned to +/- strand):  ref base T (10/25);  new base A (18/0);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 2.11e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTCGACAGTAGCCCACCTAAAATTTGTGACTGGATGGCCTGAGTAAAGTTAT  >  NC_000913/1101430‑1101922                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                  cGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCt                                                                                                                                                                                                                                                            <  1:325995/251‑1 (MQ=255)        agaattgcctttgcgctaggCGGAACAGCAGGGGAGAAGTTACTTGACGGGTAGGGTTTAAATTGCCCGTTTTCGTCCTGAAGGGTGTATTCCGAAACAAAGATTGTACCCGTCGTCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCTCGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGTCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGTCAACCAAAATAAATAAGCGCTGC‑ATg                                                                                                                                                                                                                                                      <  2:452251/230‑1 (MQ=255)            ttGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑GTAAGCCTCCGCCCA‑GTAAAAAGA‑CGGCCACCAAAAGAAAAAAAAGCTGC‑ATGaat‑att                                                                                                                                                                                                                                               >  1:632371/1‑251 (MQ=255)             tGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGatt‑att                                                                                                                                                                                                                                               <  1:301675/250‑1 (MQ=255)                        tttGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGGACCCCCTCCCATAACA‑GG‑CGGCAACCAAAAAAANAAAAAGCTGC‑TTGATT‑ATTTTTCtttttca                                                                                                                                                                                                                                    >  2:351830/1‑245 (MQ=255)                            cGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCAACCGCTCA‑GTAACAAGA‑CGGCAACCAAAAAAAAAAAAAGCTGC‑ACGAGT‑ATTTTTCTTATGAAACTgg                                                                                                                                                                                                                               >  2:653613/1‑251 (MQ=255)                             ggAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTggg                                                                                                                                                                                                                              <  1:589034/251‑1 (MQ=255)                                        ggAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑Tc                                                                                                                                                                                                                 <  2:11568/251‑1 (MQ=255)                                          agaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCgg                                                                                                                                                                                                               <  2:484419/251‑1 (MQ=255)                                                tACTTGCCGGGTAGGGTTTAAAGTGCCCGGTTGTGTCCTGAATGTTGGATACCGAAACAAAGATATTACCCGGCGGCGCTGGCAGATGGGTCAAATCTTTGTAGGTCTGAGCACGAGGCATTAATGTCGGTATGGCGGCGTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTGA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCGGC‑ATGATT‑ATTGTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAg                                                                                                                                                                                                          <  2:379836/250‑1 (MQ=255)                                                  cTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTTAGAATCCGCTCC‑GGAACATGA‑CGGCAACCCAAAAAAAAAAACGCCGC‑ATGATT‑TTTTTTCTTATGAAAATGGGG‑CCTAAAAA‑ACgggtgggttt                                                                                                                                                                                                       >  1:629485/1‑243 (MQ=255)                                                             gggTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCCCCCA‑GTAACATGA‑CGGCAACACACAGAAAAAAACGCTGC‑CTGATT‑ATTTTTCTTTTGAAACTGGGG‑GTTCAAAA‑ACGGTTGGG‑TTTTTTTTTAAAcc                                                                                                                                                                                           >  2:302380/1‑251 (MQ=255)                                                               gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCC‑GTAACATTA‑CGGGAAACAAAAGAAAAAAGCGCCGC‑ATGATT‑ATTTTTCTTTTGAAGCTGGGG‑GTTAAAAA‑TCGGTTTAG‑TTTTTTTGTaaaaccc                                                                                                                                                                                         >  1:375577/1‑247 (MQ=255)                                                                   aaaTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCCCA‑GTAACATTA‑AGGCAACCCAAAGAAATAAACGCTGC‑AAGATT‑ATTTTCCTTTTTAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTTTAAACCCGAAAc                                                                                                                                                                                     >  1:167201/1‑251 (MQ=255)                                                                          ccGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGATGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCGCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGCAAACCCGAAACCT‑GGATg                                                                                                                                                                             <  2:604250/251‑1 (MQ=255)                                                                             gTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGGTTAGAAATCC‑CGCC‑GAAACAAGA‑CGGCAACCAAAAAAAAAAAAGGCTGCCAT‑ATT‑TTTTTTCTTATGAAACTTGGG‑GTTAAAAA‑ACGGGTGNG‑TTATTTTTTAAAACCGAAACCT‑ggggtggg                                                                                                                                                                          >  2:353621/1‑245 (MQ=255)                                                                                 cGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAACATCCGCCCC‑GTAACAAGA‑CGGCAACCAAAAAAAAAAAACGCTGC‑ATGATT‑ATTTTCCTTTTGAAGCTGGGG‑CTTAAAAA‑ACGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAgg                                                                                                                                                                      >  1:45863/1‑251 (MQ=255)                                                                                  gTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑GTAAACATCCGCTCC‑GTAACATGA‑CGGCAACCAAAAGAAAAAAGGGCTGC‑AAGATT‑TTTTTCCCTATGAAGCTGGGG‑CTTAAAAA‑TCGGGTGAGGTTTTTTTTTAACCCCGAA‑CCC‑GGATgggcggggg                                                                                                                                                                     >  1:480187/1‑244 (MQ=255)                                                                                   tCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCCCCCC‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTTCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAAACCGAAAACC‑GGATGGTCGAgggt                                                                                                                                                                    >  1:9453/1‑249 (MQ=255)                                                                                     cTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑GTAAAGATCCGCCCC‑ATAAAAAGAAGGGCAAC‑AAAAGAAAAAAGAGCTTC‑ATGGTT‑TTTTTACTTATGAAACAGGGG‑GTTAAAAAATGGGTT‑AG‑GTATTTTGTAAACCCCAAAACC‑GGGTGGTCGAgggttt                                                                                                                                                                  >  2:591750/1‑247 (MQ=255)                                                                                         aTGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCg                                                                                                                                                             <  1:72579/251‑1 (MQ=255)                                                                                         aTGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑ATAAGAAACCGCTCA‑GGAAAAAGA‑CGGGCACCAAAAAAAAACAACGCTGC‑AAGATT‑ATTTTTCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGGTGGG‑GTATTTTTGAAACCCGAAAACT‑GGATGGTCGGGG‑TGTGTCCg                                                                                                                                                             >  1:523305/1‑251 (MQ=255)                                                                                               tataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTTAAAATTCCCCCA‑GTAAAAAAT‑AGGCCACCAAAAGAAAAAAACGATAC‑CAGAAT‑ATTTTTCTTTCGAAGCTGGGG‑CTTTAAAA‑TAGGTTTAG‑TTATTTTTTAAACCCGAAACCTAGGATGG‑CGGGG‑GTTTTCCGGttttt                                                                                                                                                       >  1:243672/1‑250 (MQ=255)                                                                                                ataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTAc                                                                                                                                                      <  1:351830/251‑1 (MQ=255)                                                                                                 taCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCCTTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑GTAAGCAACCCCCCA‑CTAAAATGA‑CCGCTACCAAAAGAAAAAAAAGCTGC‑ACTATTCATTTTCCT‑ATTCAGCCCGGG‑CCTAAAAA‑TCCGTTGAG‑GTATTTTGTTAACCCGAAACCC‑GGAGTGTCGAGGGTTTTCCCGGtttt                                                                                                                                                        >  2:504206/1‑249 (MQ=255)                                                                                                     gAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGAGT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑GCGGTTGAA‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCt                                                                                                                                                 <  1:655420/251‑1 (MQ=255)                                                                                                     gAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGAGG‑TTTAGCATCCGCCCA‑ATAACATGG‑CGGCAACCCAAAAAAAAAAACGCCGC‑ATGTTT‑ATTTTTCTTTTGAAGGTTGGG‑GTTTAAAA‑TCGGGTTTG‑TTTTTTTGTAAACCCGCAACCC‑GGGTGGGTGAGGGTTTTTCCGGTTT‑ACGATCt                                                                                                                                                 >  1:402340/1‑251 (MQ=255)                                                                                                           aaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑GTAAGCATTCCCCCC‑GTAAAAAGA‑CGGCAACCAAAAGAAATAAACGCGGC‑ATGTTT‑ATTTTTCTTATTAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTTTTAACCCCAAACCT‑GGGTGGGCGAGG‑TTTGTCCGGtttttcggtcttttcac                                                                                                                                           >  2:143080/1‑238 (MQ=255)                                                                                                            aGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTAGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGt                                                                                                                                         <  2:116395/251‑1 (MQ=255)                                                                                                                             gcTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCa                                                                                                                        <  1:43663/251‑1 (MQ=255)                                                                                                                                   aGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑GTAAAGATCCGCTCT‑GTAACATGA‑CGGCAACCCAAAGAAAAAAACGCCGG‑ATGATT‑ATTTTTCTTTTGAAGCTGGGG‑GTTAAAAA‑TCGGGTGAG‑TTATTTTGTAAAACCGAAAACT‑TGATGGGCGAGG‑GTTGTCCGGTTTTAAGAACTGT‑CACGTTCAACCGCAATTGACAATCGCgg                                                                                                                  >  1:258304/1‑251 (MQ=255)                                                                                                                                             aaTCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAatat                                                                                                        <  1:710440/251‑1 (MQ=255)                                                                                                                                                      aGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAatat                                                                                               <  1:653613/251‑1 (MQ=255)                                                                                                                                                       gCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAAtata                                                                                              <  1:167802/251‑1 (MQ=255)                                                                                                                                                            gAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCTGAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGGCCAGCGCGGTTGGCAATATCGACAATATAATc                                                                                           <  1:265141/249‑1 (MQ=255)                                                                                                                                                            gAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAAt                                                                                            <  2:45863/248‑1 (MQ=255)                                                                                                                                                                        ttAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTTAAAATCCGCCCA‑GGAACATCA‑CGGGAAACAAAAGAAAAAAAAGCTGC‑ATGATT‑ATTTTTCTTTTGGAACTGGGG‑CTTAAAAA‑ACGGGTGAG‑TTTTTTTTTAAACCCGAAACCT‑GGGTTGTCGCGG‑GTTGTTCGGTTTTATGATTCGT‑CACGCTCAACTGCGACTTACCCTCGCGGTTGGGCATATCGACAAAATAATCGTTGGTCACCATGc                                                                             >  1:104490/1‑251 (MQ=255)                                                                                                                                                                          aaTGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCAGCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCg                                                                            <  2:601490/250‑1 (MQ=255)                                                                                                                                                                           aTGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTTAAAATCCGCCCA‑GGAAAATGA‑AGGCAACCAAAAAAAATAAACGCCGC‑CTGATT‑ATTTTCCTTTTGAAACTGGGGGCTAAAAAA‑TGGGTTG‑G‑TTATTTTTTTAAACCGGAACCC‑GGAAGGTTGGGG‑GTTGTTCGGGTTTACGATCTGT‑TACCTTCCACTGCAAATTACCCTCGGGGGTGGCAAAATTGACaaattaaatggtggtccccatgcggg                                                                          >  1:279000/1‑227 (MQ=255)                                                                                                                                                                                ggTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑AGGGAAACAAAAGAAAAAAGCGACGC‑ATGATT‑ATTTTCCTTATGTAGCTGGGG‑CTTAAAAA‑TCGGGTGAG‑TTATTTTGTAAAACCGAAAACT‑GGATGGGCGTGG‑TTTGTCCGGTTTTATGATATGT‑CACGTTCAACTGTAATTGGCAATAGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGAGGGCCGGt                                                                     >  1:133383/1‑251 (MQ=255)                                                                                                                                                                                  tCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAAAATGA‑CCGCAAACAAAAGAAAAAAGAGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑GTTAAAAA‑TCGGTTGAG‑GTATTTTTTAAAACCGAAACCC‑GGGTGGGCGAGG‑TTTTTTCGGGTTTAAGATTTGTACACGT‑CAACTGCAAATGACCCTCCCGGGTGGCAATATCCGCAATATAAACGTTGGTCACCCTTCGGCCCgggtt                                                                   >  1:99636/1‑248 (MQ=255)                                                                                                                                                                                     ggcggcTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGttt                                                                   <  1:710097/248‑1 (MQ=255)                                                                                                                                                                                     ggcggcTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTAcc                                                                <  2:398622/251‑1 (MQ=255)                                                                                                                                                                                        ggcTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGt                                                             <  1:699941/251‑1 (MQ=255)                                                                                                                                                                                                         cgg‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTCGACAGt                                             >  1:12660/1‑250 (MQ=255)                                                                                                                                                                                                          gg‑TTAAGCATCCGCTCA‑GTAACATAA‑CGGCAACCAAAAGAAATAATCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGAG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCTAACCCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CCCGTTCACCTGCAATTGACCGTCGCGTTTGGCAATATCGACAATATAATCGTTGGTGCCCATGCGGCTCGGGTTACTGGTCTTAATATTCTACAGt                                             >  2:391621/1‑249 (MQ=255)                                                                                                                                                                                                           g‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGc                                                                                                                                                                                                                                  <  1:650539/73‑1 (MQ=255)                                                                                                                                                                                                           g‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGc                                                                                                                                                                                                                                  >  2:650539/1‑73 (MQ=255)                                                                                                                                                                                                                aGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTTAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTCGACAGTAGCCCAc                                      <  2:258304/251‑1 (MQ=255)                                                                                                                                                                                                                 gCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTCGACAGTAGCCCAcc                                     <  1:302380/251‑1 (MQ=255)                                                                                                                                                                                                                                           aaCCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTCGACAGTAGCCCACCTAAAATTTGTGACTGGATGGCCTg             >  2:720014/1‑251 (MQ=255)                                                                                                                                                                                                                                            aCCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑CTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGTCAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTCGACAGTAGCCCACCTAAAATTTGTGACTGGATGGCCTgg            >  2:194782/1‑250 (MQ=255)                                                                                                                                                                                                                                                      aTAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTCGACAGTAGCCCACCTAAAATTTGTGACTGGATGGCCTGAGTAAAGTTAt  >  1:440480/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                  CGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGG‑TTAAGCATCCGCTCA‑GTAACATGA‑CGGCAACCAAAAGAAATAAGCGCTGC‑ATGATT‑ATTTTCCTTATGAAGCTGGGG‑CTTAAAAA‑TCGGTTGAG‑TTATTTTGTAAACCCGAAACCT‑GGATGGTCGAGG‑TTTGTCCGGTTTTACGATCTGT‑CACGTTCAACTGCAATTGACCATCGCGGTTGGCAATATCGACAATATAATCGTTGGTCACCATGCGGCCCGGTTTACCGGTATTAATATTCGACAGTAGCCCACCTAAAATTTGTGACTGGATGGCCTGAGTAAAGTTAT  >  NC_000913/1101430‑1101922

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