Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	458,189	0	A	G	31.6%	56.3 / 10.1	40	D255G (GAT→GGT)	clpX	ATP‑dependent Clp protease ATP‑binding subunit ClpX
Reads supporting (aligned to +/- strand):  ref base A (5/21);  new base G (12/0);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 2.29e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCA  >  NC_000913/457946‑458432                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         aGAAAGCACAGCGTGGTATTGTCTACATCGATGACATCGACAAGATTTCTCGTACGTCAGACAACCCGTCCCTTACCCGAGCCGTTTCCGCTGACGGCGTACAGCAGGCACTGTTGAAAATGATCGAAG‑GTACGGTAACTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGTTTTCCCGTCTGGATAAAggg                                                                                                                                                                                                                                                   >  1:222434/1‑249 (MQ=255)             gTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTGCGGTCTGGGTAAAGGGGT‑TTCCCACCgt                                                                                                                                                                                                                                      >  1:79957/1‑251 (MQ=255)                   ttGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGCTGATACCTCTAAGATCCTGTTTATTTGTGGCGGAGGGGTTGGCGGGCTGGTTAAAGTGGT‑TTTCCACCGGGTAGa                                                                                                                                                                                                                                 >  1:43126/1‑250 (MQ=255)                         aCATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCg                                                                                                                                                                                                                            <  1:46664/249‑1 (MQ=255)                               aTGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCCGTCTGGATAAAGTGAT‑TTCCCCCCCGGTAGAAACCGGCTCCGGc                                                                                                                                                                                                                    >  1:44989/1‑251 (MQ=255)                                   aaTCGATAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTATGGGTTAAGTGGT‑TTTCCCCCCTGTAGGAACCCGCTCCCGCAttt                                                                                                                                                                                                                >  1:148793/1‑250 (MQ=255)                                   aaTCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCATGTAGAAACCGGCTCCGGCATTg                                                                                                                                                                                                                <  1:129066/251‑1 (MQ=255)                                               ttCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGcgcg                                                                                                                                                                                                    <  2:226534/250‑1 (MQ=255)                                                     tAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGTGAT‑TTTCCCCCGGGTAGAAACCGGCTCCGGCATTGG‑GTTTGGCGCG‑ACGGTa                                                                                                                                                                                             >  1:334835/1‑250 (MQ=255)                                                                 cccGTCAATTACCCGAGACGTTTCCGGTGCAGGCGTACAGCAGGCACTGTTGAACCTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGGTTGGGTAAAGAGAT‑TTCCCACCGTGGAGAAAAAGGCTCCGGGATTGG‑GTTTGGGGCG‑GCGGGAAAAGGGAAggacg                                                                                                                                                                                >  2:66834/1‑247 (MQ=255)                                                                  ccGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGa                                                                                                                                                                               <  2:169601/251‑1 (MQ=255)                                                                            cccGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCGGGCTGGGTAAAAAGAT‑TTTCCCCCGTGTAGAAACCCGCTCCGGCATTGG‑TTTTTGCGCG‑ACGGTAAAAGCGAAGGCCGACAaagcaagc                                                                                                                                                                     >  1:489918/1‑251 (MQ=255)                                                                               gagaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCg                                                                                                                                                                                       <  2:205788/230‑1 (MQ=255)                                                                                                           ggCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTTGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGGGCCGGTTGAACCgg                                                                                                                                      >  1:114290/1‑251 (MQ=255)                                                                                                                    taaaCTGATCGAAGTGT‑CGGTTGCTGCTGTTCCACCGCTAGGTGGGCGTAATCATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTGTATTTGTTGCGGTGCGTTTGCCGGTCTTGATAAAGTGAT‑TTCCCACCGTGTATAAACCGGCTCCGGCATGGG‑TTTTGGCGCT‑ACGGTAAAAGCGAAGTCCGACAATGCAAGCGAAGGCGAGCTGGTGGCGCAGTTTGAACCGGAAGATCTg                                                                                                                              <  2:222434/249‑1 (MQ=255)                                                                                                                     aaaCTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTGGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCTGCATTGG‑TTTGGGCGTG‑ACGGTGAAAGCGAGTTCCGACAAAGCAAGCGAAGGCGAGCGGCTGGCGCAGGTTGAATGGGATGATCTG‑At                                                                                                                           <  2:309885/251‑1 (MQ=255)                                                                                                                      aaCTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATc                                                                                                                          <  1:160752/251‑1 (MQ=255)                                                                                                                              gAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGTCGGTCTGGGTAAAGAGAT‑TTTCCCCCGTGTAGAAAACGGCTCCGGCCTTGGGTTTTGG‑GCG‑GCGGGCAAAGGGAAGTCCGACAAAGCAAGCGAAGGGGGGCTGCTGGCGGAGGGTGAACCGGAAGATCTG‑ATCAAGATTgg                                                                                                                  >  2:24473/1‑251 (MQ=255)                                                                                                                               aaG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGTAGGAATTCTTGCAGGTTGATTCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGt                                                                                                                 <  1:286881/251‑1 (MQ=255)                                                                                                                                 g‑gTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCCGGTCTGGGTAAAGGGAT‑TTTCCACCGTGTAGAAACCGGCCCCGGCATTTG‑GTTTGGCGCG‑GCGGGTAAAGCGAAGGCCCACAAAGCAAGAGGAGGCGAGCTGCTGGGGCCGGGTGGAACGGGAGATCCG‑GTCAAGTTTTGTCt                                                                                                               >  1:187445/1‑251 (MQ=255)                                                                                                                                 g‑gTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGTGAT‑TTCCCCCCGTGTAGAAACCGGCCCCGGCATTGG‑GTTTTGCGCG‑ACGGTAAAAGCGAAGGCCGACAAAGCCAGCGAAGGCGAGCGGCTGGCGGCCGTTGAACCGGGAGATCTG‑ATCAAGGTTgtgct                                                                                                               >  1:155941/1‑247 (MQ=255)                                                                                                                                                          caaagTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGGTCGTTTGCCGGTCTGGGGACAGGGATGTTCCCAC‑GTGTAGAAAACGGCTCCGGCCTTGG‑GTTTGGCGCGGAGGGAAAAAG‑GAAGTCCGACAAAGCCAGCGAAGGCGAGCCGCGGGGCCAGGTTTAAACGGGAGAGCTGGACAAAGTTTGG‑CTTATGCCTGAGTTTATTgggagtgt                                                                                       >  2:282190/5‑245 (MQ=255)                                                                                                                                                            aGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGc                                                                                     <  1:499029/251‑1 (MQ=255)                                                                                                                                                              gTGGGCGTAAACAGCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCg                                                                                   <  2:423492/251‑1 (MQ=255)                                                                                                                                                                        aCATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGCGCGTTTGTCGGTCTGGGTAAAATGAT‑TTCCCCACGTGTAGAAAACGGCTCCCGCCTTGG‑GTTTGGGGCG‑GCGGTAAAAGCGAAGTCCGGCAAAACAAACGAAGGCGAGCTGGTGGGGCAGGGTGGAACGGGAGATCTG‑ATCAAATTTTGGCTTATCCCTGGGTTTATTGGTCGTCTGTCGGTTGTCG‑cc                                                                         >  2:238488/1‑249 (MQ=255)                                                                                                                                                                             cGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTtg                                                                   <  1:24473/251‑1 (MQ=255)                                                                                                                                                                                       aTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGa                                                         <  1:114189/251‑1 (MQ=255)                                                                                                                                                                                          cTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑GCGGTAAAAGCGAAGTCCGACAAAGCAAGCGGAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTTGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCg                                                      >  2:91113/1‑251 (MQ=255)                                                                                                                                                                                             gCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCgaag                                                   <  1:66834/251‑1 (MQ=255)                                                                                                                                                                                                ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGAGTTTGCCGGGACGGCTAAAGTGAT‑TTTCCACCGTGTAGAAAACGGCTCCGGCATTTG‑TTTTGGCGCG‑ACGGGAAAAACGAAGGCCGCCAAAACAAGCGGAGGCGGGGTGGTGGGGCAGGTTGAACCGGGAAATCTG‑ATCCAGGTTGGTCTTATCCCTGCGTTTATTTGTAGTCTGCCGGTTGTCG‑CCACGTTGGATGGAATGAGCgcaga                                                  >  2:137168/1‑249 (MQ=255)                                                                                                                                                                                                    gattccTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCGGGATAAAGTGAT‑TTCCCACCGGGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGGCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCGAAGAAGctct                                            <  2:187445/247‑1 (MQ=255)                                                                                                                                                                                                     aTACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCGAAGAAGCTCTg                                           <  1:6793/251‑1 (MQ=255)                                                                                                                                                                                                      tACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGGCTGGGTAAAGTGAT‑TTTCCACCGGGTAGAAACCGGCGCCGGCCTTGG‑TTTTGGAGCG‑ACGGTAAAAACAGAGTCCGACAAAGCAAACGAAAGGGAGGCGCTGGCGCCGGGTGCAACGGGAGGTCGG‑ATCAAATTTGGTTTTTTCCCTGTGTTTTTTGGGCGTCTAACGGATGTCGGCAAC‑Ttgaatggaaagtgcgaaaaaactctg                                           >  2:422245/1‑230 (MQ=255)                                                                                                                                                                                                                     gTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGGT‑TTCCCACCGTGTAGAAACCGGCTCCGGCCTTGG‑TTTTGGCGCG‑GCGGTAAAAGCGAAGTACGGCCAAGCAAGCGAAGGCGGGATGCTGGCGCAGGCTGAACCGGAAGATCTG‑ATCAAGTTTGGACTTATCCCCGAGGTTTTTGGTCGTCTTCCGGGTGTCG‑CAAAGTTTAATGGAATGAGGGAAGAAGCTCTGATTTAGATCCTcca                             >  1:512661/1‑247 (MQ=255)                                                                                                                                                                                                                       ttATTTGTGGCGGGGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTACGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑AGCAAGTTTGGGCTTATCCCTGAGTTTATTGGTCGTCTGCCGGtt                                                                                <  1:530682/197‑1 (MQ=255)                                                                                                                                                                                                                              tgGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGCGCGAAGAAGCTCTGATTCAGATCCTCAACGAGCCGaaa                   >  2:187029/1‑250 (MQ=255)                                                                                                                                                                                                                                   gTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGcc             <  1:238488/251‑1 (MQ=255)                                                                                                                                                                                                                                     gCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGccc            <  1:91113/250‑1 (MQ=255)                                                                                                                                                                                                                                      gggttGCCGGTCTGGATAAAGTGAT‑TTCCCACCGGGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAAGGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTg          <  2:44989/248‑1 (MQ=255)                                                                                                                                                                                                                                              ggTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCa  <  1:405308/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         AGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAG‑GTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGAT‑TTCCCACCGTGTAGAAACCGGCTCCGGCATTGG‑TTTTGGCGCG‑ACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTG‑ATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCG‑CAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCA  >  NC_000913/457946‑458432

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