Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 839,209 C→T 8.3% L164L (CTG→CTA)  prs ← ribose‑phosphate pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947839,2090CT 8.3% 133.0 / 2.8 48L164L (CTG→CTA) prsribose‑phosphate pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base T (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GCGCCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGCATG  >  NC_002947/839062‑839356
                                                                                                                                                   |                                                                                                                                                   
gcgcCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCCCGGCGTGCTGGAGAGGGCGTGCGCGCACGACACCACCGATGTCCGGGGAGTCGATCATCAg                                                                                                                                                   <  2:371008/150‑1 (MQ=255)
gcgcCCTTCCACGTCGCCTATGATGTGCATAACCTCTGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAg                                                                                                                                                   <  2:751927/150‑1 (MQ=255)
 cgcCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGGCGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAgg                                                                                                                                                  >  1:412673/1‑150 (MQ=255)
  gcCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGt                                                                                                                                                 >  1:424610/1‑150 (MQ=255)
   cccTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGtt                                                                                                                                                >  1:138570/1‑150 (MQ=255)
      ttCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTtctc                                                                                                                                             <  2:38134/150‑1 (MQ=255)
      ttCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTtctc                                                                                                                                             <  2:150650/150‑1 (MQ=255)
                 cGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGt                                                                                                                                  >  2:601554/1‑150 (MQ=255)
                        gtgCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGAt                                                                                                                           <  2:410416/150‑1 (MQ=255)
                          gCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGc                                                                                                                                      <  2:708574/137‑1 (MQ=255)
                          gCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGc                                                                                                                                      >  1:708574/1‑137 (MQ=255)
                           cATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATgtc                                                                                                                        >  1:482060/1‑150 (MQ=255)
                                          gTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACggg                                                                                                         >  1:99369/1‑150 (MQ=255)
                                           ttAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGc                                                                                                        >  2:565343/1‑150 (MQ=255)
                                             aGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGa                                                                                                      <  2:138570/150‑1 (MQ=255)
                                                     cgcgGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGATAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGCTCCTCAGGTTCTCGACACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGa                                                                                              <  2:228492/150‑1 (MQ=255)
                                                                 cGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACgg                                                                                  <  1:628234/150‑1 (MQ=255)
                                                                  gatgatACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGGCTTCGATGTCGTCAACCAGGACTGGCTATCCGTAGTTGTTGTCGACggg                                                                                 >  2:720415/1‑150 (MQ=255)
                                                                          ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACCCCACCAATTTCGGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCga                                                                         >  2:292974/1‑150 (MQ=255)
                                                                          ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTtctc                                                                                                                                             <  1:555240/82‑1 (MQ=255)
                                                                          ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTtctc                                                                                                                                             >  2:555240/1‑82 (MQ=255)
                                                                          ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCga                                                                         >  1:33698/1‑150 (MQ=255)
                                                                          ccAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATACCga                                                                         >  2:377314/1‑150 (MQ=255)
                                                                           cAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCgaa                                                                        <  2:559470/150‑1 (MQ=255)
                                                                           cAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATACCgaa                                                                        <  1:377314/150‑1 (MQ=255)
                                                                                     cccAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGa                                                              <  1:367704/150‑1 (MQ=255)
                                                                                      ccAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGaa                                                             >  1:307159/1‑150 (MQ=255)
                                                                                       cAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAAt                                                            >  2:259344/1‑150 (MQ=255)
                                                                                              ttGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCgaagaa                                                                     >  1:156061/1‑134 (MQ=255)
                                                                                              ttGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCgaagaa                                                                     <  2:156061/134‑1 (MQ=255)
                                                                                                ggCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGc                                                   >  1:423608/1‑150 (MQ=255)
                                                                                                ggCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGc                                                   <  2:482060/150‑1 (MQ=255)
                                                                                                   gACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTg                                                <  1:328659/150‑1 (MQ=255)
                                                                                                      aGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAg                                             <  1:720415/150‑1 (MQ=255)
                                                                                                        cgcgTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGt                                           <  1:250339/150‑1 (MQ=255)
                                                                                                               gcgcACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGt                                    <  1:601554/150‑1 (MQ=255)
                                                                                                                cgcACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTg                                   <  2:715201/150‑1 (MQ=255)
                                                                                                                        accaccGATGTCCGGGGAGACGATCATTAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGAtg                                                                                            >  2:331572/1‑85 (MQ=255)
                                                                                                                        accaccGATGTCCGGGGAGACGATCATTAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGAtg                                                                                            <  1:331572/85‑1 (MQ=255)
                                                                                                                        accaccGATGTCCGGGGAGACGATCATTAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGAtg                                                                                            <  1:707464/85‑1 (MQ=255)
                                                                                                                        accaccGATGTCCGGGGAGACGATCATTAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGAtg                                                                                            >  2:707464/1‑85 (MQ=255)
                                                                                                                         ccaccGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACgg                          >  2:343637/1‑150 (MQ=255)
                                                                                                                           accGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTc                        >  1:534941/1‑150 (MQ=255)
                                                                                                                           accGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTc                        <  2:307159/150‑1 (MQ=255)
                                                                                                                                    cGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGGGGTTGTCGACGGCGTTAACGAAGAATCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTAc               >  1:616964/1‑150 (MQ=255)
                                                                                                                                     ggggAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCa                                              <  1:175816/118‑1 (MQ=255)
                                                                                                                                     ggggAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCa                                              >  2:175816/1‑118 (MQ=255)
                                                                                                                                        gagaCGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAAAAGAGTTTAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACa           >  2:477128/1‑150 (MQ=255)
                                                                                                                                              atcatcAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGc     <  2:424610/150‑1 (MQ=255)
                                                                                                                                               tcatcaGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGCa    >  1:207814/1‑150 (MQ=255)
                                                                                                                                                 atcaGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGGAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGCATg  >  1:617834/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGCCCTTCGACGTCGCCGATGATGTGCATAACCTCGGAGTGGTTAGCCTTTTCGCGGCGTTTGTCGATGATACCCAGATCGACGCCCAGGGACTTGGCGACAGCGCGTGCGCGCACGACACCACCGATGTCCGGGGAGACGATCATCAGGTTCTCGAAACGCTGGTCTTCGATGTCGTCAACCAGTACGGGCGAGCCGTAGATGTTGTCGACGGGGATATCGAAGAAGCCTTGAATCTGGTCAGCGTGCAGGTCGACGGTGAGTACACGGTCGATACCTACGACAGTGAGCATG  >  NC_002947/839062‑839356

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