Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,344,935 G→A T208T (ACC→ACT yjdF ← DUF2238 domain‑containing inner membrane protein YjdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,344,9350GA88.0% 54.6 / ‑2.7 26T208T (ACC→ACTyjdFDUF2238 domain‑containing inner membrane protein YjdF
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (13/9);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCGGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTCC  >  NC_000913/4344709‑4345177
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
gctgGATAAATATCGCAAAGTTCCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCTTGATGCGATATGTTTATCAGGCCCCGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCg                                                                                                                                                                                                                                <  1:3119/247‑1 (MQ=255)
          tatCGCAAGGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTa                                                                                                                                                                                                                                                <  2:91312/221‑1 (MQ=255)
          tatCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTa                                                                                                                                                                                                                                                >  1:91312/1‑221 (MQ=255)
                aaaGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAagc                                                                                                                                                                                                            >  1:291537/1‑251 (MQ=255)
                                  gTTCATGCTGAGATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCAct                                                                                                                                                                                           <  1:142029/250‑2 (MQ=255)
                                      aTGCTGATATTCCGGTCGGTGCAGTGAGTGATGTAAAAGCCTGACGTGGACAATAGCTAAAGATAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGAGTTGAAAATTTTGTAGACCGGATAAGGAATTCACGCCGCTTCCGGCATCATCAAAGCGAAATTTGTTATCCAATTATCAAGCCAAAGCGCCTTAGCTGGCGGCAATGAAAGCGAGCGTGGAATATCACcgacgt                                                                                                                                                                                      >  2:124700/1‑251 (MQ=255)
                                          tGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTCGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAAt                                                                                                                                                                                  >  1:371305/1‑251 (MQ=255)
                                                   ggtcggtGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCGGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAagc                                                                                                                                                                         >  1:5421/1‑251 (MQ=255)
                                                    gtcggtGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAgcgagcga                                                                                                                                                                                                      >  1:186569/1‑221 (MQ=255)
                                                      gtgtgccttgagtgaTGGGAAAGTCTGTCGTGGCCAATATCTAACGAAAAGGGTCTGGTGTGATATGATTAATAGGCCACGTCTGAATGTAATGCTTTGTATTTGCAAATTTTGTAGACCGGATAAGGAATTAACGCCGCATACGGCACCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCTGTAGCTGGCGGCAATGAAATCGAGCGAGGAATATCCCCGTCGTTAATGCGCCAatcaga                                                                                                                                                                      <  2:292971/243‑2 (MQ=255)
                                                        gtGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTCACGAAAATGGAGTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGACTTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCACCGTAGCTGACGGCAATGAAAGCGAGAGAGGAATATCAACGTCGTTAATACGCCAAGCAgcgc                                                                                                                                                                    >  2:323868/1‑251 (MQ=255)
                                                                                  cGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAACTTTTGTAGACCGGATAAGGAATTCACGCCGCATCAGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGCTTGCGTGTccc                                                                                                                                          >  1:167304/1‑251 (MQ=255)
                                                                                         taTAGCTAACGAAAATAGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCGGTTATCAAGCCAAAGCGCCGTAGATGGCGGAAATGAAAGCGAGCGGGGAATATCACAGTCGTTAATGCGCCACGCAGCACGCAGAACATATCGGATTGCGTGGCCCACTAGTc                                                                                                                                   >  2:141947/2‑251 (MQ=255)
                                                                                           tAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGCCCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAACGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTAGTTAATGCGCCAAGCAGCGCGCGGAACATATCGAATTGCGTGTCCCCCTGGTccc                                                                                                                                 >  2:310227/1‑249 (MQ=255)
                                                                                              cTAAAGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGGTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGAGGCAATGAAAGCGAGCGAGGAATATCACAGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGAATGCGTGTCCCACTGGACGcccc                                                                                                                              >  2:1677/1‑250 (MQ=255)
                                                                                                                 gtgtgctGTTTATCAGGCCACTGCTTAATGTAATTCTCTGGAGTTGTAAAATTTGTAGACCGGATAAGGTAATCACGCCGCATCCGGCCTCAACAAAGCGCAAGGTGTTATCCAGTTATCAAGCCAAAGCACCGGAGCTGGCGGAAATGAAAGCGAGCGAGGAATATCGCCGTCGTTAATTCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCTCACTTGTCGCTCTGGGTGCCGAGaaatgcat                                                                                                            <  2:152614/244‑6 (MQ=255)
                                                                                                                                                                 aaTTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATa                                                           <  1:72812/251‑1 (MQ=255)
                                                                                                                                                                      tGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATcg                                                      >  1:316680/1‑251 (MQ=255)
                                                                                                                                                                                              ccgcATCCGGCATCAACAAAGCGCCAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACgcagc                              >  1:136804/1‑251 (MQ=255)
                                                                                                                                                                                              ccgcATCCGGCATCAACAAAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATAGCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATACGATCCCTGACCCATTGCCAGCGCTGCCAACCACTCGAATAATTCATACATCCCGCTTATCGCCAGCGCCACgcagc                              >  1:316514/1‑251 (MQ=255)
                                                                                                                                                                                                     cGGCATCAAAAAAGCGCAAGTTGTTATACGGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGCCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACACGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAg                       <  2:5421/251‑1 (MQ=255)
                                                                                                                                                                                                            aacaaAGCGCAAGTTGTTATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGcc                >  1:327894/1‑251 (MQ=255)
                                                                                                                                                                                                               aaagcggcatgtgttATCCAGCTACCAAGCCCAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCCCCGTCGTTAATGCGCCAAGCAGCGCGCAGAACAGATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCTAGTAATGCATACATCGCGCTTATCGCCAGTGCCTCGAAGCAGACCAGGAACGccacc             <  2:136804/240‑1 (MQ=255)
                                                                                                                                                                                                                  gcgcAAGTTGTTATCCCGTTATCAAGCCAAACCGTCGTAGCTGGCGGCAATGCAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGCGTGCCGAGAAAATCATCCGCTCCCTGCCCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGTCCAGGAACGccaccaat          >  2:191632/1‑249 (MQ=255)
                                                                                                                                                                                                                           tttATCCAGTTATCAAGCCAAAGCGCCGTAGATGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGGTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTcc  <  2:165571/249‑1 (MQ=255)
                                                                                                                                                                                                                           gttATCCAGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTcc  <  1:310227/250‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
GCTGGATAAATATCGCAAAGTACCGCCCGAGCCTGTTCATGCTGATATTCCGGTCGGTGCAGTGAGTGATGTGAAAGCCTGACGTGGACAATAGCTAACGAAAATGGCCTGATGCGATATGTTTATCAGGCCACGTCTTAATGTAATGCTTTGGATTTGCAAATTTTGTAGACCGGATAAGGAATTCACGCCGCATCCGGCATCAACAAAGCGCAAGTTGTTATCCGGTTATCAAGCCAAAGCGCCGTAGCTGGCGGCAATGAAAGCGAGCGAGGAATATCACCGTCGTTAATGCGCCAAGCAGCGCGCAGAACATATCGGATTGCGTGTCCCACTGGTCGCCCTGGGTGCCGAGAAAATCATCCGCTCCCTGACCCATTGCCAGCGCTGCCCACCACTCGATTAATTCATACATCGCGCTTATCGCCAGCGCCACGCAGCAGACCAGGAACGCCACCATTTTACGTCC  >  NC_000913/4344709‑4345177

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