Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 797,006 G→A 21.7% G11S (GGC→AGC)  modC → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb797,0060GA21.7% 94.6 / 22.1 46G11S (GGC→AGC) modCmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (12/24);  new base A (3/7);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

GCTGGCCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCC  >  W3110S.gb/796937‑797074
                                                                     |                                                                    
gctggcCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTgg                                                                      <  1:207396/70‑1 (MQ=255)
    gcCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTgg                                                                      >  1:4124604/1‑66 (MQ=255)
    gcCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAAc                                                                 <  1:2351784/71‑1 (MQ=255)
    gcCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGAGCa                                                                   <  1:4718858/69‑1 (MQ=255)
      cAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAAc                                                                 <  1:385702/69‑1 (MQ=255)
        gAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCAtt                                                             <  1:788284/71‑1 (MQ=255)
           tCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCa                                                               <  1:4893782/66‑1 (MQ=255)
                 gTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCAtt                                                             >  1:3170889/1‑62 (MQ=255)
                 gTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCAtt                                                             >  1:475544/1‑62 (MQ=255)
                    aaCGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAAc                                                                 <  1:718584/55‑1 (MQ=255)
                    aaCGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTg                                                        >  1:2747920/1‑64 (MQ=255)
                    aaCGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTAtt                                                  <  1:4129574/70‑1 (MQ=255)
                    aaCGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTAtt                                                  <  1:4659871/70‑1 (MQ=255)
                    aaCGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGAGCAACCATTGCCTGACTAtt                                                  <  1:2440518/70‑1 (MQ=255)
                    aaCGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGAGCAACCATTGCCTGACTAtt                                                  <  1:201602/70‑1 (MQ=255)
                     aCGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTAtt                                                  <  1:180961/69‑1 (MQ=255)
                     aCGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCAATGCCTGACTAtt                                                  <  1:694953/69‑1 (MQ=255)
                      cgggcgggGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGc                                                                    >  1:2918099/1‑50 (MQ=255)
                      cgggcgggGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAAc                                                                 <  1:2764083/53‑1 (MQ=255)
                          cggggcgCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGaaa                                           <  1:2041869/71‑1 (MQ=255)
                               cgcTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCt                                       <  1:43203/70‑1 (MQ=255)
                                gcTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGAGCAACCATTGCCt                                                         <  1:1392379/51‑1 (MQ=255)
                                gcTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGAGCAACCATTGCCt                                                         <  1:5080905/51‑1 (MQ=255)
                                 cTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATg                                              <  1:4663372/61‑1 (MQ=255)
                                 cTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATg                                              <  1:4501948/61‑1 (MQ=255)
                                 cTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATg                                              <  1:2728311/61‑1 (MQ=255)
                                  tAATCATGCTGGAACTGAATTTTTCCCAGACGTTGAGCAACCATTGCCTGa                                                       >  1:157344/1‑51 (MQ=255)
                                  tAATCATGCTGGAACTGAATTTTTCCCAGACGTTGAGCAACCATTGCCTGa                                                       >  1:1960794/1‑51 (MQ=255)
                                   aaTCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCg                                  >  1:3640315/1‑71 (MQ=255)
                                      cATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAAc                                                                 <  1:4411283/37‑1 (MQ=255)
                                      cATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAAc                                                                 <  1:361854/37‑1 (MQ=255)
                                      cATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAAc                                                                 <  1:2791769/37‑1 (MQ=255)
                                      cATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGcc                                >  1:3408397/1‑70 (MQ=255)
                                       aTGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGaaa                                           >  1:1559301/1‑58 (MQ=255)
                                        tGCTGGAACTGAATTTTTCCCAGACGTTGGGCAAcc                                                                <  1:443200/36‑1 (MQ=255)
                                            ggAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGa                                                       >  1:687883/1‑41 (MQ=255)
                                             gAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCAt                        <  1:1342409/71‑1 (MQ=255)
                                               aCTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCa                      >  1:4971091/1‑71 (MQ=255)
                                                 tGAATTTTTCCCAGACGTTGAGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAAt                             <  1:2031860/62‑1 (MQ=255)
                                                  gAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCa                      >  1:492679/1‑68 (MQ=255)
                                                  gAATTTTTCCCAGACGTTGAGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATc                       <  1:1032835/67‑1 (MQ=255)
                                                    aTTTTTCCCAGACGTTGAGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCAc                     >  1:4951892/1‑67 (MQ=255)
                                                       tttCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCAc                     >  1:2082545/1‑64 (MQ=255)
                                                         tCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCtt            >  1:4044806/1‑71 (MQ=255)
                                                          cccAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCtt            <  1:1060978/70‑1 (MQ=255)
                                                           ccAGACGTTGGGCAACCATTGCCTGACTATTAATGAAAc                                          <  1:2801837/39‑1 (MQ=255)
                                                           ccAGACGTTGGGCAACCATTGCCTGACTATTAATGAAAc                                          <  1:4040033/39‑1 (MQ=255)
                                                               aCGTTGGGCAACCATTGCCTGACTATTAATGAAAc                                          <  1:5081135/35‑1 (MQ=255)
                                                                    gggCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTcc  <  1:3393008/70‑1 (MQ=255)
                                                                     |                                                                    
GCTGGCCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCC  >  W3110S.gb/796937‑797074

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: