Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,342,709 G→A 7.4% Q146Q (CAG→CAA yciM → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,342,7090GA 7.4% 169.5 / 2.9 54Q146Q (CAG→CAAyciMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/29);  new base A (2/2);  total (23/31)
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.54e-01

CGCGGAAGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTG  >  W3110S.gb/1342646‑1342772
                                                               |                                                               
cgcgGAAGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTa                                                           >  1:34090/1‑70 (MQ=255)
    gAAGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAg                                                                 <  1:2434747/60‑1 (MQ=255)
     aaGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATc                                                     >  1:2090437/1‑71 (MQ=255)
      aGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAAt                                                      >  1:3895841/1‑69 (MQ=255)
      aGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCt                                                    <  1:141897/71‑1 (MQ=255)
      aGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCt                                                    <  1:2823892/71‑1 (MQ=255)
        aCATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTAc                                                  >  1:4324385/1‑71 (MQ=255)
           tGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCa                                                <  1:896241/70‑1 (MQ=255)
           tGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAg                                               <  1:3189471/71‑1 (MQ=255)
           tGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAACTGCTACaaa                                                       <  1:896431/63‑1 (MQ=255)
           tGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAACTGCTACaaa                                                       <  1:153579/63‑1 (MQ=255)
            gTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCt                                                    <  1:2493414/65‑1 (MQ=255)
               cAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCt                                            >  1:807238/1‑70 (MQ=255)
                aaTCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTAc                                          <  1:1520258/71‑1 (MQ=255)
                 aTCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTAcc                                         <  1:3320242/71‑1 (MQ=255)
                   cAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCGACAGTTGCTACAAATCTACCa                                                >  1:76924/1‑62 (MQ=255)
                   cAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGtt                                                               >  1:1100043/1‑47 (MQ=255)
                    aGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCagt                                      >  1:2060513/1‑71 (MQ=255)
                    aGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCag                                       >  1:3719054/1‑70 (MQ=255)
                     gCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCt                                                    <  1:3739672/56‑1 (MQ=255)
                     gCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCagtg                                     <  1:2695861/71‑1 (MQ=255)
                         aCCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCagtgagtg                                 >  1:2421506/1‑71 (MQ=255)
                         aCCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCagtgagtg                                 >  1:3409361/1‑71 (MQ=255)
                         aCCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAACTGCTACAAAtt                                                     >  1:694467/1‑50 (MQ=255)
                         aCCGATGAAACTGACATCCGCATTGGCGCGCTGCAACAACTGCTACAAAtt                                                     >  1:2091731/1‑50 (MQ=255)
                          ccGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTgg                                <  1:2975885/71‑1 (MQ=255)
                             aTGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACa                                                         <  1:2144865/43‑1 (MQ=255)
                             aTGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACa                                                         <  1:1689341/43‑1 (MQ=255)
                             aTGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTAc                                                  <  1:2609501/50‑1 (MQ=255)
                               gAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGa                            >  1:2903623/1‑70 (MQ=255)
                                  aCTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAg                         >  1:4309314/1‑70 (MQ=255)
                                      aCTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCaa                      <  1:2636543/69‑1 (MQ=255)
                                          ccGCATTGGCGCGCTTCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGa                  <  1:2352716/69‑1 (MQ=255)
                                          ccGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGa                  <  1:2830077/69‑1 (MQ=255)
                                          ccGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATg                >  1:3415498/1‑71 (MQ=255)
                                           cGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCaa                      >  1:1238857/1‑64 (MQ=255)
                                             cATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTg             <  1:1722090/71‑1 (MQ=255)
                                              aTTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAAtt                    <  1:3029930/63‑1 (MQ=255)
                                               ttGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGa                            <  1:184868/54‑1 (MQ=255)
                                               ttGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGa                            <  1:190208/54‑1 (MQ=255)
                                               ttGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGcc           >  1:2639153/1‑71 (MQ=255)
                                                tGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGa                  >  1:3471133/1‑63 (MQ=255)
                                                 ggCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCg          <  1:3532970/70‑1 (MQ=255)
                                                   cgcgcTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGaa        >  1:2058406/1‑70 (MQ=255)
                                                    gcgcTGCAACAGTTGCTACAAATCTACCAGGCTACCagtgag                                   <  1:3423154/42‑1 (MQ=255)
                                                    gcgcTGCAACAGTTGCTACAAATCTACCAGGCTACCagtgag                                   <  1:134694/42‑1 (MQ=255)
                                                      gcTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGcc    >  1:51430/1‑71 (MQ=255)
                                                       cTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGAt                 <  1:2070543/57‑1 (MQ=255)
                                                        tGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGAt                 <  1:4531701/56‑1 (MQ=255)
                                                         gCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAAc       <  1:2174230/65‑1 (MQ=255)
                                                          cAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGc            >  1:3792040/1‑59 (MQ=255)
                                                          cAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGcc    <  1:342410/67‑1 (MQ=255)
                                                          cAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTg  <  1:3174289/69‑1 (MQ=255)
                                                            aCAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAAt                     >  1:3651372/1‑48 (MQ=255)
                                                            aCAGTTGCTACAAATCTACAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCg          <  1:1589237/58‑1 (MQ=255)
                                                               |                                                               
CGCGGAAGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTG  >  W3110S.gb/1342646‑1342772

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: