Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,343,619 C→A 8.2% L287M (CTG→ATG)  pgi‑II → glucose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,343,6190CA 8.2% 117.5 / 4.4 49L287M (CTG→ATG) pgi‑IIglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (29/16);  new base A (2/2);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

GCCGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACA  >  NC_002947/5343472‑5343766
                                                                                                                                                   |                                                                                                                                                   
gCCGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCt                                                                                                                                                    >  5:302753/1‑149 (MQ=34)
        gTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCg                                                                                                                                            >  3:289704/1‑149 (MQ=34)
           ccGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCa                                                                                                                                         >  1:197757/1‑149 (MQ=34)
           ccGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCa                                                                                                                                         >  1:34427/1‑149 (MQ=34)
                 cTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGcc                                                                                                                                    >  8:126443/1‑148 (MQ=34)
                  ttCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCaa                                                                                                                                  >  1:346282/1‑149 (MQ=34)
                                  gCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgctg                                                                                                                  >  7:144049/1‑149 (MQ=34)
                                          gCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGgcgc                                                                                                          >  5:255831/1‑149 (MQ=37)
                                           ccGCTGGTGTCGGGTTGGGATCCCGCTAACAAAACACTTCCCCATTGGGGGAGGGGGGGGGGGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGcc                                                                                                         <  4:289704/149‑1 (MQ=11)
                                                ggtttAGCGTGCGTTAACCGCAAATAAAACATCTTCCCATTGGGGGACGGGGGGGGCGGCCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTacac                                                                                                    <  7:43186/145‑1 (MQ=11)
                                                  tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCa                                                                                                  >  6:96001/1‑149 (MQ=37)
                                                        gTTCGTTACCCGCTAACAAAAAATCTCCCAATTGTGGGAGTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGAcc                                                                                            <  2:34427/149‑1 (MQ=11)
                                                                  cgcgAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTc                                                                                    >  6:357016/1‑147 (MQ=37)
                                                                        gaaAACATCTTCCCAATGAGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGACGGCCATCGGCCTGCCGATCGCCCAGGCCATCGGCACTGCCAACTTCAAGGAACTCCTGTCCGGCGCCTACACAATGGACCAGCAATTc                                                                                    >  5:353344/1‑141 (MQ=11)
                                                                           aaCATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCAcc                                                                         >  8:246432/1‑149 (MQ=37)
                                                                           aaCATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCAcc                                                                         >  8:137636/1‑149 (MQ=37)
                                                                              aTCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGtt                                                                      >  3:313556/1‑149 (MQ=37)
                                                                                       aTGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACaa                                                                >  7:384620/1‑146 (MQ=34)
                                                                                        tgtgGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGaa                                                             >  1:398438/1‑148 (MQ=34)
                                                                                           cggcgggggggggCGGCCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATg                                                         <  4:228482/143‑1 (MQ=255)
                                                                                            ggACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCTCTTCCAGACCGCACCGTTCGACAAGAACATGc                                                        >  7:321162/1‑149 (MQ=21)
                                                                                            ggACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCCCCGTTCGACAAGAACATGc                                                        >  8:136068/1‑149 (MQ=21)
                                                                                              aCTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                      <  5:357016/149‑1 (MQ=34)
                                                                                               tggggtgggGGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                      <  4:313556/146‑1 (MQ=17)
                                                                                                  ggtgggCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGc                                                  >  8:241042/1‑149 (MQ=34)
                                                                                                    tgggCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTa                                                >  2:166321/1‑149 (MQ=34)
                                                                                                     gggcggccGGTACTCGCTGTGGTCGGCGATCGGCCTGCCGATCGCCCTGGCCATCGGGACTGGCATCTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACGGCACCGTTCGACAAGAACATGCCGGTGCTAc                                               <  8:144049/149‑1 (MQ=11)
                                                                                                          ggcgCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGGCAAGAACATGCCGGTGCTACTGGcc                                          <  1:228875/149‑1 (MQ=17)
                                                                                                               tACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCctgct                                     >  3:51770/1‑149 (MQ=32)
                                                                                                                 cTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGcatgttc                                     <  1:166321/147‑7 (MQ=17)
                                                                                                                  tCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTggg                                  >  2:354958/1‑149 (MQ=32)
                                                                                                                    gCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGc                                           <  1:119932/138‑1 (MQ=32)
                                                                                                                    gCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGc                                           >  2:119932/1‑138 (MQ=32)
                                                                                                                    gCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCg                                >  3:109468/1‑149 (MQ=32)
                                                                                                                        tgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtgt                             <  5:96001/148‑1 (MQ=32)
                                                                                                                        tgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtgt                             <  7:137636/148‑1 (MQ=32)
                                                                                                                          cagggCGATCGGCAGGCCGATGGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTgg                           >  7:83519/4‑148 (MQ=11)
                                                                                                                             ggccgatggCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                       <  8:349691/142‑1 (MQ=17)
                                                                                                                             ggCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGtgtg                            >  3:147358/1‑144 (MQ=32)
                                                                                                                             ggCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGtgtg                            <  4:147358/144‑1 (MQ=32)
                                                                                                                               ccATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATgg                      >  2:172347/1‑148 (MQ=32)
                                                                                                                               ccATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGc                     <  2:346282/149‑1 (MQ=32)
                                                                                                                                cATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGc                     >  1:353996/1‑148 (MQ=32)
                                                                                                                                        tGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                       <  3:346180/138‑1 (MQ=255)
                                                                                                                                        tGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                       <  3:340343/138‑1 (MQ=255)
                                                                                                                                        tGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                       >  4:346180/1‑138 (MQ=255)
                                                                                                                                        tGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                       >  4:340343/1‑138 (MQ=255)
                                                                                                                                        tGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg            >  4:143749/1‑149 (MQ=32)
                                                                                                                                          ccGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg            <  4:109468/147‑1 (MQ=32)
                                                                                                                                                 cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAACAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCaaa    >  8:105221/1‑148 (MQ=17)
                                                                                                                                                  ccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCaaa    >  6:119321/1‑147 (MQ=32)
                                                                                                                                                  ccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCaaa    >  5:354736/1‑147 (MQ=32)
                                                                                                                                                  ccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACa  >  7:95041/1‑149 (MQ=32)
                                                                                                                                                   |                                                                                                                                                   
GCCGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACA  >  NC_002947/5343472‑5343766

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