Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 797,279 T→G 7.1% S102A (TCG→GCG)  modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb797,2790TG 7.1% 134.9 / 5.2 56S102A (TCG→GCG) modCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (30/22);  new base G (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GGCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCG  >  W3110S.gb/797210‑797346
                                                                     |                                                                   
ggCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGt                                                                     <  1:860457/70‑1 (MQ=255)
ggCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGt                                                                     <  1:710883/70‑1 (MQ=255)
ggCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGt                                                                     <  1:646610/70‑1 (MQ=255)
ggCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGt                                                                     <  1:2897598/70‑1 (MQ=255)
ggCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGt                                                                     <  1:2704573/70‑1 (MQ=255)
  cTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGa                                                                  <  1:2624627/71‑1 (MQ=255)
  cTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGa                                                                  <  1:1781459/71‑1 (MQ=255)
     tGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGaa                                                                 <  1:1156592/69‑1 (MQ=255)
     tGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGaa                                                                 <  1:2761203/69‑1 (MQ=255)
     tGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGaa                                                                 <  1:2309465/69‑1 (MQ=255)
       tttttCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAg                                                              >  1:195429/1‑70 (MQ=255)
       tttttCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAg                                                              <  1:2214937/70‑1 (MQ=255)
       tttttCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAg                                                              <  1:764949/70‑1 (MQ=255)
                aTGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGAt                                                    >  1:1817105/1‑71 (MQ=255)
                   cgcgGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAg                                                 >  1:3103079/1‑71 (MQ=255)
                   cgcgGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAg                                                 >  1:589557/1‑71 (MQ=255)
                           ttCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGAt                                                    <  1:2825009/60‑1 (MQ=255)
                           ttCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGAt                                           >  1:2569692/1‑69 (MQ=255)
                           ttCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGAt                                           >  1:3073810/1‑69 (MQ=255)
                            tCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGa                                            >  1:2984946/1‑67 (MQ=255)
                            tCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGa                                            >  1:1910773/1‑67 (MQ=255)
                             cccGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTa                                                            >  1:1441064/1‑50 (MQ=255)
                             cccGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGa                                            >  1:2315662/1‑66 (MQ=255)
                             cccGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGc                                       >  1:365755/1‑71 (MQ=255)
                                 cATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCtggt                                   >  1:788659/1‑71 (MQ=255)
                                 cATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCtggt                                   >  1:896457/1‑71 (MQ=255)
                                   ttACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTa                                                            <  1:883331/44‑1 (MQ=255)
                                     aCAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGgcg                               >  1:2783228/1‑71 (MQ=255)
                                      cAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCa                                                  >  1:2474102/1‑51 (MQ=255)
                                      cAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCa                                                  >  1:657224/1‑51 (MQ=255)
                                      cAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGgcgc                              >  1:527094/1‑71 (MQ=255)
                                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                             >  1:510117/1‑71 (MQ=255)
                                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                             >  1:3657013/1‑71 (MQ=255)
                                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                             <  1:2756624/71‑1 (MQ=255)
                                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                             >  1:2802145/1‑71 (MQ=255)
                                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                             >  1:1818285/1‑71 (MQ=255)
                                         aGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCa                                                  >  1:309441/1‑48 (MQ=255)
                                         aGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCttt                           <  1:3492804/71‑1 (MQ=255)
                                            gCGTGGCAATCTGCGCTTCGGCATGTCGAAAAGTATGGTCGAt                                                    <  1:3069982/43‑1 (MQ=255)
                                             cGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAgg                       <  1:2921882/71‑1 (MQ=255)
                                             cGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAgg                       <  1:2630556/71‑1 (MQ=255)
                                                ggCAATCTGCGCTACGGCATGTCGAAATGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCAt                    >  1:1116859/1‑71 (MQ=255)
                                                 gCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCAtt                   >  1:1776391/1‑71 (MQ=255)
                                                    aTCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCa                                                  >  1:1189607/1‑37 (MQ=255)
                                                    aTCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATTAGCTGGTGGCGCTTTTAGGCa                     <  1:2165240/66‑1 (MQ=255)
                                                    aTCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCa                     <  1:2193265/66‑1 (MQ=255)
                                                    aTCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCa                     <  1:1255820/66‑1 (MQ=255)
                                                    aTCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCa                     <  1:466442/66‑1 (MQ=255)
                                                       tGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCtgg                                    >  1:1717828/1‑48 (MQ=255)
                                                       tGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCg             >  1:2150709/1‑71 (MQ=255)
                                                       tGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCg             >  1:2140460/1‑71 (MQ=255)
                                                       tGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCg             >  1:3414509/1‑71 (MQ=255)
                                                       tGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCg             >  1:1803958/1‑71 (MQ=255)
                                                         cgcTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGtt           <  1:1553630/71‑1 (MQ=255)
                                                           cTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGc         <  1:3434153/71‑1 (MQ=255)
                                                               ggCATGGCGAAAAGCATGGTCAATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGa     <  1:388158/71‑1 (MQ=255)
                                                               ggCATGGCGAAAAGCATGGTCAATCAGTTCCATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGa     <  1:2750442/71‑1 (MQ=255)
                                                                gCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGAc    <  1:2378461/71‑1 (MQ=255)
                                                                  aTGGCGAAAAGCATGGTCAATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCg  >  1:872570/1‑71 (MQ=255)
                                                                  aTGGCGAAAAGCATGGTCAATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGATCg  >  1:3509786/1‑71 (MQ=255)
                                                                    gTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCttt                           <  1:1546484/44‑1 (MQ=255)
                                                                     tCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCAtt                   >  1:1661062/1‑51 (MQ=255)
                                                                     |                                                                   
GGCTATGTTTTTCAGGATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCG  >  W3110S.gb/797210‑797346

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