Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 704,757 G→A 8.5% 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
*AP009048704,7570GA 8.5% 146.2 / 3.6 47R131H (CGC→CAC) nagEfused N‑acetyl glucosamine specific PTS enzyme IICBA components
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base A (2/2);  total (22/25)
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.65e-01

CAGCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCC  >  AP009048/704694‑704823
                                                               |                                                                  
cAGCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTtgtg                                                             >  1:1675346/1‑71 (MQ=255)
    cTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGa                                                         <  1:435892/71‑1 (MQ=255)
        aaCCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGc                                                     >  1:373038/1‑71 (MQ=255)
             ttGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCg                                                >  1:834101/1‑71 (MQ=255)
               ggTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGa                                              >  1:3271402/1‑71 (MQ=255)
                 tCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGcc                                                    <  1:2199239/63‑1 (MQ=255)
                 tCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCGGAtt                                            >  1:464031/1‑71 (MQ=255)
                   cGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGAttct                                          >  1:3381644/1‑71 (MQ=255)
                   cGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGAttct                                          >  1:2900113/1‑71 (MQ=255)
                    gATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGAttctt                                         <  1:4106410/71‑1 (MQ=255)
                    gATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGAttctt                                         >  1:4036416/1‑71 (MQ=255)
                     atatTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGAttcttc                                        >  1:2802358/1‑71 (MQ=255)
                       atTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTg                                      >  1:991750/1‑71 (MQ=255)
                          aaaCTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCt                                   >  1:2598891/1‑71 (MQ=255)
                             cTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGt                                >  1:4311340/1‑71 (MQ=255)
                              tGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGAttc                                           <  1:2319023/59‑1 (MQ=255)
                                ccGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCt                             >  1:4436335/1‑71 (MQ=255)
                                 cGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTg                            <  1:1686370/71‑1 (MQ=255)
                                  ggACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTg                            <  1:3787474/70‑1 (MQ=255)
                                  ggACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTg                            <  1:3782100/70‑1 (MQ=255)
                                   gACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGAt                                                        >  1:2927378/1‑41 (MQ=255)
                                   gACTTCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTggc                          <  1:4130035/71‑1 (MQ=255)
                                     cTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTtgcgg                        >  1:2674870/1‑71 (MQ=255)
                                      ttCCTGAGCTTCTTCGGCGGCAAACACTTTGTGCCGATTGCCACCg                                                <  1:911994/46‑1 (MQ=255)
                                           gAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCAttt                  >  1:2761401/1‑71 (MQ=255)
                                            aGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCAttt                  >  1:2902219/1‑70 (MQ=255)
                                            aGCTGCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCAttt                  <  1:2906072/70‑1 (MQ=255)
                                             gCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCAttttt                <  1:4288030/71‑1 (MQ=255)
                                              cttcttCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTg               <  1:1567666/71‑1 (MQ=255)
                                                 cttcGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTgg              <  1:1618341/69‑1 (MQ=255)
                                                 cttcGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTgg              <  1:2666487/69‑1 (MQ=255)
                                                  ttcGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCGTGGTGCTGGCGGCCATTTTTgg              <  1:3365516/68‑1 (MQ=255)
                                                  ttcGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTggc                          <  1:4277252/56‑1 (MQ=255)
                                                   tcGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTAc          >  1:3123824/1‑71 (MQ=255)
                                                    cggcggCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGc                                     <  1:588419/43‑1 (MQ=255)
                                                    cggcggCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGt             >  1:1632643/1‑67 (MQ=255)
                                                    cggcggCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACg         >  1:3927057/1‑71 (MQ=255)
                                                      gcggcAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGt             <  1:1442411/65‑1 (MQ=255)
                                                       cggcAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCt      <  1:1558918/71‑1 (MQ=255)
                                                        ggcAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCt      <  1:4158774/70‑1 (MQ=255)
                                                         gcAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTggcggca                      <  1:3761560/53‑2 (MQ=255)
                                                         gcAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTgg                           <  1:3393216/48‑1 (MQ=255)
                                                          cAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCt                                   <  1:1744271/39‑1 (MQ=255)
                                                           aaaCGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTggcgg                        >  1:3348821/1‑49 (MQ=255)
                                                           aaaCGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGgcc  <  1:868587/71‑1 (MQ=255)
                                                              cGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTgg                           >  1:4014980/1‑43 (MQ=255)
                                                              cGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCt      <  1:3774439/64‑1 (MQ=255)
                                                               |                                                                  
CAGCCTATAACCGTTGGTCCGATATTAAACTGCCGGACTTCCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCC  >  AP009048/704694‑704823

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: