Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,268,876	0	A	T	27.8%	102.0 / 15.1	54	I238F (ATT→TTT)	kdsA	3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (8/31);  new base T (15/0);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 5.67e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGG  >  NC_000913/1268640‑1269126                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                   gcAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCatccg                                                                                                                                                                                                                                                    <  1:410474/251‑1 (MQ=255)   aaCGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGGAGGTCTGGCGGGGCTG‑TTTTTT‑TT‑AGCGCatcccgg                                                                                                                                                                                                                                                  >  2:170532/1‑248 (MQ=255)       aaaaaGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTG‑TTTTTT‑TG‑AGCGCctccgggaccc                                                                                                                                                                                                                                              >  1:461967/1‑246 (MQ=255)       aaaaaGTGATGCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCatccggatccg                                                                                                                                                                                                                                              <  1:312687/251‑1 (MQ=255)              gATTCTTTGCGATCGCGGTGCTAACTTCGGCGATGACAACCGGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCGGGT‑AACCCGCCGGTGATTTTCGACGCGACCCGCGCAC‑GGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGGCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑GTTATT‑GA‑AGCGCATCCGGATCCGGAACATg                                                                                                                                                                                                                                       <  2:143685/251‑1 (MQ=255)                 tCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCg                                                                                                                                                                                                                                     <  1:227729/250‑1 (MQ=255)                         aTCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑GTTTTT‑GT‑AGCGCATCCCGATCCGGAAAATGCG‑AAAAGTGAt                                                                                                                                                                                                                           >  1:565112/1‑251 (MQ=255)                            gcgGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTGGTTTTTT‑GA‑‑GCGCATCCGGGTCCGGGAAATGCG‑GAAAGTGATggg                                                                                                                                                                                                                        >  1:25187/1‑250 (MQ=255)                                           ggtatataacacTGGTTGTAGCTATGTTGGGCTTAAGCATTATGAAGAAAGTG‑CTGGTTAATTGGCGGGTGATTTTCGCTGTGACCCACGCCCGTG‑AATGCCGCGATTCGT‑TTGGCGCCGCTGCCGGTGGGTGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑TGAGCATCCGGATCCGGAAAATGCG‑AAATGTGAGGGTCCATCCGCGCTGCCg                                                                                                                                                                                                         <  2:412275/240‑1 (MQ=255)                                                   aaCCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCCG‑TTTTTT‑TG‑AGCGCACCCCGGTCCGGAACATGCG‑GAATGTGATGGTCCCTCCCCGCCGCCGCCGGcgt                                                                                                                                                                                                  >  2:148215/1‑248 (MQ=255)                                                     ccTGGTTGTCGATATTCTGGGCCTCAGCATTATGAAGAATGTGTTTGGT‑AACTCGCCGGTGATTTTCGGCGTGCCCCACGCCC‑TGCAGTGCCGCGATCCGT‑TTTGCGCAGCTTCCGGTGGTCTTGGTGCCCAGGTGGCTGAGCGGGCCCGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCGGCCGCTGGCTAAA‑c                                                                                                                                                                                              <  2:505727/251‑1 (MQ=255)                                                     ccTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTTTT‑GA‑AGCGCCTCCGGGTCCGGGACATTGG‑AAATGTGATGGTCCAACCCCGCCGCCCCCGGCTaaa                                                                                                                                                                                                >  1:328060/1‑250 (MQ=255)                                                         gTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTTTT‑TA‑AGCGGCTCCGGGTCCGGAAAATGCG‑GAAAGTGGTGGTACATCCGCgccgcccccggccaaaatggg                                                                                                                                                                                           >  1:557147/1‑236 (MQ=255)                                                            ttcggtatGCTGGCCTTAAGCTTTATAAAGAAATGGTTGGGT‑AACTCCCCGGGGATTTTCGTCGGGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TGGGCGCAGCTTCGGGTGGCCGTCGTGCTCAGGTGGCTGAGCTGCCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGTCTCTGGCTAAA‑CTGGAAc                                                                                                                                                                                        <  2:388885/245‑1 (MQ=255)                                                                              aGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGaa                                                                                                                                                                    <  1:374716/251‑1 (MQ=255)                                                                               gCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTG‑GTTATT‑TT‑AACGCCTCCCGGTCCGGGACAAGCGGAAAT‑AGATGGTGCCTCCGCGCTGCCGCTTGCTAAA‑CTGGAACCG‑TTTCTCAAACCGATGaa                                                                                                                                                                    >  1:566279/1‑250 (MQ=255)                                                                                 ataaTGAGGAAAGTGTCGGGT‑AGCACGACGGTGATTTTCGGCTTGCCCCCCGCGC‑TGCAACGCCGCAGTCGGT‑TTGGCGCAGCTCCCGGTGGTCTTCGAGCTCAGGTGGATGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGc                                                                                                                                                                 <  1:301831/248‑1 (MQ=255)                                                                                    aTGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTTTT‑GA‑AACGCCTCCGGGTCCGGAACATGCG‑AAAAGTGATGGGCCCTCCGGGCTGCCGCTGGCTAAAACTGGAACCG‑T‑CCTCaaaccgaataaaggggt                                                                                                                                                              >  1:61718/1‑236 (MQ=255)                                                                                     ggaagaaAGGGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGCCCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGAtt                                                                                                                                                             <  1:170532/250‑1 (MQ=255)                                                                                        agaaAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTG‑GTTATTTGAGAGCCCA‑CCG‑ATCCGGAACAAGCG‑GAAAGTGGTGGTCCATTCGCGCTGGCGCCGGCTAAA‑CTGGAAACG‑GTCCTCCAGCAGATGAAAACGATtgat                                                                                                                                                          >  2:69213/1‑251 (MQ=255)                                                                                             gtgtCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTG‑TTTATT‑TA‑AGCGCATCCCGATCCGGGACATGCG‑AAATGTGATGGGCCATCCGCGCCGCCGCCGGCTAAA‑CTGGAACCG‑TTCCTCCAACAGATGAAAGCGATTGATGGTCt                                                                                                                                                     >  1:355337/1‑251 (MQ=255)                                                                                                  tGGT‑AACTCGACGGGGATTTTCGACGTGACCCACGCAC‑TGCATGGCGGCGTGCGGT‑TTGGCGCAGCGTCCGGTGGTCGTCGTGCTCAGGGGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTg                                                                                                                                                 <  2:4410/250‑1 (MQ=255)                                                                                                           cGCCGGTGATTTTCGAGGCGACCCCCGCAC‑TGCAATGCCGCTATCAGTGTTGG‑ACAGCTGCCGGTGGCCGTCGTGCTCGGGTTGCTGAGCTGGCACGACCAGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTGATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCGGGCTAAA‑CTGGAACCG‑TTCCACAAGCATATGAAAGCGATTGATGATCTGGTGAAAGGtt                                                                                                                                          <  2:338547/249‑1 (MQ=255)                                                                                                                    ttttcgtcctgcccctcGCAC‑TGCAATGCGGGTATCTGT‑TGGGCGCATCTTCCGGTGTCCGTCGTGCTCGTGTGGCAGAGCTGGCACGAGCCGGTATGGCGTTCGGTCTGGCGGGGCTG‑TTTAGT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGGTCTGGTGAAAGGTTTTTAAGAACt                                                                                                                                <  2:216806/234‑1 (MQ=255)                                                                                                                          aCGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTG‑TTTTTT‑TA‑AACGCCTCCGGGTCCGGGAAATGCG‑AAATGTGATGGTCCATTCGCGCCGCCGCTGGCTAAA‑ATGGAACCG‑GTCCTCAAGCAGATGAAAGCGGTTGATGATCTGGGGGAAGGGTTCGGAGAAACGGATAcc                                                                                                                         >  1:323094/1‑251 (MQ=255)                                                                                                                             ggACCCACGAAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑c                                                                                                                     <  1:587450/250‑1 (MQ=255)                                                                                                                                      cAC‑TGCAAGGCAGCGATCCGT‑TTGGCGCAGGTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCt                                                                                                            <  2:328060/251‑1 (MQ=255)                                                                                                                                       aC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTG‑TTTTTT‑GG‑AGCGGATCCGGGTCCGGGAAAAGCG‑GAATGGGATGGGCCCTCCGCGCCGGCCCTGGCTAAA‑ACGGAACCG‑TTTCTCAAGAAGAAGGAAGCGGTTGATGATCTGGGGGAAGGGTTCGGAAAACTTGATaacaaaaaggaaactt                                                                                                            >  1:413515/1‑236 (MQ=255)                                                                                                                                             aaTGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTg                                                                                                      <  2:355337/251‑1 (MQ=255)                                                                                                                                                 ccGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCGGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑g                                                                                                 <  2:461967/251‑1 (MQ=255)                                                                                                                                                 ccGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTG‑TTTTTT‑TA‑AGCGCATCCGGGTCCGGGACCTGCG‑AAATGTGATGGTCCATCCCCGCCGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGGGAAAGGGTTCGAAGAAATGGATACCAA‑CAAGGAATCTTTTTTTCttt                                                                                                  >  1:325284/1‑250 (MQ=255)                                                                                                                                                   gcgATCCGT‑TTGGCGCAGCTTCCGGTGGTCTTCGTTCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑Gaa                                                                                               <  1:138098/251‑1 (MQ=255)                                                                                                                                                    cgATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTTTT‑TA‑AGCGCCCCCGGATCCCGGACATGCG‑GAAAGTGATGGGCCAACCGCGCCGCCGCTGGCTAAA‑ATTGGACCG‑TTTCTCCAACAGGTGAAAACGATTGGTGATCTGGGGAAAGGGTTCGAAGGACTGGGTACCAG‑ACAGGAATCTTTtttgttt‑gta                                                                                               >  1:99624/1‑248 (MQ=255)                                                                                                                                                     gATCCGT‑TGGGCGCGTCTTCCGGTGGTCGGCGTGCTCAGGTGGCTGAGTGGGCCCGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTGATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑Gaaaa                                                                                             <  2:118715/251‑1 (MQ=255)                                                                                                                                                               ggCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTAt                                                                                    <  1:238210/251‑1 (MQ=255)                                                                                                                                                                  gcAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGGTG‑GTTATT‑GA‑AGGGCATCCGGGACCGGGAAATGCG‑GAATGTGATGGGCCATCCGCGCCGCCGCTGGCTAAA‑CTGGGACCG‑TTTCTCAAGCAGATGAAAGCGATTGATGATCTGGGGGAAAGTTTTGAAGAAATGGATAACAG‑CCAGGAATCTTTTTTTGTT‑GAAAAATAAAGTATTa                                                                                  >  1:12890/1‑250 (MQ=255)                                                                                                                                                                         cgggtggtCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCt                                                                           <  2:325284/248‑1 (MQ=255)                                                                                                                                                                              ggtCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTa                                                                    <  2:58784/251‑1 (MQ=255)                                                                                                                                                                              ggtCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTa                                                                    <  1:259464/251‑1 (MQ=255)                                                                                                                                                                                       gCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACttttt                                                           <  1:69213/251‑1 (MQ=255)                                                                                                                                                                                          cAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACtttttt                                                          <  1:148215/249‑1 (MQ=255)                                                                                                                                                                                               ggCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGt                                                   <  2:323094/251‑1 (MQ=255)                                                                                                                                                                                                gCTGCGCTGGCACGAGCCGGTATGGCGGGAGGTCTGGGGGGGGCG‑TTTTTT‑GA‑AGCGCATCCGGAGCCGGGACCAGCG‑GAATGGGATGGTTCATCCGCGCTGCCGCTTGCTATA‑CTGTGACCG‑GTTCTCCAGGAGGTGGAAACGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATAACAGACAAG‑AATATTTTTTGCTT‑TAAAAAAAAAAGATTAACGTTCTgtggttagatttttttgtagggtg                                                   >  1:80955/1‑228 (MQ=255)                                                                                                                                                                                                 cggcgctggccagAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGc                                                 <  2:73700/239‑1 (MQ=255)                                                                                                                                                                                                 cTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTTTT‑GA‑AACGCATCCGGATCCGGAACAAGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CCGGGACCGGTTCC‑CCAGCAGATGAAAACGATTGGTGATTTGGGGGAAGGGTTCGCAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑TAAAAATAAAGAATTAACGTTCTGCG‑GTAAGACCTTTTTTATgggggc                                                 >  1:491031/1‑248 (MQ=255)                                                                                                                                                                                                   gAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTCTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCg                                               <  2:25187/251‑1 (MQ=255)                                                                                                                                                                                                    ggCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCgg                                              <  1:190639/250‑1 (MQ=255)                                                                                                                                                                                                    aGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTGGTTTTTT‑GA‑‑GCGCCTCCCGGTCCGGGACATTGG‑GAATGTGGTGGGCCATCCCCGCTGCCGCCGGGCAAA‑ATGGGACCG‑GTTCTCCAGCAGATTAAAACGGTTGGTGATCTGGGGGAAGGGTTCGAAGAAACGGGTACCAG‑CAAGTAAACTTTTTTTCTTTGAAAAA‑AAAGGATTAGCGTTTTTCGGTTAAAATTTTTTTC‑TGGGTGCCgg                                              >  2:66755/1‑251 (MQ=255)                                                                                                                                                                                                                                   ggCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCt                <  2:565112/250‑1 (MQ=255)                                                                                                                                                                                                                                       tggCTG‑TTTATT‑GA‑AGAGTATCCGGATCCGGAACATGAG‑AAATGTGATTGTCTATCAGTGCTGACGCTGGATAAA‑CTGGAACCG‑TTCATCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTAGAAGAACTGGATACCAG‑AAAGTAATATTTTTTGGTT‑GAAAAATGAAGTATTAGCGTTCTGAA‑TTAAGAATTTTTTCATGGGTGCAGGATACAAAAAAAGCCGCAGGCTGTTActattgcggg           >  2:84682/2‑242 (MQ=255)                                                                                                                                                                                                                                         gCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCg         <  1:430924/251‑1 (MQ=255)                                                                                                                                                                                                                                          cTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCgg        >  1:423150/1‑251 (MQ=255)                                                                                                                                                                                                                                            g‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑CTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGtt      >  1:347570/1‑251 (MQ=255)                                                                                                                                                                                                                                               ttATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGATTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTc    >  2:129988/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aTT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑CAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCgg  >  2:425946/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                   GCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGT‑AACTCGCCGGTGATTTTCGACGTGACCCACGCAC‑TGCAATGCCGCGATCCGT‑TTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCACGAGCCGGTATGGCGGTAGGTCTGGCGGGGCTG‑TTTATT‑GA‑AGCGCATCCGGATCCGGAACATGCG‑AAATGTGATGGTCCATCCGCGCTGCCGCTGGCTAAA‑CTGGAACCG‑TTCCTCAAGCAGATGAAAGCGATTGATGATCTGGTGAAAGGTTTCGAAGAACTGGATACCAG‑CAAGTAATCTTTTTTGCTT‑GAAAAATAAAGTATTAGCGTTCTGCG‑TTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGG  >  NC_000913/1268640‑1269126

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