Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,237,592 G→A 13.2% V119V (GTG→GTA alx → predicted inner membrane protein, part of terminus

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,237,5920GA13.2% 134.2 / 10.2 53V119V (GTG→GTAalxpredicted inner membrane protein, part of terminus
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base A (5/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

GCTGATGTTGTTCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAGTTCGACT  >  W3110S.gb/3237526‑3237659
                                                                  |                                                                   
gCTGATGTTGTTCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGc                                                                    <  1:823486/68‑1 (MQ=255)
gCTGATGTTGTTCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGc                                                                    <  1:3073739/68‑1 (MQ=255)
 cTGATGTTGTTCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTACTCgg                                                                <  1:3565057/71‑1 (MQ=255)
 cTGATGTTGTTCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTACTCgg                                                                <  1:1969327/71‑1 (MQ=255)
    aTGTTGTTCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGc                                                             >  1:2555665/1‑71 (MQ=255)
          ttCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGt                                                       <  1:1545532/71‑1 (MQ=255)
          ttCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGt                                                       <  1:1713978/71‑1 (MQ=255)
          ttCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGt                                                       <  1:324470/71‑1 (MQ=255)
               cTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGc                                                  >  1:1856950/1‑71 (MQ=255)
                   ttCTCTGTTCCGGCGTCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTAc                                              >  1:2064247/1‑71 (MQ=255)
                    tctcTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTAc                                              >  1:1670592/1‑70 (MQ=255)
                     ctctGTTCCGGCGGCATT‑‑AACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACcat                                           <  1:2217934/70‑1 (MQ=255)
                          ttCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACcatca                                         <  1:2322882/69‑1 (MQ=255)
                          ttCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATg                                       <  1:2509808/71‑1 (MQ=255)
                           tCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTACTCGGGGCGATTGTTCTGCGTACCATCATGa                                      >  1:3578715/1‑71 (MQ=255)
                           tCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTACTCGGGGCGATTGTTCTGCGTACCATCATGa                                      >  1:3009971/1‑71 (MQ=255)
                            ccGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGa                                      >  1:250997/1‑70 (MQ=255)
                              ggcggcATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCt                                   >  1:296598/1‑71 (MQ=255)
                              ggcggcATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCt                                   >  1:1300727/1‑71 (MQ=255)
                               gcggcATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCtt                                  <  1:946329/71‑1 (MQ=255)
                               gcggcATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCtt                                  <  1:2179781/71‑1 (MQ=255)
                               gcggcATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCtt                                  <  1:1848914/71‑1 (MQ=255)
                                 ggcATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCa                                <  1:3378457/71‑1 (MQ=255)
                                   cATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCa                                >  1:3128573/1‑69 (MQ=255)
                                    aTTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATg                                       >  1:2007268/1‑61 (MQ=255)
                                    aTTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACt                              >  1:989954/1‑70 (MQ=255)
                                     ttACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATg                                       <  1:1310419/60‑1 (MQ=255)
                                     ttACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATg                                       <  1:2026379/60‑1 (MQ=255)
                                      tACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATg                                       <  1:2025470/59‑1 (MQ=255)
                                       aCAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTgg                            >  1:980892/1‑69 (MQ=255)
                                       aCAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGAGTACCATCATGATCTTCACTgg                            >  1:2723883/1‑69 (MQ=255)
                                          aCGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCt                       >  1:2074119/1‑71 (MQ=255)
                                          aCGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCt                       >  1:2906220/1‑71 (MQ=255)
                                          aCGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCt                       >  1:1803602/1‑71 (MQ=255)
                                           cgccgcGTGCTGGTGTATGGCGTGCTCGGGGCGATTg                                                        <  1:1423870/37‑1 (MQ=255)
                                              cgcgTGCTGGTGTATGGCGTACTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGTAGCTGGCt                   >  1:3422109/1‑71 (MQ=255)
                                              cgcgTGCTGGTGTATGGCGTACTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGTAGCTGGCt                   >  1:3791049/1‑71 (MQ=255)
                                                cgTGCTGGTGTATGGCGTACTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGTAGCTGGCtga                 >  1:3756291/1‑69 (MQ=255)
                                                 gTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAAt                <  1:2920611/71‑1 (MQ=255)
                                                    cTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTc                                 >  1:3813256/1‑51 (MQ=255)
                                                        tgtATGGCGTGCTCTGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAg         <  1:1871110/71‑1 (MQ=255)
                                                        tgtATGGCGTGCTCTGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAg         <  1:1875456/71‑1 (MQ=255)
                                                          tATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCtt                                  <  1:3695136/44‑1 (MQ=255)
                                                          tATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCtt                                  <  1:2001072/44‑1 (MQ=255)
                                                            tGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAGTTc      >  1:729968/1‑70 (MQ=255)
                                                            tGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAGTTc      >  1:810997/1‑70 (MQ=255)
                                                            tGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAGTTc      >  1:2726153/1‑70 (MQ=255)
                                                            tGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAGTTc      >  1:1017423/1‑70 (MQ=255)
                                                            tGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATCATCTTCACTGGCAGCt                       >  1:795218/1‑53 (MQ=255)
                                                              gCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGc                        >  1:1016957/1‑50 (MQ=255)
                                                              gCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGc                        >  1:2159043/1‑50 (MQ=255)
                                                              gCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGc                        >  1:3835548/1‑50 (MQ=255)
                                                               cGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAGTTCGACt  >  1:1628010/1‑71 (MQ=255)
                                                                  |                                                                   
GCTGATGTTGTTCAGCTATTTCTCTGTTCCGGCGGCATTACAACGCCGCGTGCTGGTGTATGGCGTGCTCGGGGCGATTGTTCTGCGTACCATCATGATCTTCACTGGCAGCTGGCTAATTTCGCAGTTCGACT  >  W3110S.gb/3237526‑3237659

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: