Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,526,720 G→A 8.5% S78N (AGC→AAC)  PP_4860 → mannose‑6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,526,7200GA 8.5% 130.0 / 2.3 47S78N (AGC→AAC) PP_4860mannose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (2/2);  total (24/23)
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.63e-01

TCGCCGAACGCATCGTGCCGCTGTGGCAAGGCCCGGGCTGGAACGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGACCGGCGCAAGGACCT  >  NC_002947/5526573‑5526864
                                                                                                                                                   |                                                                                                                                                
tCGCCGAACGCATCGTGCCGCTGTGGCAAGGCCCGGGCTGGAACGCCGACATGGCACTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTcag                                                                                                                                                  >  3:72548/1‑148 (MQ=255)
 cGCCGAACGCATCGTGCCGCTGTGGCAAGGCCCGGGCTGGAACGCCGACATGGCCCCGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTcag                                                                                                                                                  <  6:297188/147‑1 (MQ=255)
   ccGAACGCATCGTGCCGCTGTGGCAAGGCCCGGGCTGGAACGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCagcagc                                                                                                                                              >  8:126789/1‑149 (MQ=255)
    cGAACGCATCGTGCCGCTGTGGCAAGGCCCGGGCTGGAACGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGcc                                                                                                                                             >  8:125987/1‑149 (MQ=255)
          cATCGTGCCGCTGTGGCAAGGCCCGGGCTGGAACGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCg                                                                                                                                       <  7:6812/149‑1 (MQ=255)
            tCGTGCCGCTGTGGCAAGGCCCGGGCTGGAACGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGc                                                                                                                                     <  1:17256/149‑1 (MQ=255)
                                ccGGGCTGGAACGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAAcg                                                                                                                 >  5:79132/1‑149 (MQ=255)
                                 cGGGCTGGAACGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAAcgc                                                                                                                >  7:91460/1‑149 (MQ=255)
                                          aCGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGcc                                                                                                       <  8:233998/149‑1 (MQ=255)
                                           cGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGccc                                                                                                      <  1:248857/149‑1 (MQ=255)
                                            gCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCACCCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCCGCCGGGCGCCGCCGAACGCGCCGCAGCCCt                                                                                                     >  4:79876/1‑149 (MQ=255)
                                            gCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCt                                                                                                     <  7:126789/149‑1 (MQ=255)
                                                 cATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTtc                                                                                                <  6:79132/149‑1 (MQ=255)
                                                    ggCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCgccgcc                                                                                                                      <  8:147264/124‑1 (MQ=255)
                                                    ggCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCgccgcc                                                                                                                      >  7:147264/1‑124 (MQ=255)
                                                      ccctgccctACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAACAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTtcgctcg                                                                                           >  8:198726/1‑149 (MQ=255)
                                                      ccctgccctACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAACAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTtcgctc                                                                                            >  8:198724/1‑148 (MQ=255)
                                                         tgccctACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTg                                                                                        >  1:250900/1‑149 (MQ=255)
                                                                      gcTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTcc                                                                           <  7:125987/149‑1 (MQ=255)
                                                                         ggACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCAc                                                                         >  7:252724/1‑148 (MQ=255)
                                                                             gCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGc                                                                    >  8:196694/1‑149 (MQ=255)
                                                                               ccAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGcc                                                                   >  5:246306/1‑148 (MQ=255)
                                                                                   cACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCa                                                              >  2:240930/1‑149 (MQ=255)
                                                                                       aaCCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACggg                                                          <  8:91460/149‑1 (MQ=255)
                                                                                        aCCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACgggg                                                         >  6:67438/1‑149 (MQ=255)
                                                                                              gCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTgg                                                    <  4:72548/148‑1 (MQ=255)
                                                                                                    gCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTc                                                                            >  2:276502/1‑118 (MQ=255)
                                                                                                    gCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTc                                                                            <  1:276502/118‑1 (MQ=255)
                                                                                                       gcgcTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCa                                          <  7:196694/149‑1 (MQ=255)
                                                                                                           tACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGa                                      >  8:80096/1‑149 (MQ=255)
                                                                                                                 gCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCa                                <  1:240930/149‑1 (MQ=255)
                                                                                                                  ccATGGCCTGTGCACGCCAGCTGTACGTGTTCAACAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCaa                               <  7:198724/149‑1 (MQ=255)
                                                                                                                  ccATGGCCTGTGCACGCCAGCTGTACGTGTTCAACAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCaa                               <  7:198726/149‑1 (MQ=255)
                                                                                                                   cATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAg                              <  8:252724/149‑1 (MQ=255)
                                                                                                                    aTGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCACgg                             >  2:255623/1‑149 (MQ=255)
                                                                                                                         cTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAga                        >  7:116664/1‑148 (MQ=255)
                                                                                                                            tgCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGcc                       <  2:250900/147‑1 (MQ=255)
                                                                                                                            tgCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCg                      <  6:246306/148‑1 (MQ=255)
                                                                                                                               aCGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGc                                                               <  3:183201/104‑1 (MQ=255)
                                                                                                                               aCGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGc                                                               >  4:183201/1‑104 (MQ=255)
                                                                                                                               aCGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTg                   >  6:83955/1‑148 (MQ=255)
                                                                                                                                       tGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGAc                <  4:77368/143‑1 (MQ=255)
                                                                                                                                       tGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGAc                >  3:77368/1‑143 (MQ=255)
                                                                                                                                        gTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCaa                          >  1:183024/1‑132 (MQ=255)
                                                                                                                                        gTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCaa                          <  2:183024/132‑1 (MQ=255)
                                                                                                                                          aCGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGACCGGCGCAAg       >  1:14811/1‑149 (MQ=255)
                                                                                                                                          aCGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGACCGGCGCAAg       >  7:124325/1‑149 (MQ=255)
                                                                                                                                            gtgtTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGAcc               >  2:35208/1‑139 (MQ=255)
                                                                                                                                            gtgtTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGAcc               <  1:35208/139‑1 (MQ=255)
                                                                                                                                               ttCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGACCGGCGCAAGGACCt  <  3:79876/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCGCCGAACGCATCGTGCCGCTGTGGCAAGGCCCGGGCTGGAACGCCGACATGGCCCTGCCCTACGAAGCGCTGGACGCCCAGCACCAACCGTTGCCGGTGCAGCGCTACCGGGCCATGGCCTGTGCACGCCAGCTGTACGTGTTCAGCAGCCGTATCGGCCAGCCGGGCGCCGCCGAACGCGCCGCAGCCCTGTTTCGCTCGCTGCAAAAGCACTTCCACGATGCCGAGCACGGGGGCTGGTTCTACAGCATCGACGCCCAAGGCAAGCCGCTGGACCGGCGCAAGGACCT  >  NC_002947/5526573‑5526864

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