Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,855,731 G→T 12.9% intergenic (‑269/‑18) ygbA ← / → mutS hypothetical protein/methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,855,7310GT12.9% 167.0 / 8.1 54intergenic (‑269/‑18)ygbA/mutShypothetical protein/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (17/30);  new base T (4/3);  total (21/33)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

TATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTA  >  W3110S.gb/2855668‑2855798
                                                               |                                                                   
tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGTAACCgg                                                              <  1:9062/71‑1 (MQ=255)
tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCgg                                                              <  1:2189242/71‑1 (MQ=255)
  tCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa                                                             <  1:1337240/70‑1 (MQ=255)
  tCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa                                                             <  1:3889494/70‑1 (MQ=255)
  tCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa                                                             <  1:1935422/70‑1 (MQ=255)
  tCATGTGCGCCTTATGTGATTACAACCAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa                                                             <  1:4637267/70‑1 (MQ=255)
    aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGAc                                                            >  1:1746979/1‑69 (MQ=255)
    aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACa                                                           >  1:5106685/1‑70 (MQ=255)
    aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACa                                                           >  1:3919250/1‑70 (MQ=255)
    aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACAt                                                          <  1:1072738/71‑1 (MQ=255)
          gcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAAcccc                                                    <  1:2470756/71‑1 (MQ=255)
            cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt                                                  <  1:95789/71‑1 (MQ=255)
             ttATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGa                                                                   <  1:5076217/53‑1 (MQ=255)
               aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGaa                                                                  >  1:2851972/1‑52 (MQ=255)
               aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATg                                                 <  1:3557281/69‑1 (MQ=255)
               aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAg                                               >  1:1965695/1‑71 (MQ=255)
                       aCAACTAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                          >  1:4070476/1‑68 (MQ=255)
                       aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                          >  1:1512268/1‑68 (MQ=255)
                       aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                          >  1:4079594/1‑68 (MQ=255)
                        cAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAg                                       <  1:4323522/70‑1 (MQ=255)
                        cAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGa                                      <  1:1109101/71‑1 (MQ=255)
                         aaCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGAc                                                            <  1:3140541/48‑1 (MQ=255)
                         aaCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaa                                     <  1:1948917/71‑1 (MQ=255)
                         aaCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaa                                     <  1:257582/71‑1 (MQ=255)
                             aaaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCATAGAAAAttt                                >  1:584102/1‑71 (MQ=255)
                               aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                               <  1:3945470/71‑1 (MQ=255)
                                aTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATaa                                                        >  1:5094849/1‑45 (MQ=255)
                                    aaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgc                          <  1:4336931/71‑1 (MQ=255)
                                      aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGa                             >  1:1392546/1‑66 (MQ=255)
                                       ccATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgccca                       <  1:4721304/71‑1 (MQ=255)
                                          tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                          <  1:382216/49‑1 (MQ=255)
                                          tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                          <  1:2908267/49‑1 (MQ=255)
                                          tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcc                         <  1:1587774/66‑1 (MQ=255)
                                            acacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccata                     <  1:5018441/68‑1 (MQ=255)
                                             cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAAt                                  >  1:888996/1‑54 (MQ=255)
                                             cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcc                 <  1:1635658/71‑1 (MQ=255)
                                             cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccat                      >  1:3563832/1‑66 (MQ=255)
                                             cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgccca                       >  1:1798739/1‑65 (MQ=255)
                                                 taattTAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggca                    <  1:85384‑M1/62‑34 (MQ=255)
                                                   aTTTAATATCAGGGAACCGGACATAACCCCATGAGTGc                                            <  1:3489432/38‑1 (MQ=255)
                                                    tttAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagt          <  1:2030762‑M1/71‑44 (MQ=255)
                                                    tttAATATCAGTAAACCGGACATAACCCatgaagaaaaa                                          >  1:3642286‑M1/1‑28 (MQ=255)
                                                    tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                               <  1:3207574/50‑1 (MQ=255)
                                                    tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                               <  1:3213705/50‑1 (MQ=255)
                                                    tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcc                 <  1:4001933/64‑1 (MQ=255)
                                                    tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgat           <  1:1546536/70‑1 (MQ=255)
                                                     ttAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagta         >  1:1119572‑M1/1‑27 (MQ=255)
                                                      tAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagtat        >  1:4048855‑M1/1‑26 (MQ=255)
                                                      tAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagtat        >  1:2371928‑M1/1‑26 (MQ=255)
                                                      tAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatga            >  1:1636736/1‑67 (MQ=255)
                                                       aaTATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAAt                                  >  1:4122962/1‑44 (MQ=255)
                                                       aaTATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATgca        <  1:4417024/70‑1 (MQ=255)
                                                        atatCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagtatc       <  1:1480321‑M1/70‑47 (MQ=255)
                                                            cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTa  >  1:3843244/1‑71 (MQ=255)
                                                               |                                                                   
TATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTA  >  W3110S.gb/2855668‑2855798

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: