Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,591,445 G→A 22.4% T18I (ACA→ATA)  hsdR ← endonuclease R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,591,4450GA22.4% 140.6 / 35.9 67T18I (ACA→ATA) hsdRendonuclease R
Reads supporting (aligned to +/- strand):  ref base G (30/22);  new base A (6/9);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCT  >  W3110S.gb/4591380‑4591513
                                                                 |                                                                    
aGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATtattat                                                                  >  1:4436333/1‑70 (MQ=255)
aGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGttatt                                                                 <  1:1140126/71‑1 (MQ=255)
  gCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTaa                                                               <  1:3774426/71‑1 (MQ=255)
    gATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATc                                                             <  1:199662/71‑1 (MQ=255)
    gATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATc                                                             <  1:1045751/71‑1 (MQ=255)
    gATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATc                                                             <  1:3070374/71‑1 (MQ=255)
    gATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATc                                                             <  1:4309851/71‑1 (MQ=255)
        gCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGt                                                                     <  1:1834008/59‑1 (MQ=255)
        gCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGt                                                                     <  1:1367912/59‑1 (MQ=255)
        gCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCAt                                                          <  1:4579815/70‑1 (MQ=255)
          aTAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGttatta                                                                >  1:1793345/1‑62 (MQ=255)
          aTAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGttatta                                                                >  1:1799198/1‑62 (MQ=255)
          aTAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGttatta                                                                >  1:3254184/1‑62 (MQ=255)
          aTAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGc                                                       >  1:677409/1‑71 (MQ=255)
          aTAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATc                                                             <  1:1395245/65‑1 (MQ=255)
          aTAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATc                                                             <  1:1724670/65‑1 (MQ=255)
          aTAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGc                                                       >  1:4674460/1‑71 (MQ=255)
            aaGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCtg                                                     >  1:4718376/1‑71 (MQ=255)
             aGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCa                                                           >  1:560379/1‑64 (MQ=255)
             aGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCtgt                                                    >  1:2194657/1‑71 (MQ=255)
             aGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCtgt                                                    <  1:2041674/71‑1 (MQ=255)
                 aaGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGtta                                                                   >  1:3105168/1‑52 (MQ=255)
                 aaGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCgg                                                <  1:3609888/71‑1 (MQ=255)
                 aaGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCgg                                                <  1:1961155/71‑1 (MQ=255)
                 aaGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCgg                                                >  1:488054/1‑71 (MQ=255)
                 aaGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCgg                                                <  1:2082757/71‑1 (MQ=255)
                 aaGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCgg                                                <  1:131396/71‑1 (MQ=255)
                    tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATc                                                             <  1:2802253/55‑1 (MQ=255)
                    tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCggg                                               <  1:1877990/69‑1 (MQ=255)
                      gTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCattattatt                                                                 >  1:1923577/1‑49 (MQ=255)
                      gTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCattattatt                                                                 >  1:1638979/1‑49 (MQ=255)
                      gTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGc                                                  >  1:680724/1‑64 (MQ=255)
                      gTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCg                                                 >  1:3321124/1‑65 (MQ=255)
                        tGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTc                                         >  1:2604135/1‑71 (MQ=255)
                        tGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTc                                         >  1:2610001/1‑71 (MQ=255)
                        tGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGt                                          <  1:1448253/70‑1 (MQ=255)
                                 tCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCt                                                      >  1:4421255/1‑49 (MQ=255)
                                 tCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAtt                                 >  1:3364947/1‑70 (MQ=255)
                                 tCAGGAATTCAAAATTGGATTTATTCATCATTATTATTaa                                                               >  1:663990/1‑40 (MQ=255)
                                  cAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAttt                                <  1:138218/70‑1 (MQ=255)
                                  cAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAttt                                <  1:4594788/70‑1 (MQ=255)
                                    ggAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAttt                                <  1:4992300/68‑1 (MQ=255)
                                     gAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAg                             >  1:4679656/1‑70 (MQ=255)
                                      aaTTCAAAATTGGATTTATTCATCATTGTTATTAATCCa                                                           >  1:78336/1‑39 (MQ=255)
                                       aTTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGcc                          >  1:1458223/1‑71 (MQ=255)
                                         tCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCt                                            <  1:932603/51‑1 (MQ=255)
                                          cAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTTTGCGGGCCTGTCCAAATATTTAAGGccc                         >  1:2692683/1‑69 (MQ=255)
                                           aaaaTTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCAt                       <  1:2183807/70‑1 (MQ=255)
                                             aaTTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGGCCAAATAttt                                <  1:3947337/59‑1 (MQ=255)
                                              aTTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAca                   <  1:4325326/71‑1 (MQ=255)
                                              aTTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAca                   <  1:4714467/71‑1 (MQ=255)
                                                  gATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatct                >  1:1387896/1‑70 (MQ=255)
                                                   aTTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGcc                          >  1:4188427/1‑59 (MQ=255)
                                                    tttATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcat             >  1:3953173/1‑71 (MQ=255)
                                                     ttATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcatc            >  1:2444677/1‑71 (MQ=255)
                                                     ttATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcatc            >  1:588547/1‑71 (MQ=255)
                                                          catcatTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAg        <  1:4543545/70‑1 (MQ=255)
                                                           atcatTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCt      <  1:1247218/71‑1 (MQ=255)
                                                           atcatTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCt      >  1:2440901/1‑71 (MQ=255)
                                                              attattattAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTa         >  1:1200942/1‑65 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAttt                                >  1:2466747/1‑42 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAttt                                >  1:126222/1‑42 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATAttt                                >  1:2556817/1‑42 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAca                   >  1:1934682/1‑55 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCttt    <  1:4499199/70‑1 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCttt    <  1:1219287/70‑1 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTc   >  1:1940811/1‑71 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTc   >  1:4136754/1‑71 (MQ=255)
                                                              atTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTc   >  1:2751967/1‑71 (MQ=255)
                                                               ttGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAca                   >  1:4691225/1‑54 (MQ=255)
                                                               ttGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCt  >  1:4161161/1‑71 (MQ=255)
                                                                 gTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTaa                              >  1:197269/1‑41 (MQ=255)
                                                                 gTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTaa                              >  1:248524/1‑41 (MQ=255)
                                                                 gTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTaa                              >  1:2510248/1‑41 (MQ=255)
                                                                 |                                                                    
AGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTTTCT  >  W3110S.gb/4591380‑4591513

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: ‑