Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,679,525 G→T A1074D (GCT→GAT)  yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,679,5250GT96.3% 84.4 / ‑5.6 27A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (15/11);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCG  >  NC_000913/2679460‑2679588
                                                                 |                                                               
gggTTTAGTGAACGCGGTTAAAAATTAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATcc                                                               >  2:452551/1‑68 (MQ=255)
 ggTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCg                                                              >  2:227250/1‑68 (MQ=255)
  gTTTAGTGAACGCGGTTAAAAATAAAAGGCTTCACGGTGGTGAATAATGCCGCTTTCGGCGAATCCgg                                                             >  2:337666/1‑68 (MQ=255)
       gagAACGCGGTTAAAAATAAAAGGCAACACGGTGGTGAATACTGACGCTTTCGGCCAATCAGGATTaa                                                        >  1:290125/3‑68 (MQ=255)
       gTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTaa                                                        >  2:237699/1‑68 (MQ=255)
        tGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                       >  2:246725/1‑68 (MQ=255)
             gcgTTTAAAAATAATAGGCATCACGTTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                                  <  2:399230/68‑1 (MQ=255)
             gcgGTTAAAAATACAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                                  <  2:598359/68‑1 (MQ=255)
             gcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCg                                                  <  1:54104/68‑1 (MQ=255)
                            agtcATCACGGTTGTGAATATTGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCATCGCGAg                                   <  1:126791/65‑1 (MQ=255)
                             ggCATCACGGTGGTGAATAATGACGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGc                                  >  2:684332/1‑68 (MQ=255)
                              gCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCa                                 >  2:717784/1‑68 (MQ=255)
                                agcACGGTGGCGAATAATGCCGCTCTCGGCCAATCCGGATTACTCGCCGTCAATTCCGCGCGAGCAGa                               <  1:655715/66‑1 (MQ=255)
                                aTCACGGTTTTTAATAATTCAGCTTTCGGCCAAGTCGTATTAATCGCCGTCAATTCAGCGCGTGCAGa                               <  1:266307/68‑1 (MQ=255)
                                aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                               <  2:248708/68‑1 (MQ=255)
                                aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                               <  1:327818/68‑1 (MQ=255)
                                aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                               <  1:266309/68‑1 (MQ=255)
                                aTCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                               <  2:800472/68‑1 (MQ=255)
                                  cACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                             <  1:452551/68‑1 (MQ=255)
                                          ggAATAATGCCGCTTTCTGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgc                     <  2:228059/67‑1 (MQ=255)
                                             ataatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa                  >  2:466192/1‑68 (MQ=255)
                                               aatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCagc                >  2:652276/1‑66 (MQ=255)
                                               aatGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc                >  2:750796/1‑68 (MQ=255)
                                                  gacGCTTTCGGCCAATCAGGATTAATAGCAGTCAAATCAGAGAGAGCAGACTCTTAATGCGAAGCCTc             >  2:278314/3‑68 (MQ=255)
                                                  gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc             >  2:355170/1‑68 (MQ=255)
                                                  gCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc             >  2:170273/1‑68 (MQ=255)
                                                             ccAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCcg  >  1:352918/1‑68 (MQ=255)
                                                                 |                                                               
GGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCG  >  NC_000913/2679460‑2679588

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