Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 407,336 A→G 24.0% D48D (GAT→GAC cobC ← predicted alpha‑ribazole‑5'‑P phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE407,3360AG24.0% 68.8 / 17.8 58D48D (GAT→GACcobCpredicted alpha‑ribazole‑5'‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base G (12/2);  total (37/21)
Fisher's exact test for biased strand distribution p-value = 6.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GGTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGT  >  minE/407275‑407398
                                                             |                                                              
ggTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAgg                                                         <  1:1127874/69‑1 (MQ=255)
 gTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGa                                                        <  1:2611610/69‑1 (MQ=255)
 gTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGa                                                        <  1:582112/69‑1 (MQ=255)
 gTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGa                                                        <  1:1253770/69‑1 (MQ=255)
 gTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGa                                                        <  1:2673964/69‑1 (MQ=255)
 gTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGa                                                        <  1:2426907/69‑1 (MQ=255)
 gTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGa                                                        <  1:2901711/69‑1 (MQ=255)
    aCTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAac                                                     <  1:1958047/69‑1 (MQ=255)
    aCTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAac                                                     <  1:2219240/69‑1 (MQ=255)
            ccAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGt                                              >  1:1679526/1‑68 (MQ=255)
            ccAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGt                                              >  1:3092368/1‑68 (MQ=255)
            ccAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGt                                              >  1:1430481/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:1422430/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:259822/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:1133003/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:3375130/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:1318976/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:2783692/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:2937462/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:2991959/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:1961536/1‑69 (MQ=255)
                  gcgcgGTATGCTGTGCCCGTGCCAGTTCACTGCATAAAACCAGGTCAAAGGAAACATCATGTAGCAGCg                                       >  1:3052592/1‑69 (MQ=255)
                            cTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGAtt                             >  1:2228324/1‑69 (MQ=255)
                            cTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGAtt                             >  1:1516843/1‑69 (MQ=255)
                                      tCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                   <  1:762777/69‑1 (MQ=255)
                                      tCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                   <  1:3157224/69‑1 (MQ=255)
                                      tCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                   <  1:3172580/69‑1 (MQ=255)
                                       ccAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                   >  1:393517/1‑68 (MQ=255)
                                       ccAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                   >  1:3331384/1‑68 (MQ=255)
                                       ccAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGc                   >  1:1780337/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:696257/1‑59 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                 >  1:1209717/1‑59 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                 >  1:2753721/1‑59 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                 >  1:3146601/1‑59 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTc                 >  1:1932976/1‑59 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc       >  1:2548901/1‑69 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc       >  1:2035066/1‑69 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc       >  1:2121989/1‑69 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc       >  1:2188714/1‑69 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc       >  1:2506150/1‑69 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgc       >  1:2555042/1‑69 (MQ=255)
                                                  cATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgagc       >  1:1305491/1‑69 (MQ=255)
                                                      aaaCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgcgcg    >  1:203799/1‑68 (MQ=255)
                                                      aaaCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgcgcg    >  1:2883198/1‑68 (MQ=255)
                                                      aaaCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACcgcgcgcg    >  1:1071600/1‑68 (MQ=255)
                                                        aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGTTTTTGCGCTTGCTCAATACCGCGCGCGGt  <  1:3209984/68‑1 (MQ=255)
                                                        aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt  <  1:1968432/68‑1 (MQ=255)
                                                        aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt  <  1:3023967/68‑1 (MQ=255)
                                                        aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt  <  1:1454125/68‑1 (MQ=255)
                                                        aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt  <  1:1313358/68‑1 (MQ=255)
                                                        aCCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGt  <  1:508049/68‑1 (MQ=255)
                                                             |                                                              
GGTCACTGAGAACCAGTCGCGCGGTATGCTGTGCCCGTTCCAGTTCACTGCATAAAACCAGATCAAAGGAAACACCATGTAGCAGCGTATGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGT  >  minE/407275‑407398

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: