Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,740 T→G 49.1% G118G (GGT→GGG araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb70,7400TG49.1% 28.3 / 89.7 49G118G (GGT→GGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (10/15);  new base G (12/12);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-01

TTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCG  >  W3110S.gb/70673‑70807
                                                                   |                                                                   
tttACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttc                                                                  <  1:4190636/71‑1 (MQ=255)
tttACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttc                                                                  <  1:3269402/71‑1 (MQ=255)
tttACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttc                                                                  <  1:2220868/71‑1 (MQ=255)
 ttACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgtt                                                                   >  1:3610591/1‑69 (MQ=255)
  tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttctt                                                                >  1:1073047/1‑71 (MQ=255)
   aCTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttctt                                                                >  1:3725550/1‑70 (MQ=255)
      ttCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGc                                                            >  1:1433774/1‑71 (MQ=255)
       tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGcc                                                           >  1:674080/1‑71 (MQ=255)
            cgcgcgCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGAt                                                      >  1:1246251/1‑71 (MQ=255)
            cgcgcgCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTAGCCAATACGGGttt                                                                   >  1:3549903/1‑58 (MQ=255)
                cgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAg                                                  <  1:3343383/71‑1 (MQ=255)
                cgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAg                                                  <  1:2241344/71‑1 (MQ=255)
                      cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCAc                                             >  1:4436835/1‑70 (MQ=255)
                       tGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCAc                                             >  1:3891406/1‑69 (MQ=255)
                        ggCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttcttt                                                               >  1:4183978/1‑50 (MQ=255)
                        ggCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttc                                                                  <  1:470404/47‑1 (MQ=255)
                        ggCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttc                                                                  <  1:2278917/47‑1 (MQ=255)
                        ggCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttc                                                                  >  1:3529399/1‑47 (MQ=255)
                           aTGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                       >  1:1711491/1‑71 (MQ=255)
                           aTGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                       >  1:856600/1‑71 (MQ=255)
                            tGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCg                                                         <  1:330888/52‑1 (MQ=255)
                              aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                 >  1:462413/1‑42 (MQ=255)
                              aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                 >  1:1841363/1‑42 (MQ=255)
                              aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttct                                                                 >  1:128787/1‑42 (MQ=255)
                              aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCGGCCGCAt                                    >  1:742126/1‑71 (MQ=255)
                              aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAt                                    >  1:1559772/1‑71 (MQ=255)
                               aTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAtt                                   <  1:2325002/71‑1 (MQ=255)
                               aTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAtt                                   >  1:1545854/1‑71 (MQ=255)
                                tGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                  <  1:3321206/71‑1 (MQ=255)
                                tGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                  <  1:1896918/71‑1 (MQ=255)
                                 ggCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                  >  1:4399850/1‑70 (MQ=255)
                                  gCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTc                                 <  1:4184192/70‑1 (MQ=255)
                                  gCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                  >  1:3337391/1‑69 (MQ=255)
                                     tAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGcc                                                           <  1:1392240/41‑1 (MQ=255)
                                      aaCTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGc                              <  1:3725043/69‑1 (MQ=255)
                                      aaCTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCg                             <  1:2384246/70‑1 (MQ=255)
                                        cTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTc                                 <  1:376475/64‑1 (MQ=255)
                                         tGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAgcg                                                <  1:2537613/48‑1 (MQ=255)
                                           gCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGt                       <  1:1638405/71‑1 (MQ=255)
                                           gCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGt                       <  1:232555/71‑1 (MQ=255)
                                            ccttCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                      <  1:2380163/68‑1 (MQ=255)
                                              gTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGc                              <  1:5076272/61‑1 (MQ=255)
                                               tCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAg                                                  <  1:1585065/40‑1 (MQ=255)
                                                 aaTATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTggg                  >  1:4435953/1‑70 (MQ=255)
                                                 aaTATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTgg                   >  1:129257/1‑69 (MQ=255)
                                                  atatTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                       <  1:2962661/48‑1 (MQ=255)
                                                  atatTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                  >  1:3997279/1‑53 (MQ=255)
                                                  atatTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCa                <  1:2113753/71‑1 (MQ=255)
                                                    atTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGAcc                          >  1:333912/1‑59 (MQ=255)
                                                    atTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaaa              >  1:3074124/1‑71 (MQ=255)
                                                     tttGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaaa              >  1:1920058/1‑70 (MQ=255)
                                                     tttGCCAATACGGGGATCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAAt             <  1:1114730/71‑1 (MQ=255)
                                                       tGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCa           >  1:1749208/1‑71 (MQ=255)
                                                          cAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAtt         <  1:2418656/70‑1 (MQ=255)
                                                           aaTACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCa                <  1:3693244/62‑1 (MQ=255)
                                                            aTACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAt                                    <  1:2015561/41‑1 (MQ=255)
                                                            aTACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAt                                    <  1:713605/41‑1 (MQ=255)
                                                             tACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc            <  1:4372172/64‑1 (MQ=255)
                                                             tACGGGGCTCTTTCGCCCGGATGAAGCGCACCAGcc                                        >  1:1006114/1‑36 (MQ=37)
                                                               cGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTgg                   >  1:4466866/1‑55 (MQ=255)
                                                               cGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGc    <  1:3778206/70‑1 (MQ=255)
                                                                gggTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAAt             <  1:3068752/60‑1 (MQ=255)
                                                                gggTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCg  >  1:3536171/1‑71 (MQ=255)
                                                                   tttctttcGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAAc      <  1:2822480/64‑1 (MQ=255)
                                                                   |                                                                   
TTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCG  >  W3110S.gb/70673‑70807

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