Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 704,757 G→A 31.2% R131H (CGC→CAC)  nagE → fused N‑acetyl glucosamine specific PTS enzyme IICBA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb704,7570GA31.2% 65.4 / 38.7 48R131H (CGC→CAC) nagEfused N‑acetyl glucosamine specific PTS enzyme IICBA components
Reads supporting (aligned to +/- strand):  ref base G (17/16);  new base A (7/8);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

GTGGCGCAGCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACG  >  W3110S.gb/704688‑704816
                                                                     |                                                           
gtggCGCAGCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGc                                                            <  1:18808/71‑1 (MQ=255)
    cgcAGCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTtg                                                        <  1:4509308/71‑1 (MQ=255)
    cgcAGCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTtg                                                        <  1:4605581/71‑1 (MQ=255)
        gCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGcc                                                    <  1:1683079/71‑1 (MQ=255)
         ccTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACg                                                             >  1:3184130/1‑61 (MQ=255)
         ccTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACg                                                             >  1:1376241/1‑61 (MQ=255)
         ccTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACg                                                             >  1:2665851/1‑61 (MQ=255)
         ccTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACg                                                             >  1:1684578/1‑61 (MQ=255)
         ccTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGcc                                                    <  1:3717884/70‑1 (MQ=255)
         ccTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCg                                                   <  1:3570225/71‑1 (MQ=255)
                ccGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACt                                                           <  1:806031/56‑1 (MQ=255)
                ccGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCa                                            <  1:3924374/71‑1 (MQ=255)
                 cGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCAc                                           >  1:5142431/1‑71 (MQ=255)
                   ttGGTCCTATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCg                                                   <  1:1373362/61‑1 (MQ=255)
                   ttGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCg                                                   <  1:3758749/61‑1 (MQ=255)
                   ttGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCg                                         >  1:4161154/1‑71 (MQ=255)
                        ccGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCGGAttc                                    >  1:1639830/1‑71 (MQ=255)
                         cGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGTCAAACGCTTTGTGCCgt                                                  >  1:1101123/1‑55 (MQ=255)
                         cGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCg                                         >  1:465899/1‑65 (MQ=255)
                          gATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGc                                                            <  1:1847348/45‑1 (MQ=255)
                            tatTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCt                                                           <  1:2485458/44‑1 (MQ=255)
                             atTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGc                                                            <  1:1787041/42‑1 (MQ=255)
                             atTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGa                                                  >  1:693595/1‑52 (MQ=255)
                             atTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGa                                                  >  1:3653492/1‑52 (MQ=255)
                             atTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGa                                                  >  1:1429319/1‑52 (MQ=255)
                              ttAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCttt                                                         <  1:1141622/44‑1 (MQ=255)
                               tAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCg                                                   >  1:1971369/1‑49 (MQ=255)
                               tAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGc                                              >  1:3009528/1‑54 (MQ=255)
                                aaaCTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCtt                                                          >  1:1489195/1‑41 (MQ=255)
                                aaaCTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCtt                                                          >  1:263427/1‑41 (MQ=255)
                                aaaCTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCtt                                                          >  1:209400/1‑41 (MQ=255)
                                 aaCTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCt                            <  1:3800049/70‑1 (MQ=255)
                                 aaCTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCGGATTCTTCTGCCt                            <  1:3180562/70‑1 (MQ=255)
                                   cTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGAtt                                                <  1:1060346/48‑1 (MQ=255)
                                   cTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCAc                                           <  1:3330389/53‑1 (MQ=255)
                                   cTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGt                         >  1:321426/1‑71 (MQ=255)
                                    tGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCAcc                                          <  1:304556/53‑1 (MQ=255)
                                     gCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCg                                         >  1:3041205/1‑53 (MQ=255)
                                     gCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCg                                         >  1:3238983/1‑53 (MQ=255)
                                              ccTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGa                                                  >  1:1062949/1‑35 (MQ=255)
                                              ccTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTg                        >  1:577820/1‑61 (MQ=255)
                                              ccTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCAc                                           >  1:46239/1‑42 (MQ=255)
                                                  aGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCAtttt          >  1:5140439/1‑71 (MQ=255)
                                                  aGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGCCCAtttt          <  1:229903/71‑1 (MQ=255)
                                                           ggcggcAAACGCTTTGTGCCGCTTGCCACCGGAttctt                                  >  1:4644921/1‑38 (MQ=255)
                                                           ggcggcAAACGCTTTGTGCCGATTGCCACCGGAttctt                                  >  1:4153359/1‑38 (MQ=255)
                                                           ggcggcAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACg  <  1:2851356/70‑1 (MQ=255)
                                                            gcggcAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGtt     >  1:2072704/1‑66 (MQ=255)
                                                            gcggcAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGtt     >  1:1122535/1‑66 (MQ=255)
                                                               gcAAACACTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTgg                    <  1:1083325/48‑1 (MQ=255)
                                                                  aaCGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTggc                   <  1:2811732/46‑1 (MQ=255)
                                                                  aaCGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTggc                   <  1:647135/46‑1 (MQ=255)
                                                                     gCTTTGTGCCGGTTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCAttttt         <  1:4141715/53‑1 (MQ=255)
                                                                     gCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCAttttt         <  1:4047081/53‑1 (MQ=255)
                                                                     |                                                           
GTGGCGCAGCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACG  >  W3110S.gb/704688‑704816

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: