Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,629,281 G→A I414I (ATC→ATT 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,2810GA90.0% 64.0 / 1.4 30I414I (ATC→ATTettAenergy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (19/8);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

GGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGAT  >  NC_000913/4629051‑4629519
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
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)
                                                                                                                                           gTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTa                                                                                 <  1:325301/251‑1 (MQ=255)
                                                                                                                                            ttCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGCCCAGAGATCATAAGGAACAGGGGCGCTTTa                                                                                 >  2:127821/1‑250 (MQ=255)
                                                                                                                                             tCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTAc                                                                                <  1:306818/250‑1 (MQ=255)
                                                                                                                                                       gTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTc                                                                     >  1:379901/1‑251 (MQ=255)
                                                                                                                                                              cTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAgtttcaccgtt                                                                                                                                                      <  1:169199/163‑1 (MQ=255)
                                                                                                                                                              cTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAgtttcacagtt                                                                                                                                                      >  2:169199/1‑163 (MQ=255)
                                                                                                                                                               tGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTCAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGCCCgat                                                             >  1:202666/1‑251 (MQ=255)
                                                                                                                                                                gATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGATGCCGATCTTCCTAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGCCCAGAGAACATCCGGCACAGGGTCGATTTACCACCGCCGTTCGGACCgatt                                                            >  2:180523/1‑250 (MQ=255)
                                                                                                                                                                  tCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATACCAGTTTCACCGTTTCACCTAAAGTGATGGAGCCGCTGTCCGGCTGTTCCTGACCAGAGAACATAAGGAACAGGGTCGATTTACCCGCGCCGTTAGCACAgatgat                                                          >  1:247337/1‑251 (MQ=255)
                                                                                                                                                                        ccTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTTCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAAACGATGCCAGTTTCACCGTTTCACCTAAAGTGCTGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATAAGGAACAGGGGCGCTTTACCCGCGCCGATCGGAccagtgttcccga                                                     >  2:212583/1‑239 (MQ=255)
                                                                                                                                                                              aagttaaagCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGCTTTACCCGCGCCGTTCGGACCGATGATCCCGACGAttgc                                              >  1:242762/1‑248 (MQ=255)
                                                                                                                                                                                  taaagCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGAGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGAttgcccct                                          >  1:4543/1‑244 (MQ=255)
                                                                                                                                                                                   aaagCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGCGATCATACGTAACAGGGTCGATTTACCCGCGCCGTTCGGCCCGATGAACCcgacaattgcccctt                                         >  2:66881/1‑240 (MQ=255)
                                                                                                                                                                                                  ggCGCGGCTTGGCATCTCGGTGTTGCCGGTCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGCCGATTGCCCCTTTCGGGTTCGACAAGc                          >  2:60507/1‑251 (MQ=255)
                                                                                                                                                                                                            ggCATCTCGGTGTTGCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGtcatc                <  1:60507/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCCGATCTTCATAATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGATGCCAGTTTCACCGTTTCACCTAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCGCCGTTCGGACCGATGATCCCGACGATTGCCCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGGCGAt  <  1:127821/251‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
GGGAACTCCAGCAGGGCGTTTTCCAGCGCGCGCAGGGTTTCGATATCCAGGTCGTTGGTTGGTTCGTCGAGCAGCAGCATGTTGCCGCCAACCTGCAGCAGCTTCGCCAGATGCAGACGACCGCGCTCACCACCGGAGAGTTCACCAACGCGTTTACCCTGATCAACCCCTTTAAAGTTAAAGCGGCCAACGTAGGCGCGGCTTGGCATCTCGGTGTTGCCGATCTTCATGATATCCAGCCCGCCGGAAACTTCTTCCCAAACGGTTTTGCTGTTATCCATTGAGTCACGGAACTGATCAACCGACGCCAGTTTCACCGTTTCACCCAAAGTGATGGTGCCGCTGTCCGGCTGTTCCTGACCAGAGATCATACGGAACAGGGTCGATTTACCCGCACCGTTCGGACCGATGATCCCGACGATCGCTCCTTTCGGGATCGAGAAGCTCAGGTCATCAATCAGCAGACGAT  >  NC_000913/4629051‑4629519

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