Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,342,710 T→C 7.5% L147L (TTG→CTG)  yciM → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,342,7100TC 7.5% 137.0 / 3.6 53L147L (TTG→CTG) yciMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/29);  new base C (2/2);  total (22/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 = 8.60e-01

CGCGGAAGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCT  >  W3110S.gb/1342646‑1342780
                                                                |                                                                      
cgcgGAAGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTa                                                                   >  1:34090/1‑70 (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)
                                                                ttGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGa                                    <  1:984170/37‑1 (MQ=255)
                                                                ttGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCt  <  1:1391465/71‑1 (MQ=255)
                                                                |                                                                      
CGCGGAAGACATGTTCAATCAGCTGACCGATGAAACTGACTTCCGCATTGGCGCGCTGCAACAGTTGCTACAAATCTACCAGGCTACCAGTGAGTGGCAGAAAGCAATTGATGTTGCCGAACGCCTGGTGAAGCT  >  W3110S.gb/1342646‑1342780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: