Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,297,593 C→T 50.0% A62A (GCG→GCA cobT ← nicotinate‑nucleotide dimethylbenzimidazole‑P phophoribosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,297,5930CT50.0% ‑4.8 / 59.5 66A62A (GCG→GCAcobTnicotinate‑nucleotide dimethylbenzimidazole‑P phophoribosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (23/10);  new base T (23/10);  total (46/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

ACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGT  >  minE/1297528‑1297647
                                                                 |                                                      
aCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGc                                                      >  1:2843745/1‑68 (MQ=255)
aCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCt                                                     >  1:2955996/1‑69 (MQ=255)
aCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCt                                                     >  1:1993747/1‑69 (MQ=255)
aCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCt                                                     >  1:2515818/1‑69 (MQ=255)
aCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCt                                                     >  1:2596375/1‑69 (MQ=255)
aCTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCt                                                     >  1:815835/1‑69 (MQ=255)
  ttcttTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCttt                                                   <  1:229431/69‑1 (MQ=255)
  ttcttTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCttt                                                   <  1:445581/69‑1 (MQ=255)
                 aGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    >  1:2746550/1‑69 (MQ=255)
                 aGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    >  1:337751/1‑69 (MQ=255)
                 aGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    >  1:123809/1‑69 (MQ=255)
                 aGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    >  1:1922150/1‑69 (MQ=255)
                 aGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    >  1:2486604/1‑69 (MQ=255)
                 aGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCg                                     >  1:835296/1‑68 (MQ=255)
                 aGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACAAAACCAGTACCGCTTTTTTGCCCACATGCg                                     >  1:216003/1‑68 (MQ=255)
                       cccTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    <  1:2527165/63‑1 (MQ=255)
                       cccTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    <  1:911597/63‑1 (MQ=255)
                       cccTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    <  1:2058282/63‑1 (MQ=255)
                       cccTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    <  1:692473/63‑1 (MQ=255)
                       cccTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCgg                                    <  1:1699751/63‑1 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:1205275/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:44520/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:3133027/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:3253395/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:335379/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:847932/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:3377417/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:772762/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:3413984/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:752279/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:749145/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:568629/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:379709/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:498439/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:402387/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:1450062/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:1990422/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:2086851/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:2340991/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:2769156/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAccc                     >  1:2711702/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCAATTCAAccc                     >  1:405480/1‑69 (MQ=255)
                                ccAGACGCCGTGATCGGCACACATAACCAGTACTGCTTTTTTGCCCACATGAGGTATGCCATTCAAc                       >  1:3344968/1‑67 (MQ=255)
                                  aGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAAcccc                    >  1:1165675/1‑68 (MQ=255)
                                  aGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAAcccc                    >  1:1652305/1‑68 (MQ=255)
                                  aGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAAcccc                    >  1:2637117/1‑68 (MQ=255)
                                  aGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAAcccc                    >  1:3136575/1‑68 (MQ=255)
                                        cGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCAtt              <  1:3029737/68‑1 (MQ=255)
                                            aTCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGc                               >  1:2953610/1‑47 (MQ=255)
                                            aTCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGc                               >  1:3357721/1‑47 (MQ=255)
                                            aTCAGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGc                               >  1:2037064/1‑47 (MQ=255)
                                              cGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCa       >  1:1798852/1‑69 (MQ=255)
                                              cGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCa       >  1:3249698/1‑69 (MQ=255)
                                              cGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCa       >  1:2946529/1‑69 (MQ=255)
                                                   cacaTAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGt  <  1:710276/69‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCTGTATGCCATTCAACCCCgg                  <  1:467270/50‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAAcccctg                  <  1:1798707/50‑3 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCgg                  <  1:689244/50‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCgg                  <  1:2115229/50‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCgg                  <  1:698927/50‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCgg                  <  1:1746553/50‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCgg                  <  1:1681633/50‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCgg                  <  1:3311493/50‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCgg                  <  1:1277719/50‑1 (MQ=255)
                                                      aTAACCAGTACTGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCgg                  <  1:2942966/50‑1 (MQ=255)
                                                         aCCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCa                         <  1:730128/40‑1 (MQ=255)
                                                                ccGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCAt               >  1:2202410/1‑43 (MQ=255)
                                                                 |                                                      
ACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTACCGCTTTTTTGCCCACATGCGGTATGCCATTCAACCCCGGCATTCCTGCCAGTTGT  >  minE/1297528‑1297647

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