Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,079 C→T 26.3% R231R (CGC→CGT araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,0790CT26.3% 61.8 / 23.6 38R231R (CGC→CGTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (17/11);  new base T (4/6);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

TCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTT  >  W3110S.gb/71011‑71149
                                                                    |                                                                      
tCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTAt                                                                      <  1:2367943/71‑1 (MQ=255)
tCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCAt                                                                      <  1:1687706/71‑1 (MQ=255)
  gcgTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCAt                                                                      >  1:354906/1‑69 (MQ=255)
  gcgTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCAt                                                                      >  1:2539083/1‑69 (MQ=255)
          cacaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCAt                                                                      <  1:544574/61‑1 (MQ=255)
            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:2450344/71‑1 (MQ=255)
                                 gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                     <  1:1667996/71‑1 (MQ=255)
                                 gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                     <  1:2407224/71‑1 (MQ=255)
                                        gCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtttt                                                 >  1:1901910/1‑52 (MQ=255)
                                        gCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCtttt                                                 >  1:1819664/1‑52 (MQ=255)
                                        gCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCattt                                                 >  1:2393743/1‑52 (MQ=255)
                                               aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCa                       >  1:2341636/1‑71 (MQ=255)
                                               aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCAc                                           >  1:606764/1‑51 (MQ=255)
                                               aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCAc                                           >  1:1816994/1‑51 (MQ=255)
                                               aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCAc                                           >  1:1787661/1‑51 (MQ=255)
                                               aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                       >  1:1282635/1‑71 (MQ=255)
                                               aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                       >  1:1422799/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)
                                                                    |                                                                      
TCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTT  >  W3110S.gb/71011‑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: