Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,085 T→C 15.1% S233S (AGT→AGC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,0850TC15.1% 133.6 / 16.4 53S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (12/33);  new base C (3/5);  total (15/38)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

ACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGA  >  W3110S.gb/71022‑71153
                                                               |                                                                    
acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGa                                                               <  1:1143843/71‑1 (MQ=255)
acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGa                                                               <  1:976018/71‑1 (MQ=255)
acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGa                                                               <  1:850258/71‑1 (MQ=255)
acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGa                                                               <  1:764750/71‑1 (MQ=255)
acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGa                                                               <  1:764704/71‑1 (MQ=255)
acaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGa                                                               <  1:338763/71‑1 (MQ=255)
             cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                  <  1:1917346/71‑1 (MQ=255)
             cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                  <  1:1914764/71‑1 (MQ=255)
             cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                  <  1:1727232/71‑1 (MQ=255)
             cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                  <  1:157899/71‑1 (MQ=255)
             cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAg                                                  <  1:1428065/71‑1 (MQ=255)
              agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                                 <  1:743226/71‑1 (MQ=255)
              agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                                 <  1:1872395/71‑1 (MQ=255)
              agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                                 <  1:763274/71‑1 (MQ=255)
              agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGc                                                 <  1:877517/71‑1 (MQ=255)
              agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGc                                                 <  1:1386258/71‑1 (MQ=255)
              agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGc                                                 <  1:1051562/71‑1 (MQ=255)
                   ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAc                                            <  1:1758083/71‑1 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCAccc                                          >  1:760272/1‑71 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCAccc                                          >  1:471082/1‑71 (MQ=255)
                     aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCAccc                                          >  1:1602958/1‑71 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:997656/71‑1 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:553292/71‑1 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:36041/71‑1 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:354256/71‑1 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:870853/71‑1 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:264247/71‑1 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:1212126/71‑1 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:1969664/71‑1 (MQ=255)
                      gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:1503306/71‑1 (MQ=255)
                       ggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:838431/70‑1 (MQ=255)
                       ggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:1622252/70‑1 (MQ=255)
                         aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGAt                                      >  1:1583036/1‑71 (MQ=255)
                         aTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGAt                                      >  1:1392684/1‑71 (MQ=255)
                          ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAc                                            <  1:206755/64‑1 (MQ=255)
                          ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAc                                            <  1:955641/64‑1 (MQ=255)
                          ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAc                                            <  1:1309820/64‑1 (MQ=255)
                          ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAc                                            <  1:1338875/64‑1 (MQ=255)
                          ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAc                                            <  1:535245/64‑1 (MQ=255)
                          ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAc                                            <  1:1823029/64‑1 (MQ=255)
                          ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAc                                            <  1:193480/64‑1 (MQ=255)
                                             cgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:613887/48‑1 (MQ=255)
                                               aGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGc                <  1:48730/71‑1 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATg                                     >  1:184539/1‑36 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGcc                                   >  1:1183222/1‑38 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGcc                                   >  1:141168/1‑38 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGcc                                   >  1:50724/1‑38 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGcc                                   >  1:925927/1‑38 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATg            >  1:1466822/1‑61 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATg            >  1:1303884/1‑61 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa  >  1:1227908/1‑71 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa  >  1:66433/1‑71 (MQ=255)
                                                             aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa  >  1:1173631/1‑71 (MQ=255)
                                                               |                                                                    
ACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGA  >  W3110S.gb/71022‑71153

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: