Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,581 C→T 41.0% 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,5810CT41.0% 41.9 / 72.6 61V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (24/12);  new base T (13/12);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

AAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAG  >  W3110S.gb/70514‑70647
                                                                   |                                                                  
aaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGc                                                                 >  1:4232607/1‑71 (MQ=255)
 aGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGcc                                                                >  1:1912352/1‑71 (MQ=255)
  ggTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCt                                                                   <  1:494026/67‑1 (MQ=255)
      atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCg                                                           >  1:1999385/1‑71 (MQ=255)
      atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGa                                                              >  1:3075074/1‑68 (MQ=255)
      atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCg                                                           >  1:4002579/1‑71 (MQ=255)
      atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCg                                                           >  1:3648492/1‑71 (MQ=255)
       tatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCgg                                                          >  1:65406/1‑71 (MQ=255)
       tatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCg                                                           <  1:2709883/70‑1 (MQ=255)
         ttCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAtttgtttg                                                                  >  1:2018826/1‑61 (MQ=255)
          tctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGt                                                        >  1:3764322/1‑70 (MQ=255)
           ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGa                                                      >  1:1834618/1‑71 (MQ=255)
                 ctcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtttt                                                >  1:4172449/1‑71 (MQ=255)
                  tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                               >  1:78619/1‑71 (MQ=255)
                  tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                               >  1:1460316/1‑71 (MQ=255)
                  tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                               >  1:2567922/1‑71 (MQ=255)
                   cACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGc                                              >  1:4448574/1‑71 (MQ=255)
                    aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCt                                             >  1:3885100/1‑71 (MQ=255)
                    aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAttt                                                 <  1:2545916/67‑1 (MQ=255)
                     ccATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                            >  1:2470306/1‑71 (MQ=255)
                        ttCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTc                                         >  1:2922449/1‑71 (MQ=255)
                          cgcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTccc                                       >  1:2782514/1‑71 (MQ=255)
                           gcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCg                                      <  1:4530415/71‑1 (MQ=255)
                             ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                    <  1:3233533/71‑1 (MQ=255)
                              gTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGata                                                    >  1:4097260/1‑54 (MQ=255)
                               tCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGa                                                              <  1:3046507/43‑1 (MQ=255)
                               tCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                  >  1:4329723/1‑71 (MQ=255)
                               gcAGGGGGTGGGGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGata                                                    <  1:2951537/52‑1 (MQ=39)
                                 aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                <  1:3948649/71‑1 (MQ=255)
                                 aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTg                                               >  1:1490332/1‑56 (MQ=255)
                                 aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                  >  1:2451252/1‑69 (MQ=255)
                                  gggggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATAtttt                                                <  1:2125186/54‑1 (MQ=255)
                                   ggggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                              <  1:2220800/71‑1 (MQ=255)
                                   ggggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                                >  1:4713743/1‑69 (MQ=255)
                                       tggtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAtt                          <  1:2759515/71‑1 (MQ=255)
                                         gtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGatat                                                   <  1:3166044/44‑1 (MQ=255)
                                         gtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                        <  1:1009540/71‑1 (MQ=255)
                                          tgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGc                                              <  1:1132433/48‑1 (MQ=255)
                                            aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                     >  1:4445163/1‑71 (MQ=255)
                                            aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGATATtcatca                     <  1:3092981/71‑1 (MQ=255)
                                            aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                     >  1:1174818/1‑71 (MQ=255)
                                              aaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGtt                                          <  1:2619438/48‑1 (MQ=255)
                                               aaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                    <  1:815298/53‑1 (MQ=255)
                                               aaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTc                                         >  1:1028808/1‑48 (MQ=255)
                                                   aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgt              >  1:62753/1‑71 (MQ=255)
                                                     ggACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcg            >  1:3767000/1‑71 (MQ=255)
                                                      gACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                     <  1:1683277/45‑1 (MQ=255)
                                                       aCGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgt           >  1:932821/1‑70 (MQ=255)
                                                       aCGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgt           >  1:770123/1‑70 (MQ=255)
                                                         gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgaga                              <  1:2942242/49‑1 (MQ=255)
                                                         gagaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                     >  1:1760108/1‑58 (MQ=255)
                                                          agaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATc       >  1:932485/1‑71 (MQ=255)
                                                          agaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATc       >  1:1399277/1‑71 (MQ=255)
                                                           gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc      <  1:4449213/71‑1 (MQ=255)
                                                           gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc      >  1:541159/1‑71 (MQ=255)
                                                            aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCg     <  1:1350868/71‑1 (MQ=255)
                                                            aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCg     <  1:3322878/71‑1 (MQ=255)
                                                            aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                    <  1:1866819/40‑1 (MQ=255)
                                                            aaTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                    <  1:4248309/56‑1 (MQ=255)
                                                              tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                                  >  1:286748/1‑40 (MQ=255)
                                                               ttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCg                                      >  1:3548519/1‑35 (MQ=255)
                                                               ttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGc                                     >  1:2090079/1‑36 (MQ=255)
                                                                  tCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAg  >  1:3516586/1‑68 (MQ=255)
                                                                   |                                                                  
AAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAG  >  W3110S.gb/70514‑70647

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: