Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,085 T→C 25.7% 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,0850TC25.7% 55.2 / 20.8 35S233S (AGT→AGCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (17/9);  new base C (4/5);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGA  >  W3110S.gb/71023‑71153
                                                              |                                                                    
caTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGa                                                               <  1:1615258/70‑1 (MQ=255)
 aTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAg                                                                   <  1:2290798/65‑1 (MQ=255)
             agcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAg                                                  >  1:420663/1‑70 (MQ=255)
               cagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCAc                                               <  1:336776/71‑1 (MQ=255)
                agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACt                                              <  1:1016852/71‑1 (MQ=255)
                agTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACt                                              <  1:1846743/71‑1 (MQ=255)
                 gTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACTCATTAGTCAGGCGAAGCTGCTTTTGAGCACt                                              <  1:1791634/70‑1 (MQ=255)
                  ttAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGcaccac                                            <  1:1500006/71‑1 (MQ=255)
                    aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCa                                                                    >  1:1094183/1‑45 (MQ=255)
                    aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAg                                                                   >  1:884728/1‑46 (MQ=255)
                    aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGaa                                                              >  1:426704/1‑51 (MQ=255)
                    aGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAg                                                             >  1:707937/1‑52 (MQ=255)
                     gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                         <  1:1667996/71‑1 (MQ=255)
                     gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                         <  1:2450344/71‑1 (MQ=255)
                     gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                         <  1:2407224/71‑1 (MQ=255)
                            gCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtttt                                                     >  1:1819664/1‑52 (MQ=255)
                            gCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtttt                                                     >  1:1901910/1‑52 (MQ=255)
                            gCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCattt                                                     >  1:2393743/1‑52 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCa                           >  1:2341636/1‑71 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCAc                                               >  1:1787661/1‑51 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCAc                                               >  1:1816994/1‑51 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCAc                                               >  1:606764/1‑51 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                           >  1:2066462/1‑71 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                           >  1:2062802/1‑71 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                           >  1:1422799/1‑71 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                           >  1:1282635/1‑71 (MQ=255)
                                     gCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCAcc                         <  1:1632733/71‑1 (MQ=255)
                                        ggCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtc                      >  1:2072344/1‑71 (MQ=255)
                                            cgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCgtcggtc                  >  1:1440003/1‑71 (MQ=255)
                                               ggACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCa               <  1:1363749/71‑1 (MQ=255)
                                                     aCGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGcc                                   >  1:738010/1‑45 (MQ=255)
                                                     aCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg         <  1:591748/71‑1 (MQ=255)
                                                     aCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg         <  1:732405/71‑1 (MQ=255)
                                                        cATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtt      >  1:355955/1‑71 (MQ=255)
                                                            aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa  >  1:918248/1‑71 (MQ=255)
                                                              |                                                                    
CATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGA  >  W3110S.gb/71023‑71153

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: