Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,046,442 A→G 7.1% K134E (AAG→GAG)  PP_4450 → phosphoglycerate mutase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,046,4420AG 7.1% 161.5 / 2.7 56K134E (AAG→GAG) PP_4450phosphoglycerate mutase family protein
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base G (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CTCGCACCCGCGCGTTGCGTCGATACCACAGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGATACGAGATTG  >  NC_002947/5046295‑5046583
                                                                                                                                                   |                                                                                                                                             
cTCGCACCCGCGCGTTGCGTCGATACCACAGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGaa                                                                                                                                              <  3:180439/149‑1 (MQ=255)
   gCACCCGCGCGTTGCGTCGATACCACAGTCGCTCAACGACGGGATTGGGTAGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAag                                                                                                                                             <  3:276932/147‑1 (MQ=255)
     aCCCGCGCGTTGCGTCGATACCACAGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAagcagc                                                                                                                                         <  2:240029/149‑1 (MQ=255)
                        aCCACAGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCg                                                                                                                       >  8:157456/1‑148 (MQ=255)
                        aCCACAGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGCGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTATGTATATCGAGCACAGGCTATCCTGAAGCAGCTTGACGGGCACTCGACCg                                                                                                                       >  5:17912/1‑148 (MQ=255)
                          cacaGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCgg                                                                                                                          >  4:248496/1‑143 (MQ=255)
                          cacaGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCgg                                                                                                                          <  3:248496/143‑1 (MQ=255)
                          cacaGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCa                                                                                                                     >  6:25359/1‑148 (MQ=255)
                            caGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCaaa                                                                                                                   >  7:168904/1‑148 (MQ=255)
                                  cTCAATGACGGGATTGGGTGGTTGCGTATTGGGCGAGAGAAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATAGTgtc                                                                                                            <  3:228351/149‑1 (MQ=255)
                                   tCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTgtcg                                                                                                           >  3:238208/1‑149 (MQ=255)
                                        gacgGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCttt                                                                                                      >  6:20657/1‑149 (MQ=255)
                                        gacgGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCtt                                                                                                       <  1:205637/148‑1 (MQ=255)
                                        gacgGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCtt                                                                                                       >  2:205637/1‑148 (MQ=255)
                                          cgGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTc                                                                                                    <  5:25359/149‑1 (MQ=255)
                                          cgGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTc                                                                                                    <  7:157456/149‑1 (MQ=255)
                                          cgGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTc                                                                                                    <  4:238208/149‑1 (MQ=255)
                                             gATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCa                                                                                                 >  3:163965/1‑149 (MQ=255)
                                               ttGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACg                                                                                               >  8:225930/1‑149 (MQ=255)
                                               ttGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACg                                                                                               <  3:35427/149‑1 (MQ=255)
                                                        tGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGtttt                                                                                      >  4:242113/1‑149 (MQ=255)
                                                          cGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGttttt                                                                                     <  7:225930/148‑1 (MQ=255)
                                                          cGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGttttt                                                                                     <  1:47380/148‑1 (MQ=255)
                                                          cGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGttttt                                                                                     <  5:20657/148‑1 (MQ=255)
                                                          cGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTg                                                                                    <  3:173852/149‑1 (MQ=255)
                                                              ttGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATg                                                                                <  3:242113/149‑1 (MQ=255)
                                                                 ggCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCg                                                                             >  2:310921/1‑149 (MQ=255)
                                                                   cGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTc                                                                           <  4:163965/149‑1 (MQ=255)
                                                                         cagagcCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTgg                                                                     <  8:168904/149‑1 (MQ=255)
                                                                             gcCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTAc                                                                 <  6:139772/149‑1 (MQ=255)
                                                                              ccTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACt                                                                >  2:38444/1‑149 (MQ=255)
                                                                                    ttttGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAgcg                                                          <  1:310921/149‑1 (MQ=255)
                                                                                     tttGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTgtc                                                                                                            >  6:31672/1‑98 (MQ=255)
                                                                                     tttGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTgtc                                                                                                            <  5:31672/98‑1 (MQ=255)
                                                                                              gCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCag                                                <  3:290976/149‑1 (MQ=255)
                                                                                                        tGAGTCATTCCTTGATTTTGTATATCGAGCATAGGCGTTCCTGGAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAAt                                                                                 >  7:271738/1‑106 (MQ=255)
                                                                                                        tGAGTCATTCCTTGATTTTGTATATCGAGCATAGGCGTTCCTGGAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAAt                                                                                 <  8:271738/106‑1 (MQ=255)
                                                                                                        tGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGGAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAAt                                                                                 <  4:99808/106‑1 (MQ=255)
                                                                                                        tGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGGAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAAt                                                                                 >  3:99808/1‑106 (MQ=255)
                                                                                                               ttCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCAGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGc                               >  2:281185/1‑149 (MQ=255)
                                                                                                                 ccTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTa                             >  8:77776/1‑149 (MQ=255)
                                                                                                                  cTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGcc                                                                              <  8:82651/99‑1 (MQ=255)
                                                                                                                  cTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGcc                                                                              >  7:82651/1‑99 (MQ=255)
                                                                                                                  cTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTCCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTAt                            >  7:94345/1‑149 (MQ=255)
                                                                                                                               atCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACt                                                                                                               <  8:221760/53‑1 (MQ=255)
                                                                                                                                tCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTTGACAGGCACGCTATGGCTGATTGGCGga              <  7:77776/149‑1 (MQ=255)
                                                                                                                                   aGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGga              <  4:274706/146‑1 (MQ=255)
                                                                                                                                   aGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGga              >  3:274706/1‑146 (MQ=255)
                                                                                                                                     cacaGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGTATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGATAc         >  7:214604/1‑149 (MQ=255)
                                                                                                                                        aGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTg                                          <  7:217366/113‑1 (MQ=255)
                                                                                                                                        aGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTg                                          >  8:217366/1‑113 (MQ=255)
                                                                                                                                          gCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTAGAATGCCGTCGCCTGGCTACTCGAg                                                            >  6:308261/1‑93 (MQ=255)
                                                                                                                                          gCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTAGAATGCCGTCGCCTGGCTACTCGAg                                                            <  5:308261/93‑1 (MQ=255)
                                                                                                                                            gTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGATACgaga     <  4:94972/146‑1 (MQ=255)
                                                                                                                                            gTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGATACgaga     >  3:94972/1‑146 (MQ=255)
                                                                                                                                            gTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGATACGAGAtt   >  1:8380/1‑148 (MQ=255)
                                                                                                                                             ttCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGATACGAGATTg  >  1:199759/1‑148 (MQ=255)
                                                                                                                                                   aaGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGAt           <  6:231422/133‑1 (MQ=255)
                                                                                                                                                   aaGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGAt           >  5:231422/1‑133 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CTCGCACCCGCGCGTTGCGTCGATACCACAGTCGCTCAACGACGGGATTGGGTGGATGCGTATTGGGCGAGAGCAGAGCCTGCATTTTGTGACGGCGAAGGTGCTGAGTCATTCCTTGATTTTGTATATCGAGCACAGGCGTTCCTGAAGCAGCTTGCCGGGCACTCGGCCGCAAATACTGTCGTCTTTTCGCACGGACAGTTTTTGAATGCCGTCGCCTGGCTACTCGAGCGCAGACCGCAGAGCCTGGACAGGCACGCTATGGCTGATTGGCGGAGATACGAGATTG  >  NC_002947/5046295‑5046583

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: ‑