Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,581 C→T 48.3% V65V (GTC→GTT araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb70,5810CT48.3% 12.2 / 129.7 87V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (32/13);  new base T (19/23);  total (51/36)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGA  >  W3110S.gb/70513‑70646
                                                                    |                                                                 
aaaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAtttgtttg                                                                 >  1:1086315/1‑71 (MQ=255)
 aaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                >  1:2479581/1‑71 (MQ=255)
 aaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGc                                                                >  1:742542/1‑71 (MQ=255)
   ggTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGc                                                                >  1:899989/1‑69 (MQ=255)
     ttATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGAc                                                            >  1:3827837/1‑71 (MQ=255)
      tataTTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCg                                                              <  1:4372163/68‑1 (MQ=255)
       atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCg                                                          >  1:1194530/1‑71 (MQ=255)
        tatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCgg                                                         >  1:4212310/1‑71 (MQ=255)
          ttCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGt                                                       <  1:2303495/71‑1 (MQ=255)
          ttCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGt                                                       <  1:2202885/71‑1 (MQ=255)
           tctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTc                                                                   >  1:2804429/1‑58 (MQ=255)
            ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGAc                                                            >  1:611210/1‑64 (MQ=255)
            ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGa                                                     >  1:2818491/1‑71 (MQ=255)
               aaTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGatat                                                  >  1:4066323/1‑71 (MQ=255)
                   tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTg                                              >  1:4126967/1‑71 (MQ=255)
                   tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTg                                              >  1:231509/1‑71 (MQ=255)
                   tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                              >  1:1415606/1‑71 (MQ=255)
                   tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                              >  1:1747619/1‑71 (MQ=255)
                   tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                              >  1:1409908/1‑71 (MQ=255)
                    cACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGc                                             >  1:154270/1‑71 (MQ=255)
                      ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCt                                                                  <  1:844950/48‑1 (MQ=255)
                      ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCt                                                                  <  1:1796947/48‑1 (MQ=255)
                      ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCg                                                          >  1:1249046/1‑56 (MQ=255)
                           cgcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGat                                                    >  1:3957765/1‑57 (MQ=255)
                            gcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTccc                                      <  1:3334361/70‑1 (MQ=255)
                              ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                   <  1:4734203/71‑1 (MQ=255)
                                tCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                 <  1:4004444/71‑1 (MQ=255)
                                tCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                 <  1:3378815/71‑1 (MQ=255)
                                tCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                 <  1:428794/71‑1 (MQ=255)
                                 cAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCg                                                          >  1:4785052/1‑45 (MQ=255)
                                     gggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagag                            <  1:3308493/71‑1 (MQ=255)
                                     gggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagag                            <  1:2815840/71‑1 (MQ=255)
                                        tggtgAAAAATCATGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                         <  1:548862/71‑1 (MQ=255)
                                        tggtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                 >  1:3637114/1‑63 (MQ=255)
                                          gtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtc                        <  1:2486193/70‑1 (MQ=255)
                                          gtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                       >  1:2821926/1‑71 (MQ=255)
                                           tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                           <  1:2336231/66‑1 (MQ=255)
                                           tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                      <  1:1072940/71‑1 (MQ=255)
                                           tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                      <  1:3542572/71‑1 (MQ=255)
                                           tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                      <  1:958812/71‑1 (MQ=255)
                                           tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                      >  1:4215317/1‑71 (MQ=255)
                                           tgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGc                                             <  1:2988693/48‑1 (MQ=255)
                                           tgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCa                                  >  1:2465749/1‑59 (MQ=255)
                                            gAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                       >  1:2506675/1‑69 (MQ=255)
                                            gAAAAATCAGGGACGAGAATTTGTCTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                      >  1:2606408/1‑70 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                   <  1:1714694/56‑1 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                         >  1:500004/1‑66 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                             <  1:2197593/62‑1 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                    >  1:2757251/1‑71 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                    >  1:505984/1‑71 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                    >  1:4330523/1‑71 (MQ=255)
                                             aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatc                     <  1:3582243/70‑1 (MQ=255)
                                               aaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                           >  1:1430075/1‑46 (MQ=255)
                                               aaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                           >  1:4097621/1‑46 (MQ=255)
                                                aaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCa                                  >  1:4770601/1‑54 (MQ=255)
                                                aaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAt                          <  1:1806364/62‑1 (MQ=255)
                                                aaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                 >  1:4671905/1‑71 (MQ=255)
                                                 aTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                    >  1:736740/1‑51 (MQ=255)
                                                 aTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                 >  1:2756761/1‑70 (MQ=255)
                                                 aTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                 >  1:4517035/1‑70 (MQ=255)
                                                  tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGc                                             <  1:3072124/41‑1 (MQ=255)
                                                  tCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg               >  1:67439/1‑71 (MQ=255)
                                                  tCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCt                                            <  1:3973179/42‑1 (MQ=255)
                                                  tCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCt                                            <  1:2014939/42‑1 (MQ=255)
                                                   cAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGtt                                         <  1:1506740/44‑1 (MQ=255)
                                                   cAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACgg              >  1:4204227/1‑71 (MQ=255)
                                                   cAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                 >  1:3432851/1‑52 (MQ=255)
                                                   cAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                 >  1:279090/1‑52 (MQ=255)
                                                    aGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTg                                              >  1:1502745/1‑38 (MQ=255)
                                                        aCGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGt                                          <  1:3852877/38‑1 (MQ=255)
                                                        aCGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatc                     <  1:3332367/59‑1 (MQ=255)
                                                          gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                   <  1:701827/43‑1 (MQ=255)
                                                          gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtc            <  1:3183129/66‑1 (MQ=255)
                                                          gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCa        >  1:2570468/1‑70 (MQ=255)
                                                          gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCAt       >  1:2052899/1‑71 (MQ=255)
                                                          gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCAt       >  1:4252426/1‑71 (MQ=255)
                                                           agaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATc      <  1:2591248/71‑1 (MQ=255)
                                                           agaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                    <  1:3811120/41‑1 (MQ=255)
                                                            gaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCa                                  >  1:4658804/1‑42 (MQ=255)
                                                            gaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATc      <  1:3233280/70‑1 (MQ=255)
                                                            gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgtc         >  1:2533267/1‑67 (MQ=255)
                                                            gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc     >  1:756115/1‑71 (MQ=255)
                                                            gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc     >  1:2501307/1‑71 (MQ=255)
                                                            gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc     >  1:2506985/1‑71 (MQ=255)
                                                             aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                 <  1:4705465/42‑1 (MQ=255)
                                                             aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACg               >  1:2859829/1‑60 (MQ=255)
                                                              aTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgt             <  1:765310/61‑1 (MQ=255)
                                                               tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                 <  1:2833157/40‑1 (MQ=255)
                                                               tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgagaga                           >  1:4276282/1‑46 (MQ=255)
                                                               tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                 >  1:1374777/1‑56 (MQ=255)
                                                               tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGa  <  1:4607510/71‑1 (MQ=255)
                                                               tttgtttgACGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTa                 >  1:1845539/1‑56 (MQ=255)
                                                               tttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                   >  1:2197296/1‑38 (MQ=255)
                                                                 tGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                    >  1:1375243/1‑51 (MQ=255)
                                                                    |                                                                 
AAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGA  >  W3110S.gb/70513‑70646

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: