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 / 91.9 51G118G (GGT→GGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (12/14);  new base G (17/8);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

ACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACG  >  W3110S.gb/70676‑70804
                                                                |                                                                
aCTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                              >  1:2915094/1‑69 (MQ=255)
 cTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                              <  1:919720/68‑1 (MQ=255)
 cTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                              <  1:812073/68‑1 (MQ=255)
 cTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                              <  1:3066668/68‑1 (MQ=255)
 cTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                              <  1:1591321/68‑1 (MQ=255)
                  attGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAgcg                                             <  1:683260/66‑1 (MQ=255)
                   cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAgcgc                                            <  1:187601/68‑1 (MQ=255)
                     ggCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAgcg                                             <  1:3506467/65‑1 (MQ=255)
                         tGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCa                                        <  1:3033033/66‑1 (MQ=255)
                          gAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgtt                                                                <  1:2055661/41‑1 (MQ=255)
                          gAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                    >  1:3253491/1‑69 (MQ=255)
                           aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttcttt                                                            >  1:3684609/1‑44 (MQ=255)
                           aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                    >  1:2718381/1‑68 (MQ=255)
                           aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                    >  1:2713933/1‑68 (MQ=255)
                           aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                    >  1:3043254/1‑68 (MQ=255)
                           aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATg                                                  >  1:2488833/1‑54 (MQ=255)
                           aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGc                                   >  1:2947305/1‑69 (MQ=255)
                           aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGc                                   >  1:2918782/1‑69 (MQ=255)
                           aaTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGc                                   >  1:3144348/1‑69 (MQ=255)
                            aTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCg                                                      <  1:1928220/49‑1 (MQ=255)
                            aTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCg                                                      >  1:2908024/1‑49 (MQ=255)
                               gCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAt                                 >  1:3407128/1‑67 (MQ=255)
                               gCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAt                                 >  1:3556043/1‑67 (MQ=255)
                                    aCTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCg                          >  1:795543/1‑69 (MQ=255)
                                    aCTGGCAGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCg                          >  1:789759/1‑69 (MQ=255)
                                     cTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCg                          <  1:430752/68‑1 (MQ=255)
                                          cGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTg                     <  1:3159218/68‑1 (MQ=255)
                                          cGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTg                     <  1:3163689/68‑1 (MQ=255)
                                          cGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTg                     <  1:1632694/68‑1 (MQ=255)
                                          cGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTg                     <  1:71011/68‑1 (MQ=255)
                                          cGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTg                     <  1:1628696/68‑1 (MQ=255)
                                          cGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTg                     <  1:738095/68‑1 (MQ=255)
                                           gTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCt                      <  1:753614/66‑1 (MQ=255)
                                              aaTATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTg                 >  1:88/1‑68 (MQ=255)
                                              aaTATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTg                 >  1:1232700/1‑68 (MQ=255)
                                              aaTATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTg                 >  1:3250035/1‑68 (MQ=255)
                                              aaTATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTg                 >  1:2454020/1‑68 (MQ=255)
                                                  tttGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaa            >  1:3181367/1‑69 (MQ=255)
                                                  tttGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCa             >  1:3408252/1‑68 (MQ=255)
                                                  tttGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCa             >  1:1352675/1‑68 (MQ=255)
                                                   ttGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaa            <  1:1496774/68‑1 (MQ=255)
                                                    tGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAAt          <  1:3344437/69‑1 (MQ=255)
                                                    tGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaaa           <  1:360902/68‑1 (MQ=255)
                                                    tGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaaa           <  1:2864446/68‑1 (MQ=255)
                                                     gCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc         >  1:1770794/1‑69 (MQ=255)
                                                      ccAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc         <  1:504289/68‑1 (MQ=255)
                                                        aaTACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAtt      <  1:1701486/69‑1 (MQ=255)
                                                        aaTACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAtt      <  1:3676244/69‑1 (MQ=255)
                                                          tACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAg                            >  1:635758/1‑45 (MQ=255)
                                                          tACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTaa    >  1:1886227/1‑69 (MQ=255)
                                                          tACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTa     >  1:2865369/1‑68 (MQ=255)
                                                           aCGGGTTTCTTTGGCCCGGTTGAAGCGCACCAGCCGCATTTc                              <  1:149461/42‑1 (MQ=38)
                                                           aCGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGa                         >  1:2639299/1‑47 (MQ=255)
                                                           aCGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAAc   <  1:1282296/69‑1 (MQ=255)
                                                           aCGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc         >  1:1111728/1‑63 (MQ=255)
                                                           aCGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAAc   >  1:2747374/1‑69 (MQ=255)
                                                           aCGGGGTTCTTTCGCCCGGATGAAGCGAACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAAc   >  1:180639/1‑69 (MQ=255)
                                                            cGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACg  >  1:1277361/1‑69 (MQ=255)
                                                                |                                                                
ACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACG  >  W3110S.gb/70676‑70804

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: