Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 797,287 T→C 6.9% S104S (AGT→AGC modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb797,2870TC 6.9% 158.2 / 2.8 58S104S (AGT→AGCmodCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (37/17);  new base C (2/2);  total (39/19)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCA  >  W3110S.gb/797226‑797357
                                                             |                                                                      
aTGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGAt                                                               >  1:1817105/1‑71 (MQ=255)
   cgcgGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAg                                                            >  1:589557/1‑71 (MQ=255)
   cgcgGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAg                                                            >  1:3103079/1‑71 (MQ=255)
           ttCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGAt                                                               <  1:2825009/60‑1 (MQ=255)
           ttCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGAt                                                      >  1:3073810/1‑69 (MQ=255)
           ttCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGAt                                                      >  1:2569692/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:896457/1‑71 (MQ=255)
                 cATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCtggt                                              >  1:788659/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:2802145/1‑71 (MQ=255)
                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                                        <  1:2756624/71‑1 (MQ=255)
                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                                        >  1:1818285/1‑71 (MQ=255)
                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                                        >  1:3657013/1‑71 (MQ=255)
                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                                        >  1:510117/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:466442/66‑1 (MQ=255)
                                    aTCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCa                                <  1:1255820/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:1803958/1‑71 (MQ=255)
                                       tGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCg                        >  1:3414509/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)
                                                       gAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTa                                    >  1:269592/1‑43 (MQ=255)
                                                       gAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAAc                          >  1:1048220/1‑53 (MQ=255)
                                                       gAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTAc        <  1:236899/71‑1 (MQ=255)
                                                        aaaaGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTAcc       <  1:3371696/71‑1 (MQ=255)
                                                         aaaGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTAc        >  1:2122529/1‑69 (MQ=255)
                                                            gTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTa         >  1:180352/1‑65 (MQ=255)
                                                            gTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTa         >  1:3230676/1‑65 (MQ=255)
                                                            gTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTa         >  1:3104655/1‑65 (MQ=255)
                                                            gTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTa         >  1:1871120/1‑65 (MQ=255)
                                                             tATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCa  <  1:1671795/71‑1 (MQ=255)
                                                             |                                                                      
ATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCA  >  W3110S.gb/797226‑797357

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: