Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,851,115 A→C 8.9% A67A (GCA→GCC)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,1150AC 8.9% 128.3 / 2.7 43A67A (GCA→GCC) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (16/23);  new base C (2/2);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCG  >  NC_002947/1850969‑1851263
                                                                                                                                                  |                                                                                                                                                    
gAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCcaggcag                                                                                                                                                     <  5:328787/148‑1 (MQ=255)
  aGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGcc                                                                                                                                                  <  7:59034/149‑1 (MQ=255)
        gTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCCTGCCTGGaa                                                                                                                                             <  6:285033/148‑1 (MQ=255)
         ttGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTTCGAGTTTGCTCTTGAGGTTTAGAAAAAGCGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAAt                                                                                                                                           <  1:260906/149‑1 (MQ=255)
           ggACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATc                                                                                                                                          <  3:41564/148‑1 (MQ=255)
            gACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGc                                                                                                                                        >  8:119970/1‑149 (MQ=255)
              ccATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCg                                                                                                                                      <  5:146202/149‑1 (MQ=255)
                     ttAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCt                                                                                                                               <  7:119970/148‑1 (MQ=255)
                        aGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCCTGCCTGGaagctg                                                                                                                                         >  8:53598/1‑132 (MQ=255)
                        aGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCCTGCCTGGaagctg                                                                                                                                         <  7:53598/136‑5 (MQ=255)
                         gCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCt                                                                                                                           >  2:124527/1‑149 (MQ=255)
                                     cACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTc                                                                                                               <  1:124527/149‑1 (MQ=255)
                                     cACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTc                                                                                                               <  3:288580/149‑1 (MQ=255)
                                           aTCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAg                                                                                                         <  8:6630/149‑1 (MQ=255)
                                              gcgcAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGa                                                                                                       <  8:187073/148‑1 (MQ=255)
                                              gcgcAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGa                                                                                                       <  8:124515/148‑1 (MQ=255)
                                               cgcAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACt                                                                                                     >  7:10322/1‑149 (MQ=255)
                                               cgcAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACt                                                                                                     >  4:312498/1‑149 (MQ=255)
                                                 cAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTcc                                                                                                   >  3:65578/1‑149 (MQ=255)
                                                    cgcTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCtg                                                                                                >  4:308408/1‑149 (MQ=255)
                                                      cTGAAAGTAGCCTGCGAGTTTGGTCTTGAGGTAGTGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCtggt                                                                                              <  4:65578/149‑1 (MQ=255)
                                                      cTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCtggt                                                                                              >  6:160679/1‑149 (MQ=255)
                                                         aaaGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCtggtg                                                                                             <  2:88631/147‑1 (MQ=255)
                                                              agCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCgg                                                                                      <  3:312498/149‑1 (MQ=255)
                                                                 cTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCa                                                                                   <  1:50734/149‑1 (MQ=255)
                                                                    cGAGTTTGCTCTTCAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGCATGGTACGTCCAGCTTCCACTCCGGCCTGGAAATCGCCGAAATTCTGGATGACCTCAAGATCGCCCAGGACTCCCAGCTGGCTGAGGTCAta                                                                                 >  3:185487/1‑147 (MQ=255)
                                                                     gAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGCCCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTa                                                                               >  7:258122/1‑149 (MQ=255)
                                                                              ctTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGt                                                                      >  4:180653/1‑149 (MQ=255)
                                                                                          aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCg                                                                                                              <  1:35344/97‑1 (MQ=255)
                                                                                          aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCg                                                                                                              >  2:35344/1‑97 (MQ=255)
                                                                                          aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGcgcg                                                                 <  6:144204/142‑1 (MQ=255)
                                                                                          aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGcgcg                                                                 >  5:144204/1‑142 (MQ=255)
                                                                                          aaaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGTGAGGGc                                                            >  6:35323/1‑147 (MQ=255)
                                                                                           aaaaGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAg                                                         <  3:180653/149‑1 (MQ=255)
                                                                                                ggCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCAc                                                    <  1:76082/149‑1 (MQ=255)
                                                                                                                aTTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTcagc                                    >  1:308874/1‑149 (MQ=255)
                                                                                                                  tCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTcagc                                    >  4:127535/1‑147 (MQ=255)
                                                                                                                          ggATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTg                          <  5:160679/149‑1 (MQ=255)
                                                                                                                                tACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCgg                    <  3:308408/149‑1 (MQ=255)
                                                                                                                                  cGTCCAGCTGCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGTCCGGATGGTCGGCCAGGTCTGCCTTGTGGCTGCTGTCAGCTACCGCTCGGTGCGCGTGGGCATGGGCACCCTCGCCGAGGTCAGCCAGAGGTTTGGCCCggtg                  <  4:185487/149‑1 (MQ=255)
                                                                                                                                      cAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAg                                          >  7:276031/1‑121 (MQ=255)
                                                                                                                                      cAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAg                                          <  8:276031/121‑1 (MQ=255)
                                                                                                                                       gatttCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAg                                          <  6:266926/117‑1 (MQ=255)
                                                                                                                                       gatttCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAg                                          >  5:266926/4‑120 (MQ=255)
                                                                                                                                                  aggcCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCACCCAGAGGTTTGGCCCGGTGGTGTCCAAGCTGATCg  >  8:147159/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAAGCATGGTTGGACCATGTGGTAAGCGTCGATTCGGCACTGGATCGCGCAGCGCTGAAAGAAGCCTGCGAGTTTGCTCTTGAGGTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCG  >  NC_002947/1850969‑1851263

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