Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,127,320 G→A 6.1% A188A (GCC→GCT wcaI ← predicted glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,127,3200GA 6.1% 193.2 / 2.3 66A188A (GCC→GCTwcaIpredicted glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (2/2);  total (34/32)
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.28e-01

AAAACGGGCAATTTCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCGCTACGTTCG  >  W3110S.gb/2127254‑2127386
                                                                  |                                                                  
aaaaCGGGCAATTTCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCtt                                                                <  1:3696735/71‑1 (MQ=255)
           tttCTGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATAGCTTTATTcatcatt                                                     <  1:1536878/71‑2 (MQ=255)
           tttCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTAtt                                                            <  1:3134967/64‑1 (MQ=255)
           tttCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTAtt                                                            <  1:1203628/64‑1 (MQ=255)
           tttCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTAtt                                                            <  1:762747/64‑1 (MQ=255)
           tttCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCatcatc                                                     >  1:1893115/1‑71 (MQ=255)
             tCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCatcatc                                                     >  1:2738547/1‑69 (MQ=255)
             tCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCatcatc                                                     >  1:1249812/1‑69 (MQ=255)
             tCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCg                                                    >  1:3370913/1‑70 (MQ=255)
             tCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCg                                                    >  1:701238/1‑70 (MQ=255)
             tCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCg                                                    >  1:2023058/1‑70 (MQ=255)
             tCCGACCAGTTGGGGAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGaa                                                  >  1:836940/1‑71 (MQ=255)
             ctcGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCg                                                    <  1:1013119/68‑1 (MQ=255)
                gACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCg                                                    <  1:130844/67‑1 (MQ=255)
                gACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCAACGCCTTTTTCGATGGCTTTATTCATCATCg                                                    <  1:1888212/67‑1 (MQ=255)
                       tGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAAcgc                                               >  1:821308/1‑65 (MQ=255)
                       tGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAAcgc                                               >  1:1998439/1‑65 (MQ=255)
                       tGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATc                                         >  1:335888/1‑71 (MQ=255)
                        ggggAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCg                                        >  1:1222052/1‑71 (MQ=255)
                         gggAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCg                                        <  1:226880/70‑1 (MQ=255)
                         gggAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCg                                        <  1:1256289/70‑1 (MQ=255)
                          ggAAGAAGATGACGTTTTCCGCCGCCACGCCTTTCTCGATGGCTTTATTCATCATCGAACGCGAAATCGTg                                      <  1:2103355/71‑1 (MQ=255)
                            aagaagATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGga                                    <  1:2998829/71‑1 (MQ=255)
                             agaagaTGACGTTTTCCGCCGCCACGCCTTTTTCGATg                                                                    >  1:886398/1‑38 (MQ=255)
                             agaagaTGACGTTTTCCGCCGCCACGCCTTTTTCGATg                                                                    >  1:890911/1‑38 (MQ=255)
                              gaagaTGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGga                                    <  1:1622359/69‑1 (MQ=255)
                               aagaTGACGTTTTCCGCCGCCACGCCTTTTTCGATAGCTTTATTCATCATTGAGCGCGAAATCGTGGAGAc                                 >  1:3047927/1‑71 (MQ=255)
                               aagaTGACGTTTTCCGCCGCCACGCCTTTTTCGATAGCTTTATTCATCATTGAGCGCGAAATCGTGGAGAc                                 >  1:1627096/1‑71 (MQ=255)
                                agaTGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTa                                                              >  1:1189693/1‑41 (MQ=255)
                                agaTGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACg                                <  1:210928/71‑1 (MQ=255)
                                  aTGACGTTTTCCGCCGCCACGCCTTTTTCGATAGCTTTa                                                              <  1:895715/39‑1 (MQ=255)
                                   tGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCtt                                                                <  1:2399447/36‑1 (MQ=255)
                                    gACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTAt                            >  1:113172/1‑71 (MQ=255)
                                     aCGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATc                           <  1:2772857/71‑1 (MQ=255)
                                      cGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCg                          <  1:1441299/71‑1 (MQ=255)
                                          ttCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGt                      >  1:633918/1‑71 (MQ=255)
                                           tCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGt                      <  1:2427481/70‑1 (MQ=255)
                                             cgccgcCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTa                    <  1:3601231/70‑1 (MQ=255)
                                                 gccACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAg               <  1:1586479/71‑1 (MQ=255)
                                                   cACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTc             <  1:36195/71‑1 (MQ=255)
                                                    aCGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTcc            <  1:2126487/71‑1 (MQ=255)
                                                    aCGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTcc            >  1:1681011/1‑71 (MQ=255)
                                                         tttttCGATGGCTTTATTCATCATCGAACGCGTAATCGTGGAGACg                                <  1:3279211/46‑1 (MQ=255)
                                                         tttttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGga                                    <  1:1040343/42‑1 (MQ=255)
                                                         tttttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGga                                    <  1:2100174/42‑1 (MQ=255)
                                                         tttttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGAc                                 <  1:3546108/45‑1 (MQ=255)
                                                         tttttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGt                      <  1:947510/56‑1 (MQ=255)
                                                           tttCGATTGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGtt                              >  1:349433/1‑46 (MQ=255)
                                                           tttCGATTGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGtt                              >  1:2256123/1‑46 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCg                                        >  1:2940538/1‑36 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCg                                        >  1:688488/1‑36 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGa                         >  1:3014834/1‑51 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTa                    >  1:2840733/1‑56 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATg                  <  1:3798252/58‑1 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCa                >  1:2987136/1‑60 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCa                >  1:1717177/1‑60 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCg           >  1:1105452/1‑65 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCg           >  1:3200745/1‑65 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCGCTACg      <  1:671017/70‑1 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCGCTACg      <  1:246363/70‑1 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCGCTACGt     >  1:3651317/1‑71 (MQ=255)
                                                           tttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCGCTACGt     >  1:1588828/1‑71 (MQ=255)
                                                           tttCGACGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGa                         >  1:497854/1‑51 (MQ=255)
                                                            ttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGga                                    >  1:583586/1‑39 (MQ=255)
                                                            ttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGa                         >  1:2666767/1‑50 (MQ=255)
                                                            ttCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGa                         >  1:649948/1‑50 (MQ=255)
                                                             tCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCGCTACGtt    >  1:2568970/1‑70 (MQ=255)
                                                               gATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCGCTACGTTCg  <  1:2884425/70‑1 (MQ=255)
                                                                  |                                                                  
AAAACGGGCAATTTCCGACCAGTTGGGGAAGAAGATGACGTTTTCCGCCGCCACGCCTTTTTCGATGGCTTTATTCATCATCGAACGCGAAATCGTGGAGACGTTATCGACGTTATGCAGTCCGCTACGTTCG  >  W3110S.gb/2127254‑2127386

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