Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,629,188 G→A L445L (CTC→CTT)  ettA ← energy‑dependent translational throttle protein EttA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,629,1880GA92.6% 69.2 / ‑3.1 28L445L (CTC→CTT) ettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (12/13);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

AAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGAC  >  NC_000913/4628943‑4629411
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               
aaGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTc                                                                                                                                                                                                                            >  1:369634/1‑251 (MQ=255)
          cGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGt                                                                                                                                                                                                                  <  1:20357/251‑1 (MQ=255)
             agaaCTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCTGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTa                                                                                                                                                                                                               <  1:361071/251‑1 (MQ=255)
                           ccTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAAcccc                                                                                                                                                                                                 >  1:376296/1‑251 (MQ=255)
                            cgtCATCCTGGTGGTTCAGGATGGGCGTGGCGATATGTTCGAGGAACCAACTGTCGTGAGAGATAACCATATCTCAGCACGGGATTGCCAGCAGGTCGGTTTCCAGCGCGAGCAGGGTTTCGATATCATTGTCGTTGTTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCAAGATGCAGACGACCGCGCTCAACTCCGGAAAGTGCACCAACGCGTTTACCCTGATCTACCCCt                                                                                                                                                                                                <  2:201956/249‑1 (MQ=255)
                                   cTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTttaaagt                                                                                                                                                                                         >  1:224877/1‑251 (MQ=255)
                                    tGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCACACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGCTCAACCCCTttaaagtt                                                                                                                                                                                        >  2:157081/1‑251 (MQ=255)
                                        agtCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCACACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTttaaagttaaag                                                                                                                                                                                    <  1:157081/249‑1 (MQ=255)
                                         gtccacggggtgcgtgcCGATAAGGTTGAGGAACCAACGGTCGGGCTAGAGTACCATCTCACAGCCCGGGTACTCCAGCAGTGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTGCGTCGAGCAGCAGCATGTTGCCGCCAACTTGCAGCAGCTGCGCCAGATGCAGACGACAGCGCTCACCTCCGGAAAGTTCACCAATGCGTTGACCCTGATCAACCCCTTGAAAGTTAAAGc                                                                                                                                                                                   <  2:186645/242‑1 (MQ=255)
                                          gccatgaggtgctgggCGAGACGGTCGAGGAACCAACGGGCGTGCGAGATAACCAGCGCACAGCCCGGGAACTCTAGCAGGGCTGTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTGCGTCGAGCATCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCTCCTCCGGAAAGTTCACCAACTCGTTTTCCCTGATCAACCCCTTTAAAGTTAAAGCg                                                                                                                                                                                  <  2:376296/237‑1 (MQ=255)
                                                    gCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAg                                                                                                                                                                        <  1:368715/251‑1 (MQ=255)
                                                    gCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAg                                                                                                                                                                        <  1:128954/251‑1 (MQ=255)
                                                                   cGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCAtc                                                                                                                                                         <  1:200900/251‑1 (MQ=255)
                                                                         aCCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGgtg                                                                                                                                                   <  1:70396/251‑1 (MQ=255)
                                                                                                            gggAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGCTTCGTCGAGCAGCAGCATGTTCCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGCCCGCGCTCACCACCGGACAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAACGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCCTCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGacac                                                                                                                >  1:225991/1‑248 (MQ=255)
                                                                                                                    cagcagGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGAAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGAAAAGTTCACCAACGCGTTTACCCTGCTAAACCCCTTTAACGTTAAAGAGGCACACGTAGGCGCGGCTTGGCATATCGGTGGTGCCGATATCAATAAAATCCAGCCCGACCGAAActtcttgcc                                                                                                        >  2:383717/1‑248 (MQ=255)
                                                                                                                         tggCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGGTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACACTGATCAACCCCTTTAAAGTTAAAGCGGCAAACGTAGGCGCGGCTTGGCCTCCAGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGATACTTCCTCCCAAACg                                                                                                   >  2:158035/2‑251 (MQ=255)
                                                                                                                                  ccAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTg                                                                                          >  1:378355/1‑251 (MQ=255)
                                                                                                                                                                cGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCCACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTAACCAACGCGTTTACCCTGATCAACCCCTTTCAAGTTACAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGGTGCCAGTCTTCATAATATAAATCCCGACGGAAACTTCTTCCCAAACGGTTTTCCTCTTATCCATTGAGTAACgaacatgatcaccc                                                            >  2:323311/1‑238 (MQ=255)
                                                                                                                                                                cGTTGGGTGTTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAAcc                                                            <  1:215701/251‑1 (MQ=255)
                                                                                                                                                                       tggttCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCACACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGcc                                                      >  1:36452/1‑250 (MQ=255)
                                                                                                                                                                                gAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcaccg                                            <  1:379496/251‑1 (MQ=255)
                                                                                                                                                                                   cagcagcaTGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAgtttcaccgttt                                         >  1:88048/1‑251 (MQ=255)
                                                                                                                                                                                    agcagcaTGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAACGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGCCACTTCTTCCCAAACGGTTTTGCTGTTATCCCTTGAGTCACGGAACTGATCAACCAATGCCAgtttcaccgtttc                                        >  2:66263/1‑251 (MQ=255)
                                                                                                                                                                                         caTGTTGCCGCCAGCCTGCAGCAGCTGCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTGCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCa                                   <  1:180798/251‑1 (MQ=255)
                                                                                                                                                                                                       ccTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCg                     <  1:66263/251‑1 (MQ=255)
                                                                                                                                                                                                                    cGCCAGATGCAGACGACCGCGCTCACCACCGGAAAGATCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAAATGATCAACCGATGCCAGTTTAACCGTTTCAACTAAAGTGATGGTGCCGCTGTCCGGCTGtt        >  1:69615/1‑251 (MQ=255)
                                                                                                                                                                                                                          attcagtcgacCGCGCTCACCACCGGAAAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCTGCTTGGCTTCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGTTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGAc  <  1:368042/244‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                               
AAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATCGCACAGCCCGGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGAC  >  NC_000913/4628943‑4629411

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