Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,344,252 T→C 10.9% F498L (TTC→CTC)  pgi‑II → glucose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,344,2520TC10.9% 142.5 / 7.9 55F498L (TTC→CTC) pgi‑IIglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (19/30);  new base C (3/3);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GCGCTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCT  >  NC_002947/5344106‑5344391
                                                                                                                                                  |                                                                                                                                           
gcgcTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTc                                                                                                                                           >  7:188008/1‑149 (MQ=32)
  gcTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCg                                                                                                                                          >  8:345936/1‑148 (MQ=32)
      aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTa                                                                                                                                        >  1:201037/1‑146 (MQ=32)
      aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTa                                                                                                                                        >  8:60758/1‑146 (MQ=32)
      aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTACat                                                                                                                                     >  2:223408/1‑148 (MQ=11)
       aGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTACatg                                                                                                                                    <  1:223408/149‑3 (MQ=11)
       aGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTACatg                                                                                                                                    <  1:74181/149‑3 (MQ=11)
          ccTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCg                                                                                                                                 <  1:307702/149‑1 (MQ=32)
          ccTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCg                                                                                                                                 >  2:173778/1‑149 (MQ=32)
              aaCGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTAt                                                                                                                             >  2:167591/1‑149 (MQ=32)
                 gAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTg                                                                                                                          >  1:22971/1‑149 (MQ=32)
                   aGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTggg                                                                                                                        <  2:85360/149‑1 (MQ=32)
                   aGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTggg                                                                                                                        <  4:36311/149‑1 (MQ=32)
                   aGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTggg                                                                                                                        <  2:119616/149‑1 (MQ=32)
                   aGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTggg                                                                                                                        <  1:227723/149‑1 (MQ=32)
                             gAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTg                                                                                                                          >  6:335683/1‑137 (MQ=32)
                             gAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTg                                                                                                                          <  5:335683/137‑1 (MQ=32)
                              aaaaGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGcc                                                                                                             >  4:286506/1‑149 (MQ=32)
                                      gCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGAt                                                                                                       <  1:167591/147‑1 (MQ=32)
                                         ccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTg                                                                                                  <  2:292812/149‑1 (MQ=32)
                                         ccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTg                                                                                                  >  5:186241/1‑149 (MQ=32)
                                          cGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTgg                                                                                                 <  7:345936/149‑1 (MQ=32)
                                           gCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTggg                                                                                                <  1:173778/149‑1 (MQ=32)
                                           gCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTggg                                                                                                <  2:201037/149‑1 (MQ=32)
                                            cacaAGGTGATCCCGGGTAACCGCCCTAGCAACACGTTGGTGGTCGAGCTCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTgggg                                                                                               <  7:256211/149‑1 (MQ=17)
                                             acaAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGc                                                                                              <  2:22971/149‑1 (MQ=32)
                                                        ccGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgg                                                                                    <  7:60758/148‑1 (MQ=32)
                                                         cGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCtt                                                                                                           >  6:198504/1‑124 (MQ=32)
                                                         cGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCtt                                                                                                           <  5:198504/124‑1 (MQ=32)
                                                          gggTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggc                                                                                  >  4:820/1‑148 (MQ=32)
                                                            gTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaag                                                                               <  3:286506/149‑1 (MQ=32)
                                                             tAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaagg                                                                              <  2:272013/149‑1 (MQ=32)
                                                               aCCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaagga                                                                             <  3:284645/148‑1 (MQ=32)
                                                                   cccGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaaggagctgg                                                                        <  2:75742/149‑1 (MQ=34)
                                                                         cAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGagctgggcaaggagctgggcaagg                                                                  <  3:820/149‑1 (MQ=35)
                                                                           acacGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGc                                                                >  8:314578/1‑149 (MQ=37)
                                                                                  ggtggtCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTAc                                                          <  8:188008/148‑1 (MQ=34)
                                                                                      gtCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAgc                                                      <  7:314578/148‑1 (MQ=34)
                                                                                        cGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGCATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaaggagctgg                                                                        <  4:104314/128‑1 (MQ=17)
                                                                                        cGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGCATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaaggagctgg                                                                        >  3:104314/1‑128 (MQ=17)
                                                                                            cgcATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGt                                               <  2:371736/149‑1 (MQ=34)
                                                                                               aTCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCgg                                            >  7:238177/1‑149 (MQ=34)
                                                                                                                  gTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGcc                         <  5:295835/149‑1 (MQ=34)
                                                                                                                        tggtggCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGa                                   >  4:32337/1‑133 (MQ=34)
                                                                                                                        tggtggCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGa                                   <  3:32337/133‑1 (MQ=34)
                                                                                                                              ctctGTACGAACACCTTGTGCTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGaa                                  <  4:352922/125‑1 (MQ=11)
                                                                                                                              ctctGTACGAACACCTTGTGCTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGaa                                  >  3:352922/4‑128 (MQ=11)
                                                                                                                              ccATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTc              <  1:60078/148‑1 (MQ=34)
                                                                                                                              ccATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTc              >  2:60078/1‑148 (MQ=34)
                                                                                                                               tctGTACGAACACCTTGTGCTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCa            >  2:370488/3‑149 (MQ=11)
                                                                                                                               cATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCa            <  8:238177/149‑1 (MQ=34)
                                                                                                                                      gAGCACAAGGTGTTAGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGATCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAgg     <  7:182662/149‑1 (MQ=21)
                                                                                                                                          acaAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGc    <  4:181215/146‑1 (MQ=34)
                                                                                                                                          acaAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGc    >  3:181215/1‑146 (MQ=34)
                                                                                                                                          acaAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCt  >  2:262513/1‑148 (MQ=34)
                                                                                                                                                  |                                                                                                                                           
GCGCTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCT  >  NC_002947/5344106‑5344391

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: