Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 797,294 G→A 6.3% D107N (GAT→AAT)  modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb797,2940GA 6.3% 192.5 / 2.0 63D107N (GAT→AAT) modCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (40/19);  new base A (2/2);  total (42/21)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

ATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGC  >  W3110S.gb/797226‑797359
                                                                    |                                                                 
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:2569692/1‑69 (MQ=255)
           ttCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGAt                                                        >  1:3073810/1‑69 (MQ=255)
            tCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGa                                                         >  1:1910773/1‑67 (MQ=255)
            tCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGa                                                         >  1:2984946/1‑67 (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)
                     aCAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGgcg                                            >  1:2783228/1‑71 (MQ=255)
                      cAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCa                                                               >  1:657224/1‑51 (MQ=255)
                      cAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCa                                                               >  1:2474102/1‑51 (MQ=255)
                      cAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGgcgc                                           >  1:527094/1‑71 (MQ=255)
                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                                          >  1:3657013/1‑71 (MQ=255)
                       aaaGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCt                                          >  1:510117/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)
                         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:3414509/1‑71 (MQ=255)
                                       tGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCg                          >  1:2150709/1‑71 (MQ=255)
                                       tGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCg                          >  1:2140460/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)
                                                       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:3230676/1‑65 (MQ=255)
                                                            gTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTa           >  1:1871120/1‑65 (MQ=255)
                                                            gTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTa           >  1:3104655/1‑65 (MQ=255)
                                                            gTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTa           >  1:180352/1‑65 (MQ=255)
                                                             tATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCa    <  1:1671795/71‑1 (MQ=255)
                                                              aTGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAg   >  1:1639142/1‑71 (MQ=255)
                                                               tGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGAc                 >  1:808518/1‑56 (MQ=255)
                                                               tGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGc  >  1:600331/1‑71 (MQ=255)
                                                               tGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGc  <  1:3144716/71‑1 (MQ=255)
                                                               tGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGc  <  1:3140306/71‑1 (MQ=255)
                                                               tGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGc  >  1:1772897/1‑71 (MQ=255)
                                                                    |                                                                 
ATGCGCGGCTGTTCCCGCATTACAAAGTGCGTGGCAATCTGCGCTACGGCATGTCGAAAAGTATGGTCGATCAGTTCGATAAGCTGGTGGCGCTTTTAGGCATTGAACCGTTGCTTGACCGTTTACCAGGCAGC  >  W3110S.gb/797226‑797359

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: