Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,051,238 T→C intergenic (+286/+147) insB → / ← cspH IS1 transposase InsAB'/stress protein, member of the CspA‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,051,2380TC98.9% 323.4 / ‑4.3 87intergenic (+286/+147)insB/cspHIS1 transposase InsAB'/stress protein, member of the CspA‑family
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (72/14);  total (72/15)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AGTTTCAAGTGGTACTATAACTGGCTGTTAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTTG  >  W3110S.gb/1051210‑1051297
                            |                                                           
agattCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:393915/4‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2843847/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2454630/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:1648362/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:165712/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2841573/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:1764056/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:283061/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2797084/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:787696/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:600216/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2727037/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2573503/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2109681/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2553140/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:2408491/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:736877/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:816102/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:1038197/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:1074106/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:1347970/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:496827/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:473731/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAgaga                              >  1:3022678/1‑60 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGatatgg                         >  1:1129208/1‑63 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGatatag                         >  1:1087530/1‑63 (MQ=255)
aGTTTCAAGTGGTACTATAACTGGCTGTCAAATGGAATTTATAACTATTAGGTGCAgaga                              >  1:1697607/1‑60 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:2849554/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:2094051/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:762982/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:1102932/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:1952406/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:1936055/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:1898269/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:2981251/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:2917690/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:1817642/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:926854/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:1545798/65‑1 (MQ=255)
 gTTTCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATCCg                        <  1:396225/65‑1 (MQ=255)
  tttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCg                        <  1:588281/64‑1 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:760350/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:430409/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:293157/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:2929699/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:1168551/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:1317080/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:309727/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:173060/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:1042622/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:2784477/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:1893716/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:1944251/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAgg                                      >  1:612276/1‑49 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGc                       >  1:1859868/1‑64 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGc                       >  1:1820083/1‑64 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGc                       >  1:1234738/1‑64 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGc                       >  1:1964765/1‑64 (MQ=255)
   ttCAAGTGGTACTATAACTGGCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGc                       >  1:1137444/1‑64 (MQ=255)
    tCAAGTGGTACTATAACTGGCTGTTAAATGGCAtt                                                   <  1:401672/35‑1 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:679216/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:953096/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:930482/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:825377/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:793080/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:568593/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:623571/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2414534/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:1057972/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:1101401/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:1454177/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:1548371/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:1577068/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:1589995/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:1870800/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2075760/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2223282/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:330273/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2520073/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2576819/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2622440/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2736275/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2825106/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2825240/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2851112/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2915693/1‑65 (MQ=255)
                       gCTGTCAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTtg  >  1:2987663/1‑65 (MQ=255)
                            |                                                           
AGTTTCAAGTGGTACTATAACTGGCTGTTAAATGGCATTTATAACTATTAGGTGCAGAGATATTCGCTTAAATGGATTAGTTTCTTTG  >  W3110S.gb/1051210‑1051297

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: