Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 407,349 C→T 22.4% G44D (GGT→GAT)  cobC ← predicted alpha‑ribazole‑5'‑P phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE407,3490CT22.4% 84.2 / 16.2 62G44D (GGT→GAT) cobCpredicted alpha‑ribazole‑5'‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base C (33/15);  new base T (12/2);  total (45/17)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGC  >  minE/407287‑407413
                                                              |                                                                
ccAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGt                                                             >  1:3092368/1‑68 (MQ=255)
ccAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGt                                                             >  1:1430481/1‑68 (MQ=255)
ccAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGt                                                             >  1:1679526/1‑68 (MQ=255)
     cgcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTagcagc                                                       <  1:2381994/69‑1 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCGAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:2741610/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTagcagc                                                       <  1:2646383/68‑1 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTagcagc                                                       >  1:3287816/1‑68 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:2937462/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:2991959/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:2783692/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:259822/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:1961536/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:3375130/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:1133003/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:1318976/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:1422430/1‑69 (MQ=255)
      gcgcgGTATGCTGTGCCCGTGCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                                      >  1:3052592/1‑69 (MQ=255)
                cTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGAtt                                            >  1:1516843/1‑69 (MQ=255)
                cTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGAtt                                            >  1:2228324/1‑69 (MQ=255)
                          tCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                                  <  1:762777/69‑1 (MQ=255)
                          tCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                                  <  1:3172580/69‑1 (MQ=255)
                          tCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                                  <  1:3157224/69‑1 (MQ=255)
                           ccAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                                  >  1:1780337/1‑68 (MQ=255)
                           ccAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                                  >  1:393517/1‑68 (MQ=255)
                           ccAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                                  >  1:3331384/1‑68 (MQ=255)
                             aGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                <  1:357473/68‑1 (MQ=255)
                                      cATAAAGCCAGATCAAAGGAAACACCATGTagcagc                                                       >  1:443141/1‑36 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGAtt                                            >  1:2389980/1‑47 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                >  1:3146601/1‑59 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                >  1:696257/1‑59 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                >  1:1209717/1‑59 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                >  1:2753721/1‑59 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                >  1:1932976/1‑59 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc                      >  1:2506150/1‑69 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc                      >  1:2548901/1‑69 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc                      >  1:2555042/1‑69 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc                      >  1:2188714/1‑69 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc                      >  1:2121989/1‑69 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc                      >  1:2035066/1‑69 (MQ=255)
                                      cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgagc                      >  1:1305491/1‑69 (MQ=255)
                                          aaaCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgcgcg                   >  1:1071600/1‑68 (MQ=255)
                                          aaaCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgcgcg                   >  1:2883198/1‑68 (MQ=255)
                                          aaaCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgcgcg                   >  1:203799/1‑68 (MQ=255)
                                            aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGTTTTTGCGCTTGCTCAATACCGCGCGCGGt                 <  1:3209984/68‑1 (MQ=255)
                                            aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt                 <  1:3023967/68‑1 (MQ=255)
                                            aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt                 <  1:1313358/68‑1 (MQ=255)
                                            aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt                 <  1:508049/68‑1 (MQ=255)
                                            aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt                 <  1:1454125/68‑1 (MQ=255)
                                            aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt                 <  1:1968432/68‑1 (MQ=255)
                                                   cAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                <  1:1870717/46‑1 (MQ=255)
                                                   cAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                <  1:1556760/46‑1 (MQ=255)
                                                   cAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                <  1:3078195/46‑1 (MQ=255)
                                                   cAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                <  1:274669/46‑1 (MQ=255)
                                                   cAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                                <  1:703643/46‑1 (MQ=255)
                                                          aaCACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCTCGCGCGGTCAGGGGGGTGGgcgc  >  1:3280246/1‑69 (MQ=255)
                                                          aaCACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAggg            >  1:1628202/1‑59 (MQ=255)
                                                          aaCACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTggg     >  1:1033839/1‑66 (MQ=255)
                                                          aaCACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGgcgc  >  1:2052010/1‑69 (MQ=255)
                                                          aaCACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGgcgc  >  1:2268726/1‑69 (MQ=255)
                                                          aaCACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGgcgc  >  1:2593462/1‑69 (MQ=255)
                                                          aaCACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGgcgc  >  1:2681197/1‑69 (MQ=255)
                                                          aaCACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGgcgc  >  1:589815/1‑69 (MQ=255)
                                                              |                                                                
CCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGC  >  minE/407287‑407413

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