Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,591,162 C→A 75.0% V112V (GTG→GTT hsdR ← endonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,591,1620CA75.0% 29.6 / 24.4 32V112V (GTG→GTThsdRendonuclease R
Reads supporting (aligned to +/- strand):  ref base C (6/2);  new base A (14/10);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGG  >  W3110S.gb/4591099‑4591226
                                                               |                                                                
aGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCt                                                           >  1:1103907/1‑71 (MQ=255)
 gCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTg                                                          <  1:1539937/71‑1 (MQ=255)
 gCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTg                                                          <  1:2398713/71‑1 (MQ=255)
   cAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGt                                                        <  1:2454529/71‑1 (MQ=255)
       cGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                    <  1:2610648/71‑1 (MQ=255)
        ggAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                    >  1:948328/1‑70 (MQ=255)
        ggAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                    >  1:2361493/1‑70 (MQ=255)
         gAACCCGAGTCTCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCaa                                                  <  1:951857/71‑1 (MQ=255)
                 gTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGt                                          <  1:2017410/71‑1 (MQ=255)
                    gCAGGCACATCTGGGCATCGTTGAGATCGTTATGATTTTCGTGAACCGCCTGGTTACCAATGCGGCGTAAt                                       >  1:2479294/1‑71 (MQ=255)
                    gCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTaa                                        >  1:1744808/1‑70 (MQ=255)
                         cacaTCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTtgt                                   >  1:1943741/1‑70 (MQ=255)
                          acaTCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCaa                                                  >  1:201384/1‑54 (MQ=255)
                            aTCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGaaa                               <  1:33268/71‑1 (MQ=255)
                             tCTGGGCATCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGa                                 <  1:1149815/68‑1 (MQ=255)
                               tGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCaa                                                  >  1:2551243/1‑49 (MQ=255)
                                     tCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                    <  1:1762684/41‑1 (MQ=255)
                                     tCGTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTAcc                                                    <  1:2457444/41‑1 (MQ=255)
                                       gTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGa                    >  1:122314/1‑71 (MQ=255)
                                       gTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATa                             >  1:888893/1‑62 (MQ=255)
                                       gTTGAGATCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGa                    >  1:1970420/1‑71 (MQ=255)
                                          gagaTCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGAtgt                 >  1:2163616/1‑71 (MQ=255)
                                              tCGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGtc             <  1:904704/71‑1 (MQ=255)
                                               cGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTaa                                        >  1:2751523/1‑43 (MQ=255)
                                               cGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATa                             >  1:1724685/1‑54 (MQ=255)
                                               cGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGtca            >  1:2503049/1‑71 (MQ=255)
                                               cGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGtca            >  1:2259708/1‑71 (MQ=255)
                                               cGTTATGATATTCGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGtca            >  1:2691632/1‑71 (MQ=255)
                                                      atatTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATcaaca      <  1:266145/70‑1 (MQ=255)
                                                        atTCGTGAACCGCCTGGTTACCAATGCGGCGTAAt                                       >  1:1787631/1‑35 (MQ=255)
                                                           cGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATa                             >  1:2646505/1‑42 (MQ=255)
                                                           cGTGAACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAgg  >  1:2683221/1‑69 (MQ=255)
                                                               |                                                                
AGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGG  >  W3110S.gb/4591099‑4591226

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: