Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,082 T→C 26.5% I232I (ATT→ATC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,0820TC26.5% 52.4 / 21.6 34I232I (ATT→ATCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (16/9);  new base C (4/5);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTT  >  W3110S.gb/71023‑71149
                                                           |                                                                   
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:1816994/1‑51 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCAc                                           >  1:606764/1‑51 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCAc                                           >  1:1787661/1‑51 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                       >  1:1422799/1‑71 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                       >  1:1282635/1‑71 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                       >  1:2066462/1‑71 (MQ=255)
                                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                       >  1:2062802/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)
                                                           |                                                                   
CATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTT  >  W3110S.gb/71023‑71149

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: