Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,581 C→T 62.2% 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,5810CT62.2% 25.6 / 40.3 37V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (10/4);  new base T (16/7);  total (26/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

ATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCC  >  W3110S.gb/70522‑70643
                                                           |                                                              
atTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCggg                                                     >  1:416283/1‑71 (MQ=255)
 ttCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGt                                                    >  1:1864100/1‑71 (MQ=255)
   ctcAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGa                                                  >  1:1163717/1‑71 (MQ=255)
         ctcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATAtttt                                            >  1:1049349/1‑71 (MQ=255)
          tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGAc                                                         >  1:1725896/1‑57 (MQ=255)
            aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATAttt                                             <  1:523963/67‑1 (MQ=255)
            aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCt                                         <  1:875853/71‑1 (MQ=255)
               aTTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGtt                                      >  1:1794906/1‑71 (MQ=255)
               aTTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTg                                        >  1:415277/1‑69 (MQ=255)
                  cgcgGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTccc                                   >  1:1506064/1‑71 (MQ=255)
                         aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                              >  1:1190935/1‑69 (MQ=255)
                           ggggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGcc                                <  1:817227/65‑1 (MQ=255)
                               tggtgAAAAATCAGGGACGAGAATATGTTTGCCGACCg                                                       <  1:1222472/38‑1 (MQ=38)
                                ggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                              >  1:986961/1‑62 (MQ=255)
                                 gtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                    >  1:484479/1‑71 (MQ=255)
                                 gtgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCTCCAGgag                           >  1:84290/1‑64 (MQ=255)
                                  tgAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGat                                                 <  1:866607/41‑1 (MQ=255)
                                   gAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCt                                         >  1:822957/1‑48 (MQ=255)
                                    aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGt                                       <  1:2062128/49‑1 (MQ=255)
                                    aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                 >  1:1774629/1‑71 (MQ=255)
                                    aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatca                 >  1:1382873/1‑71 (MQ=255)
                                    aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatc                  >  1:1774966/1‑70 (MQ=255)
                                    aaaaaTCAGGGACGAGAATTTGTCTGCCGACCGGGTGatat                                               >  1:590347/1‑41 (MQ=255)
                                     aaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                >  1:303415/1‑71 (MQ=255)
                                        aTCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTAc             >  1:513304/1‑71 (MQ=255)
                                         tCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGtt                                      <  1:106557/45‑1 (MQ=255)
                                              gACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCa                               >  1:1074445/1‑47 (MQ=255)
                                                cGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgat                           <  1:204919/49‑2 (MQ=255)
                                                cGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgag                           >  1:526015/1‑49 (MQ=255)
                                                cGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGgag                           >  1:528645/1‑49 (MQ=255)
                                                   gaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtca                    >  1:2109754/1‑53 (MQ=255)
                                                   gaATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATcc  >  1:2200375/1‑71 (MQ=255)
                                                    aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGga                            <  1:1806472/44‑1 (MQ=255)
                                                    aaTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcat                   <  1:1878909/53‑1 (MQ=255)
                                                      tttgtttgCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                              <  1:495492/40‑1 (MQ=255)
                                                      tttGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAc                >  1:478873/1‑54 (MQ=255)
                                                        tGTTTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatc                  >  1:438849/1‑50 (MQ=255)
                                                        tGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAg                              <  1:608335/38‑1 (MQ=255)
                                                           |                                                              
ATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCC  >  W3110S.gb/70522‑70643

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: