Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,135,190 T→A 6.9% T115T (ACA→ACT wcaA ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,135,1900TA 6.9% 188.3 / 4.8 72T115T (ACA→ACTwcaApredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (23/44);  new base A (3/2);  total (26/46)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCA  >  W3110S.gb/2135124‑2135259
                                                                  |                                                                     
cAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGa                                                                   >  1:2425311/1‑71 (MQ=255)
cAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGa                                                                   >  1:2961435/1‑71 (MQ=255)
 aaGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGAc                                                                  >  1:1329944/1‑71 (MQ=255)
  aGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGAcc                                                                 <  1:2086894/71‑1 (MQ=255)
  aGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGAcc                                                                 <  1:2381118/71‑1 (MQ=255)
  aGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGAcc                                                                 <  1:3165853/71‑1 (MQ=255)
  aGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGAcc                                                                 <  1:3810651/71‑1 (MQ=255)
   gCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCa                                                                <  1:881718/71‑1 (MQ=255)
   gCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCa                                                                <  1:1601400/71‑1 (MQ=255)
        ccGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGc                                                           <  1:1056566/71‑1 (MQ=255)
            ttGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:3179371/71‑1 (MQ=255)
            ttGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:587598/71‑1 (MQ=255)
            ttGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:2822031/71‑1 (MQ=255)
             tGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:1741630/70‑1 (MQ=255)
              gggAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:2596861/69‑1 (MQ=255)
               ggAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTAt                                                    <  1:3598729/71‑1 (MQ=255)
               ggAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTAt                                                    <  1:1691423/71‑1 (MQ=255)
                 aaTAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:1413761/66‑1 (MQ=255)
                 aaTAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:1205747/66‑1 (MQ=255)
                 aaTAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:1507449/66‑1 (MQ=255)
                 aaTAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCATGTGTGACCAGTTGCTGtt                                                       <  1:3312913/66‑1 (MQ=255)
                 aaTAGACTTCGCACTGGCAAACGTAGTCGTTAGCGTACAAAAAGGC‑‑GTGTGACCAGTTGCTGtt                                                       <  1:1481309/64‑1 (MQ=255)
                   tAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGc                                                >  1:46591/1‑71 (MQ=255)
                     gACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCgtgt                                                                       >  1:711380/1‑46 (MQ=255)
                      aCTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAg                                             <  1:1315291/71‑1 (MQ=255)
                        ttCGCCCTGGCAAACGTAGTCGTTGGCGTACAAAAAGGCGTGTGTGAcc                                                                 <  1:1129520/49‑1 (MQ=255)
                        ttCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCgtgtgt                                                                     <  1:2417016/45‑1 (MQ=255)
                        ttCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTg                                                            <  1:697765/54‑1 (MQ=255)
                        ttCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAgg                                            >  1:969533/1‑70 (MQ=255)
                          cGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGa                                           <  1:3368841/69‑1 (MQ=255)
                          cGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGa                                           <  1:615749/69‑1 (MQ=255)
                           gCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGa                                        >  1:1355289/1‑71 (MQ=255)
                              cTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGc                                                           <  1:3840031/49‑1 (MQ=255)
                              cTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGc                                                           <  1:338681/49‑1 (MQ=255)
                               tGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGttt                                                      <  1:2055785/53‑1 (MQ=255)
                                ggCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTa                                                     >  1:2201219/1‑53 (MQ=255)
                                 gCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAgg                                            <  1:1058865/61‑1 (MQ=255)
                                 gCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCa                                  <  1:1113332/71‑1 (MQ=255)
                                 gCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCa                                  <  1:679623/71‑1 (MQ=255)
                                     aCGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACg                              <  1:2194685/71‑1 (MQ=255)
                                      cGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGa                                        >  1:1819417/1‑60 (MQ=255)
                                      cGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACgg                             >  1:1485782/1‑71 (MQ=255)
                                       gTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGt                            <  1:2536052/71‑1 (MQ=255)
                                        tAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTg                                                         <  1:1322025/41‑1 (MQ=255)
                                        tAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGttt                                                      <  1:1578229/44‑1 (MQ=255)
                                          gTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTg                          <  1:1049570/70‑1 (MQ=255)
                                          gTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTg                          <  1:334401/70‑1 (MQ=255)
                                          gTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTg                          <  1:2982058/70‑1 (MQ=255)
                                          gTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTg                          <  1:2821969/70‑1 (MQ=255)
                                          gTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTg                          <  1:2601216/70‑1 (MQ=255)
                                           tCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTggg                        >  1:139457/1‑71 (MQ=255)
                                           tCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTggg                        >  1:13755/1‑71 (MQ=255)
                                           tCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTgg                         >  1:1631598/1‑70 (MQ=255)
                                            cGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGtt                                                       <  1:1440080/39‑1 (MQ=255)
                                             gTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACgg                             <  1:445807/64‑1 (MQ=255)
                                              ttAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTa                                                     >  1:1455684/1‑39 (MQ=255)
                                              ttAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTa                                                     >  1:679079/1‑39 (MQ=255)
                                              ttAGCGTACAAAAAGGCGTGAGTGACCAGTTGCTGTTTa                                                     >  1:2481801/1‑39 (MQ=255)
                                               tAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGgtgt                     >  1:1080462/1‑70 (MQ=255)
                                                aGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTc                    <  1:299384/70‑1 (MQ=255)
                                                aGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTc                    >  1:3061220/1‑70 (MQ=255)
                                                aGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTc                    >  1:2963329/1‑70 (MQ=255)
                                                aGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGACGAGGAAGACGCTCAGACGGTTGGGTGTc                    >  1:805988/1‑70 (MQ=255)
                                                    tACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtc               >  1:1271341/1‑71 (MQ=255)
                                                    tACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtc               <  1:3221442/71‑1 (MQ=255)
                                                      cAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt             >  1:3481438/1‑71 (MQ=255)
                                                      cAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt             >  1:1888866/1‑71 (MQ=255)
                                                      cAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgt             >  1:1929333/1‑71 (MQ=255)
                                                      cAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcg              >  1:468884/1‑70 (MQ=255)
                                                        aaaaGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtca           <  1:2922105/71‑1 (MQ=255)
                                                        aaaaGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTCCATtcatcg           <  1:48072/71‑2 (MQ=255)
                                                        aaaaGGCGTGAGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGGCGGTTGGGTGTCCATtcatcg           <  1:2219764/71‑2 (MQ=255)
                                                          aaGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGt                            >  1:1553456/1‑52 (MQ=255)
                                                          aaGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGt                            >  1:2676393/1‑52 (MQ=255)
                                                                 gtgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTcgtcatcgtcat    >  1:3818907/1‑69 (MQ=255)
                                                                  tgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCa                  >  1:3771630/1‑54 (MQ=255)
                                                                  tgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATtcgtca           <  1:3489899/61‑1 (MQ=255)
                                                                  tgtgACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGtcatca  >  1:2102978/1‑70 (MQ=255)
                                                                  |                                                                     
CAAGCTTGCCGGTTGGGAATAGACTTCGCCCTGGCAAACGTAGTCGTTAGCGTACAAAAAGGCGTGTGTGACCAGTTGCTGTTTATGGGCGAGGAAGACGCTCAGACGGTTGGGTGTCCATTCGTCATCGTCATCA  >  W3110S.gb/2135124‑2135259

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: