Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 503,377 C→T 40.2% D123D (GAC→GAT modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE503,3770CT40.2% 26.6 / 40.7 62D123D (GAC→GATmodCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (30/7);  new base T (20/5);  total (50/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 = 1.00e+00

TGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGACAGC  >  minE/503310‑503445
                                                                   |                                                                    
tGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGAcc                                                                     >  1:960578/1‑69 (MQ=255)
tGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGAcc                                                                     >  1:1730368/1‑69 (MQ=255)
tGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGAcc                                                                     >  1:1480831/1‑69 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:3089367/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:2890937/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:3357491/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:2688468/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:2659078/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:1013049/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:351790/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:398647/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:3179181/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:2110851/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:403689/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:556696/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:1894357/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:1803388/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:2033319/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:778866/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:1797770/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:1350225/1‑68 (MQ=255)
  tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCg                                                                    >  1:657882/1‑68 (MQ=255)
              gTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACGGTTTACCAGGCAGc                                                       >  1:1608448/1‑69 (MQ=255)
              gTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAg                                                        <  1:1407640/68‑1 (MQ=255)
              gTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAg                                                        <  1:1589267/68‑1 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGTATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:2235254/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTg                                                    >  1:1687571/1‑66 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:3295457/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:3176598/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:3286337/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:1154109/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:91326/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:3353645/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:887635/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:1022398/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:341047/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:808043/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:3412462/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:3416448/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:2422189/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:2228065/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:1576391/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:1947547/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:1672224/1‑68 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGt                                                   >  1:2740620/1‑67 (MQ=255)
                    cAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTAAACCGTTGCTTGATCGTTTACCAGGCAGCCTGTc                                                  >  1:2538125/1‑68 (MQ=255)
                                cTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGaaaa                                      <  1:260573/68‑1 (MQ=255)
                                           tttAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGc                          <  1:2735593/69‑1 (MQ=255)
                                             tAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGc                          <  1:2036277/67‑1 (MQ=255)
                                                     gAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGa                        >  1:3159544/1‑61 (MQ=255)
                                                     gAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCgg                >  1:424841/1‑69 (MQ=255)
                                                     gAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCgg                >  1:1021312/1‑69 (MQ=255)
                                                     gAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCgg                >  1:1261038/1‑69 (MQ=255)
                                                     gAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCgg                >  1:2029489/1‑69 (MQ=255)
                                                     gAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCgg                >  1:2178564/1‑69 (MQ=255)
                                                     gAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCgg                >  1:2538058/1‑69 (MQ=255)
                                                     gAACCGTTGCATAACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCgg                >  1:876026/1‑69 (MQ=255)
                                                          gTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCttt           >  1:1292704/1‑69 (MQ=255)
                                                                tGATCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGAc     <  1:2614489/69‑1 (MQ=255)
                                                                tGATCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGAc     <  1:76776/69‑1 (MQ=255)
                                                                tGATCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGAc     <  1:1346596/69‑1 (MQ=255)
                                                                tGATCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGAc     <  1:1178163/69‑1 (MQ=255)
                                                                tGATCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGAc     <  1:2752769/69‑1 (MQ=255)
                                                                tGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGAc     <  1:2693024/69‑1 (MQ=255)
                                                                tGACCGTTTACCAGGAAGCCTGTCCGGAGGCGATAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGAc     <  1:1452780/69‑1 (MQ=255)
                                                                   ccGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGTCGGGCTTTGCTGACAGc  <  1:2133763/68‑1 (MQ=255)
                                                                   |                                                                    
TGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGCCTGTCCGGAGGCGAAAAACAGCGCGTGGCGATTGGTCGGGCTTTGCTGACAGC  >  minE/503310‑503445

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: