Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 798,003 C→T 11.1% A343V (GCG→GTG)  modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb798,0030CT11.1% 129.8 / 6.9 45A343V (GCG→GTG) modCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base T (3/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGC  >  W3110S.gb/797946‑798069
                                                         |                                                                  
gTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCC‑‑TGGCTGTACGCGCAAATTAAAAg                                                        <  1:384067/68‑1 (MQ=255)
 tATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAg                                                        <  1:2586021/69‑1 (MQ=255)
   tCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtg                                                      >  1:900095/1‑69 (MQ=255)
   tCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtg                                                      >  1:1335399/1‑69 (MQ=255)
     aGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtgtg                                                    >  1:3617389/1‑69 (MQ=255)
     aGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtgtg                                                    >  1:3510726/1‑69 (MQ=255)
     aGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtgtg                                                    >  1:2076748/1‑69 (MQ=255)
     aGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtgtg                                                    >  1:1721135/1‑69 (MQ=255)
     aGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAgtgtg                                                    >  1:1668456/1‑69 (MQ=255)
        ccGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCg                                                 <  1:2421934/69‑1 (MQ=255)
            gggCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATaa                                             <  1:1654794/69‑1 (MQ=255)
                        aaCTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTaa                                                           >  1:1689694/1‑43 (MQ=255)
                        aaCTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTaa                                 <  1:1978688/69‑1 (MQ=255)
                        aaCTGGCGATCAAACATGGCCTGTGGCTGTACGCGCAAATTaa                                                           >  1:73819/1‑43 (MQ=255)
                           tGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                              <  1:1847560/69‑1 (MQ=255)
                           tGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                              <  1:2230017/69‑1 (MQ=255)
                           tGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                              <  1:1517789/69‑1 (MQ=255)
                            ggCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAg                                                        >  1:79091/1‑42 (MQ=255)
                                 tCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTg                         >  1:1890610/1‑68 (MQ=255)
                                  cAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                  >  1:1218179/1‑58 (MQ=255)
                                  cAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTa                                  <  1:348720/58‑1 (MQ=255)
                                  cAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGc                       >  1:3299150/1‑69 (MQ=255)
                                   aaaCCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGc                       >  1:2005733/1‑68 (MQ=255)
                                     aCCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCta                     <  1:1152171/68‑1 (MQ=255)
                                     aCCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCta                     <  1:445066/68‑1 (MQ=255)
                                     aCCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCta                     <  1:1579645/68‑1 (MQ=255)
                                     aCCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCta                     >  1:2088453/1‑68 (MQ=255)
                                      ccTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtata                   <  1:3595800/69‑1 (MQ=255)
                                      ccTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtata                   <  1:1651837/69‑1 (MQ=255)
                                      ccTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtata                   <  1:1816309/69‑1 (MQ=255)
                                      ccTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtata                   <  1:603911/69‑1 (MQ=255)
                                      ccTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtata                   <  1:1399477/69‑1 (MQ=255)
                                      ccTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCtata                   <  1:3192266/69‑1 (MQ=255)
                                            cTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATg               >  1:1044450/1‑67 (MQ=255)
                                            cTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATg               >  1:249788/1‑67 (MQ=255)
                                                   tGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc       <  1:2364169/68‑1 (MQ=255)
                                                   tGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGc       <  1:2678611/68‑1 (MQ=255)
                                                   tGTACGCGCAAATTAAAAGTGTGGCGATAACCGCCTGATTAAATc                              <  1:2061766/45‑1 (MQ=255)
                                                     tACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGa              <  1:2911212/59‑1 (MQ=255)
                                                      aCGTGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                              >  1:3360212/1‑42 (MQ=255)
                                                      aCGTGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                              <  1:3452807/42‑1 (MQ=255)
                                                      aCGTGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                              >  1:2523429/1‑42 (MQ=255)
                                                      aCGTGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                              <  1:1249599/42‑1 (MQ=255)
                                                      aCGTGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATc                              >  1:637500/1‑42 (MQ=255)
                                                      aCGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATg   >  1:3674414/1‑69 (MQ=255)
                                                       cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:3081889/1‑69 (MQ=255)
                                                       cgcgCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGc  >  1:564974/1‑69 (MQ=255)
                                                         |                                                                  
GTATCAGCCCGTGGGCCAGGGATGAACTGGCGATCAAACCTGGCCTGTGGCTGTACGCGCAAATTAAAAGTGTGTCGATAACCGCCTGATTAAATCAGGTGGCTATAAATGAACTGGGCAATGC  >  W3110S.gb/797946‑798069

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