Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,109 T→C 45.5% 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,1090TC45.5% 3.7 / 43.7 33T241T (ACT→ACCaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (9/9);  new base C (7/8);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAG  >  W3110S.gb/71039‑71176
                                                                      |                                                                   
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)
     gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                                <  1:1667996/71‑1 (MQ=255)
     gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCg                                                                <  1:2450344/71‑1 (MQ=255)
     gggATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCg                                                                <  1:2407224/71‑1 (MQ=255)
                   aaGCTGGCGCGAGGACCAACGTATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCa                                                  >  1:2341636/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)
                   aaGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCa                                                  >  1:1422799/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)
                                            agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGtttt                           <  1:131161/66‑1 (MQ=255)
                                            agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa                         >  1:1612797/4‑71 (MQ=255)
                                            agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa                         <  1:173070/68‑1 (MQ=255)
                                            agccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa                         >  1:575689/4‑71 (MQ=255)
                                            aGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGa                         >  1:918248/1‑71 (MQ=255)
                                              ccAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACg                       <  1:906767/70‑1 (MQ=255)
                                                 ggCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                    <  1:1230924/71‑1 (MQ=255)
                                                  gCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                    >  1:1636325/1‑70 (MQ=255)
                                                  gCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                    >  1:1194476/1‑70 (MQ=255)
                                                  gCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATc                    >  1:2548575/1‑70 (MQ=255)
                                                   cGAAGCTGCTTTTGATCACCACCCGGATGCCTATCGCCACCGTCTGTCGct                                      <  1:1476818/51‑2 (MQ=255)
                                                        ctgctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTa             <  1:2341048/71‑1 (MQ=255)
                                                        ctgctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTa             <  1:1687699/71‑1 (MQ=255)
                                                           ctTTTGAGCATCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAtt           <  1:1323404/70‑1 (MQ=255)
                                                           ctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt            <  1:1956097/69‑1 (MQ=255)
                                                           ctTTTGAGCACTACACGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAt            <  1:2091573/69‑1 (MQ=255)
                                                           ctTTTGAGCACCACCCGTATGCCTATCGCCACCGTCGGTCGCaa                                     >  1:1168515/1‑44 (MQ=38)
                                                                   caccacCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTcgcggg  <  1:1094333/71‑3 (MQ=255)
                                                                      |                                                                   
AGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAG  >  W3110S.gb/71039‑71176

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: