Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,109 T→C 59.0% T241T (ACT→ACC araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,1090TC59.0% 20.2 / 44.2 39T241T (ACT→ACCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (6/10);  new base C (9/14);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTA  >  W3110S.gb/71052‑71178
                                                         |                                                                     
cGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTa                                                          >  1:956233/1‑71 (MQ=255)
  tCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATc                                                        <  1:1094685/71‑1 (MQ=255)
    ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                                  <  1:763568/59‑1 (MQ=255)
    ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                                  <  1:761552/59‑1 (MQ=255)
    ttAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                                  <  1:141067/59‑1 (MQ=255)
     tAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGc                                                      >  1:1198161/1‑70 (MQ=255)
      aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGTTGCTTTTGAGCACCACCCGGATGCCTATCGc                                                      <  1:1593218/69‑1 (MQ=255)
      aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCg                                                       <  1:156642/68‑1 (MQ=255)
      aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGcc                                                     <  1:830762/70‑1 (MQ=255)
      aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                                                    >  1:694422/1‑71 (MQ=255)
       aGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCAc                                                   <  1:1459624/71‑1 (MQ=255)
          tGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgt                                                >  1:656082/1‑71 (MQ=255)
            gcgcGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcg                                              <  1:94546/71‑1 (MQ=255)
            gcgcGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCgtcg                                              <  1:1483129/71‑1 (MQ=255)
              gcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCgg                                                                 <  1:1268265/50‑1 (MQ=255)
              gcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCgtcggt                                            >  1:426200/1‑71 (MQ=255)
                gAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCgg                                                                 <  1:1547279/48‑1 (MQ=255)
                 aGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                                  <  1:39652/46‑1 (MQ=255)
                 aGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGc                                         <  1:1492563/71‑1 (MQ=255)
                      cAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                                  <  1:151285/41‑1 (MQ=25)
                      cAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGt                                    >  1:1668924/1‑71 (MQ=255)
                         cGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTAGCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTgg                                 >  1:1356362/1‑71 (MQ=255)
                             ttAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                             >  1:1244100/1‑71 (MQ=255)
                                  cAGGCGAAGCTGCTTTTGAGGACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGa                        <  1:1291301/71‑1 (MQ=255)
                                   aGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTg                                  >  1:395242/1‑60 (MQ=255)
                                    ggCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTgg                                 >  1:1378415/1‑60 (MQ=255)
                                    ggCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                      <  1:1638559/71‑1 (MQ=255)
                                     gCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                      >  1:1225972/1‑70 (MQ=255)
                                       gAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAAc                   >  1:221299/1‑71 (MQ=255)
                                          gctgctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAActc                 >  1:1753503/1‑70 (MQ=255)
                                              ctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGc                                         <  1:1669041/42‑1 (MQ=255)
                                              ctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACTATCAACTCTAt              <  1:685885/69‑1 (MQ=255)
                                              ctTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt              <  1:621605/69‑1 (MQ=255)
                                               ttttGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                      <  1:1663611/60‑1 (MQ=255)
                                                tttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACg                         >  1:1306160/1‑56 (MQ=255)
                                                 ttGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGAACAACTCTATTtctc         >  1:1341189/1‑71 (MQ=255)
                                                 ttGAGCACCACCCGCATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTtct          <  1:1640706/70‑1 (MQ=255)
                                                     gCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTcgcg      <  1:1470505/70‑1 (MQ=255)
                                                        cTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa  <  1:834093/71‑1 (MQ=255)
                                                         |                                                                     
CGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTA  >  W3110S.gb/71052‑71178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: