Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,085,352 C→G 15.4% A58A (GCG→GCC)  yeeE ← UPF0394 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,085,3520CG15.4% 66.2 / 5.3 26A58A (GCG→GCC) yeeEUPF0394 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base G (2/2);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACC  >  NC_000913/2085205‑2085492
                                                                                                                                                   |                                                                                                                                            
gATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACg                                                                                                                                             <  1:562530/149‑1 (MQ=255)
  tCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAAGGCa                                                                                                                                           >  1:532525/1‑149 (MQ=255)
  tCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAAGGCa                                                                                                                                           >  2:410208/1‑149 (MQ=255)
    aaTCCCTCACCCGCGCGATACCAGGTACCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAGGTCAGTAGGCCCGCCTGAATTAACGaaaa                                                                                                                                         >  1:339971/1‑149 (MQ=255)
        ccTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCGTCGTAGGTCCGTAGGCGCGCCTGAATTAACGCAACAAcc                                                                                                                                     >  1:266885/1‑149 (MQ=255)
               cgcgcgATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGc                                                                                                                              >  2:485337/1‑149 (MQ=255)
                     aTACCAGGGCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTAGAACAGGACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTa                                                                                                                        <  2:309529/149‑1 (MQ=255)
                              cccGGTGACACATCCACCGGCCAGAACAATTCCCAGCCCGAAGATATACCCACGTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAACCCCCCACGCTTTGTACCGAAATTg                                                                                                               <  2:511367/149‑1 (MQ=255)
                                  gTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAAt                                                                                                           >  2:557841/1‑149 (MQ=255)
                                  gTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATGGCAAt                                                                                                           >  2:222048/1‑149 (MQ=255)
                                       acaTCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCa                                                                                                      <  1:485337/149‑1 (MQ=255)
                                                                    cGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCCATTATTTTTCACg                                                                         >  2:217135/1‑149 (MQ=255)
                                                                       aGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAAGGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGata                                                                      <  1:410208/149‑1 (MQ=255)
                                                                       aGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAAGGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGata                                                                      <  2:532525/149‑1 (MQ=255)
                                                                                 aCCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATcgc                                                            >  1:546549/1‑149 (MQ=255)
                                                                                     ataACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGATATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaaa                                                        <  2:511732/149‑1 (MQ=255)
                                                                                           gTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAaccacc                                                  <  1:582111/149‑1 (MQ=255)
                                                                                               cAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                              >  1:569591/1‑149 (MQ=255)
                                                                                                      cGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCaaa                                       <  1:625610/149‑1 (MQ=255)
                                                                                                              ccccGTCTCCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCa                               <  1:70735/149‑1 (MQ=255)
                                                                                                                 cGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCAccc                            >  2:602896/1‑147 (MQ=255)
                                                                                                                    cTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGGAAAACATGCGATTATTTTTCACGATCTACATATCGCGACAAAAACCTGTCAGCCACAAACGCCCACgacgc                         >  2:73681/1‑145 (MQ=255)
                                                                                                                     ttCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCa                        <  1:350362/149‑1 (MQ=255)
                                                                                                                                       cGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCagag      >  2:223905/1‑149 (MQ=255)
                                                                                                                                         ccTGAATTAACGCAAAACCCCCCACGCTTTGTCCCGAAATTGTAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGca    <  2:208103/149‑1 (MQ=255)
                                                                                                                                           tGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGCAAACCACCTGTCAGGCAAAAACGCCCACGCTGAATGACAAATCCAAGCAGAGcaca  >  2:583926/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACATCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACC  >  NC_000913/2085205‑2085492

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