Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,164,007 C→A 7.6% T10T (ACC→ACA yheO → predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,164,0070CA 7.6% 211.0 / 2.3 65T10T (ACC→ACAyheOpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (44/16);  new base A (2/3);  total (46/19)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

AAGTTAACGACAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCTGAAATCCTAC  >  W3110S.gb/4163939‑4164073
                                                                    |                                                                  
aaGTTAACGACAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAg                                                                  >  1:1137872/1‑71 (MQ=255)
    tAACGACAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTg                                                                >  1:2068589/1‑69 (MQ=255)
     aaCGACAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCTCTTTTAACCAACGAAACCa                                                                   >  1:3987608/1‑65 (MQ=255)
     aaCGACAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAAcc                                                                    >  1:1364960/1‑64 (MQ=255)
     aaCGACAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGt                                                             >  1:3561885/1‑71 (MQ=255)
       cGACAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGt                                                             >  1:3041349/1‑69 (MQ=255)
         aCAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGa                                                         >  1:176305/1‑71 (MQ=255)
         aCAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACAAGTGAGTTGGa                                                         <  1:905334/71‑1 (MQ=255)
            ggCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGAtt                                                       >  1:4245375/1‑70 (MQ=255)
                 cTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGAtt                                                       <  1:2406808/65‑1 (MQ=255)
                 cTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACt                                                   <  1:3554938/69‑1 (MQ=255)
                 cTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTg                                                  <  1:1067274/70‑1 (MQ=255)
                 cTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACAAGTGAGTTGGATTTACTg                                                  >  1:2825467/1‑70 (MQ=255)
                  tGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGa                                                >  1:3519496/1‑71 (MQ=255)
                  tGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGa                                                >  1:2553089/1‑71 (MQ=255)
                   gAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAAcc                                                                    >  1:2367248/1‑50 (MQ=255)
                   gAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGAt                                               >  1:1902444/1‑71 (MQ=255)
                   gAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGAt                                               >  1:924950/1‑71 (MQ=255)
                    aaaaGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGa                                                <  1:2989365/69‑1 (MQ=255)
                    aaaaGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGCATTTACTGGa                                                >  1:1722028/1‑69 (MQ=255)
                    aaaaGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACAAGTGAGTTGGATTTACTGGAt                                               >  1:590173/1‑70 (MQ=255)
                         gAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAAcc                                                                    <  1:1173665/44‑1 (MQ=255)
                           gtgtTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTcc                                       >  1:2318889/1‑71 (MQ=255)
                            tgtTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCa                                             <  1:2355362/64‑1 (MQ=255)
                             gtTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTAc                                                    >  1:391234/1‑56 (MQ=255)
                             gtTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCtt                                     >  1:896473/1‑71 (MQ=255)
                              ttttttttCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCttt                                    >  1:699101/1‑71 (MQ=255)
                                ttttttCATGTCCAGGTCGCTTTTAACCAACGAAACCAGt                                                                 >  1:64818/1‑40 (MQ=255)
                                ttttttCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGt                                         <  1:789329/64‑1 (MQ=255)
                                ttttttCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTc                                   >  1:2895242/1‑70 (MQ=255)
                                  ttttCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGAttt                                                      <  1:2420990/49‑1 (MQ=255)
                                   tttCATGTCCAGGTCGCTTTTAACCAACGAAACCAGt                                                                 >  1:2155756/1‑37 (MQ=255)
                                   tttCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGAtt                                                       >  1:1263126/1‑47 (MQ=255)
                                    ttCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGAc                                >  1:1073750/1‑69 (MQ=255)
                                     tCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCtt                                     >  1:2205504/1‑63 (MQ=255)
                                       aTGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATc                                              >  1:2891627/1‑52 (MQ=255)
                                       aTGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTcc                                       >  1:3732281/1‑59 (MQ=255)
                                               gTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCa                                             <  1:2233737/45‑1 (MQ=255)
                                               gTCGCTTTTAACCAACGAAACCAGTGAGTTGGACTTACTGGATCa                                             <  1:611000/45‑1 (MQ=255)
                                                 cGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCtt                                     >  1:1124321/1‑51 (MQ=255)
                                                 cGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGCCCGATTTTGa                  >  1:2254081/1‑70 (MQ=255)
                                                 cGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGa                  >  1:4233396/1‑70 (MQ=255)
                                                  gCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATc                                              >  1:3635373/1‑41 (MQ=255)
                                                  gCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGata                >  1:2564890/1‑71 (MQ=255)
                                                  gCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGata                >  1:1556459/1‑71 (MQ=255)
                                                   cTTTTAACCAACGAAACAAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGAtttt                    <  1:223919/66‑1 (MQ=255)
                                                    ttttAACCAACGAAACCAGTGAGTTGGATTTACTg                                                  <  1:1670997/35‑1 (MQ=255)
                                                    ttttAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGa                  <  1:3429401/67‑1 (MQ=255)
                                                     tttAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTc             >  1:99843/1‑71 (MQ=255)
                                                     tttAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTc             >  1:3049554/1‑71 (MQ=255)
                                                      ttAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGa                                 >  1:1183810/1‑50 (MQ=255)
                                                      ttAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAg                             >  1:2608726/1‑54 (MQ=255)
                                                      ttAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATAtt              <  1:244165/69‑1 (MQ=255)
                                                      ttAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCt            >  1:2104362/1‑71 (MQ=255)
                                                      ttAACCAACGAAACAAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATAtt              <  1:38190/69‑1 (MQ=255)
                                                       tAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGAc                           >  1:850850/1‑55 (MQ=255)
                                                          ccAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCg                                  <  1:439619/45‑1 (MQ=255)
                                                          ccAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGAcc                               >  1:2022999/1‑48 (MQ=255)
                                                          ccAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGAc                                >  1:2064268/1‑47 (MQ=255)
                                                           cAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTcc                                       <  1:3979619/39‑1 (MQ=255)
                                                           cAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCTGaaa        >  1:182760/1‑70 (MQ=255)
                                                             aCGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCTGAAATcc     >  1:478422/1‑71 (MQ=255)
                                                             aCGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCTGAAATcc     >  1:3590484/1‑71 (MQ=255)
                                                             aCGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCTGAAATcc     >  1:798533/1‑71 (MQ=255)
                                                              cGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGAtt                      <  1:3748835/53‑1 (MQ=255)
                                                              cGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCTGAAATcc     >  1:3169499/1‑70 (MQ=255)
                                                                aaaCCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCTGAAATCCTAc  >  1:1623246/1‑71 (MQ=255)
                                                                  aCCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGata                <  1:1144172/55‑1 (MQ=255)
                                                                    |                                                                  
AAGTTAACGACAGGCTCCTGAAAAGGAGTGTTTTTTTTCATGTCCAGGTCGCTTTTAACCAACGAAACCAGTGAGTTGGATTTACTGGATCAACGTCCTTTCGACCAGACCGATTTTGATATTCTGAAATCCTAC  >  W3110S.gb/4163939‑4164073

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: