Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,248,4240CA37.9% 35.2 / 10.8 29V37F (GTT→TTT) lysPlysine:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (16/2);  new base A (0/11);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 2.25e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAATCTCTTTTCTGG  >  NC_000913/2248174‑2248670
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
aCCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGcaac                                                                                                                                                                                                                                                        <  1:408935/251‑1 (MQ=255)
       gagctgCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAg                                                                                                                                                                                                                                                  <  1:86042/250‑1 (MQ=255)
                       aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAAGCAACAAAAAGACCTGTACCGATGGAAc                                                                                                                                                                                                                                 >  1:282288/1‑251 (MQ=255)
                         cGGCGATAGTCCCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTGCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCACCCGTACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCg                                                                                                                                                                                                                               <  2:57418/251‑1 (MQ=255)
                             gATAGTCACCGCCCAGTTGTCCCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCACTCAGCATATACG‑AGAGCAATGCCACGCCCGGCCCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACACAAAGACCTGTACCGATGGAACCGCCa                                                                                                                                                                                                                            >  2:461288/1‑250 (MQ=255)
                             gATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCa                                                                                                                                                                                                                            >  1:254469/1‑250 (MQ=255)
                                         ccAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCg                                                                                                                                                                                                               >  1:399280/1‑251 (MQ=255)
                                                  ccAGTAGTCCCAGCCCGGGGCGAAGCAAAAGCCTTCTTAAACATATTTCTGACCGAAAGTGGAAAACAAACAGGAACCCGGAATATAAGCAGCCAGTCCACCGAGACTGGACA‑TCAGGAAGAAAACCAGCGGGCAAAACAGAATATGCG‑AGAGAAACCCCCCCCCCGGCCCTGCAGGA‑AAAACCGTGGCGCCAGAGACAAAAAAAAGCCCTGTACCGGGGGAACCGCCAATGGCAATCATCGTCAGGTGAc                                                                                                                                                                                                      <  2:78898/251‑1 (MQ=255)
                                                   cAGTAGTCCCGGCCCCGCGCGAAGCAAAACCTTTCTCAAACATATTTCTGCCCGAAAGGGGAAAACGAACCGAAACCCGGAATATAGGCAGCCAGTTCCCCGGGACTTGTCA‑TCAGGAAGTAACCCACCAGGCCAACCAGCATATAAG‑AGAGAAAGGCCCCGCCCGGCCCTGCCGGA‑GAAATCTTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGAcg                                                                                                                                                                                                     <  1:385868/251‑1 (MQ=255)
                                                                            aaaccccTTCTTCAACATAGTTCTGACCGAAAGTGGAAAACGACCCGAAACCCGGAATATATGCAGCCAGTTCACCGAGACGGGTAA‑ACAGAAAGAAAACCATCGGGCCAATCAGCATATACG‑AGAGAAACCCCCCGCCGGGCCCTGCCGGA‑AAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGccc                                                                                                                                                                            <  1:9381/246‑1 (MQ=255)
                                                                              aaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCACACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAGTCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCgg                                                                                                                                                                          >  2:35620/1‑251 (MQ=255)
                                                                               caccTTCTACAACATATTTCTGACCGAAAGAGGAAAACGAGCCGAAACCCGTCAAATATACAGCCAGTTACCCGAGACCGGTCA‑ACAGGAATTAAACCAACAGCCCAACAACCATATACG‑AAAACAAACCCCCCCCCGGACC‑GCCTGAAGAAATCGTTCCGCCGAAGACAAAAAAAAGCCATTTCCAGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                                                                                                                         <  1:37232/249‑1 (MQ=255)
                                                                                   ttcttcAACACAGTACTGACCGTGCGAGGCCAACGAACCGGAAACCGGCATATATGCAGCAAGCTCACCGAGCCTGGTTA‑ACAGGAAGCAAACCATCAGGCCACTCAGCATATACG‑AGAGCAATGCCCCGCCCGGCCCTGCCTGA‑GAACTCGTTACGCCAGAGGCAACCAAAAGACCTGTCCCGCTGGAACCGCCCATGGCACTCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTCAGCCCGGCGCtt                                                                                                                                                                     >  2:287102/1‑251 (MQ=255)
                                                                                       aaaataTAGTTCTACCCGTAAGGGACAAAAGACCCGGAACCCGGAAAATAGGCAGCCATTCCCCCGAGACTGGTAA‑CCAGGAAGTAAACCACCAGGCCAACCAGAAAAGAGG‑AGAACAAGCCCCCGCCCGGCCCTGCCTGA‑AAAATCGTGGCCCCAGAGGCAAAAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                                                                                                                 <  2:73212/246‑1 (MQ=255)
                                                                                       aaaactttgtccTGACCGAAAGTGGAAAACGAACCGAAAACCGGAATATATGAACCAAGTTCACCGAGA‑TGGGCAATCAGAAAGAAACCCATCAGCCCAATCAGAATATACG‑AGAGCAATCCCCCGCCCGGCCCTGCTTAA‑AAATTCGTTGCCCAAGAGGCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                                                                                                                 <  1:244075/240‑1 (MQ=255)
                                                                                               gTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                                                                                                                           >  2:221914/1‑249 (MQ=255)
                                                                                                   tGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTc                                                                                                                                                     >  2:413576/1‑251 (MQ=255)
                                                                                                    gACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCg                                                                                                                                                    >  2:97251/1‑251 (MQ=255)
                                                                                                                aaaacaaacCGAAACCCGACATATAGGCAGCCAGTTCACCGAGATTGGTAA‑AAAGGAAGAAACCCACCAGCCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGAAAAAAAAAGCCCTGTACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaa                                                                                                                                        <  1:379223/250‑1 (MQ=255)
                                                                                                                           gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                                                                                                                             >  2:403169/1‑251 (MQ=255)
                                                                                                                                 ggCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAACTGCTATCCATCtt                                                                                                                        >  2:265464/1‑250 (MQ=255)
                                                                                                                                            aGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcg                                                                                                            >  1:309160/1‑251 (MQ=255)
                                                                                                                                                    cACCGGGAGTGGTCA‑TAAGAAAGAAAACCATCAGGCCAATCAGCATA‑ACGGAGAGCAAGCCCCCGCCCGGCCCTGCCTGA‑AAAATCGTTGCCCCAGAGCCAAAAAAAAGCCCTGTACAGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTa                                                                                                    <  1:413576/251‑1 (MQ=255)
                                                                                                                                                              ggTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCATTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAtt                                                                                          >  2:233733/1‑251 (MQ=255)
                                                                                                                                                                       ggAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGACCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATcc                                                                                  >  2:399595/1‑251 (MQ=255)
                                                                                                                                                                                                        gagcacAATCCCCCGCCCGGACCTGCCGAA‑GAAATCGTGCCGCAAGAGACAAAAAAAAGACCTTTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCt                                                  <  1:153766/246‑1 (MQ=255)
                                                                                                                                                                                                                             ccggccgaa‑gaAACCGTGGCGCAAGAAGCAAAAAAAAGACCTGAACCGAAGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTGATAAGACACCGTCATGATCGCCCAGGGATtat                             <  2:226142/243‑1 (MQ=255)
                                                                                                                                                                                                                                   atgaaaaaTCGTGGCCCCAGAGGAAAAAAAAAGCCGTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTa                        <  1:80856/246‑1 (MQ=255)
                                                                                                                                                                                                                                               tGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAAtc            >  1:247024/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         gCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAATCTCTTTTCTgg  >  1:99261/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      
ACCAGCTGAGCTGCAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCA‑TCAGGAAGTAAACCATCAGGCCAATCAGCATATACG‑AGAGCAATGCCCCGCCCGGACCTGCCTGA‑GAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCAATCCTTCCGGCAGAAATACACGGTTATAAGACACCTTCATGATCGCCCAGGGATTATAAGTAAAGCAGCCAATCTCTTTTCTGG  >  NC_000913/2248174‑2248670

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: ‑