Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,853,207 G→A 7.2% intergenic (+52/‑72) relA → / → mazG ATP:GTP 3'‑pyrophosphotransferase/modified nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,853,2070GA 7.2% 162.2 / 2.0 55intergenic (+52/‑72)relA/mazGATP:GTP 3'‑pyrophosphotransferase/modified nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

ATCGAGATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACC  >  NC_002947/1853060‑1853354
                                                                                                                                                    |                                                                                                                                                       
aTCGAGATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGt                                                                                                                                                        >  3:304184/1‑149 (MQ=255)
     gaTTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGa                                                                                                                                                    >  5:337965/1‑148 (MQ=255)
      aTTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGttt                                                                                                                                                      <  5:162463/145‑1 (MQ=255)
      aTTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGttt                                                                                                                                                      >  6:162463/1‑145 (MQ=255)
         ccAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAAt                                                                                                                                                 <  1:80108/148‑1 (MQ=255)
         ccAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCt                                                                                                                                               >  7:301621/1‑149 (MQ=255)
             gCCTGGACGCGCTGGGGCGCCTGCAGGGGCGGATCTCGCAGTTGCCGAACATCATCCAGACGCGGCGTAATCGTACACCTTGAGCCCGGGGCGTATGCTTAGCGGGTAAACCCGCTGCTACAGGATGGAAGCAAGTTTGAAATCTGCAc                                                                                                                                           >  4:18242/1‑149 (MQ=255)
                 ggACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTa                                                                                                                                       >  1:143222/1‑149 (MQ=255)
                      cgcTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCtgt                                                                                                                                  <  6:89701/149‑1 (MQ=255)
                        cTGGGGCGCCTGCTGGGGCGGATCTCGGAGTTGCGGAACATCATCGAGACGCGGCGTAATCGTAGCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAAGCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTgg                                                                                                                                <  7:4104/149‑1 (MQ=255)
                           gggCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAg                                                                                                                             >  5:179223/1‑149 (MQ=255)
                            ggCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGc                                                                                                                            <  5:264345/149‑1 (MQ=255)
                               gcCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCggg                                                                                                                         <  6:337965/149‑1 (MQ=255)
                                       gggCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCAcc                                                                                                                                          <  6:24226/124‑1 (MQ=255)
                                       gggCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCAcc                                                                                                                                          >  5:24226/1‑124 (MQ=255)
                                           ggATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCcgcg                                                                                                               >  6:104314/1‑147 (MQ=255)
                                           ggATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTACTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCcgcg                                                                                                               >  8:202862/1‑147 (MQ=255)
                                                tcGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGc                                                                                                                                             <  1:329032/112‑1 (MQ=255)
                                                  gCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGcca                                                                                                      >  4:110157/1‑149 (MQ=255)
                                                  gCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTCTGGGAGCGGGCAAGCCCGCGAAGAAGcc                                                                                                       <  1:200636/148‑1 (MQ=255)
                                                          cGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGTAAGCCCGCGAAGAAGCCACCACCg                                                                                               >  2:257564/1‑147 (MQ=255)
                                                                          gcggcgTAATCGTAACCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGa                                                                           <  5:104314/149‑1 (MQ=255)
                                                                              cgTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGTAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTa                                                                       <  1:257564/149‑1 (MQ=255)
                                                                              cgTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTa                                                                       <  6:179223/149‑1 (MQ=255)
                                                                                  aTCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACAcc                                                                   <  4:304184/149‑1 (MQ=255)
                                                                                        ccccTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCTCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTgac                                                             <  1:262985/149‑1 (MQ=255)
                                                                                           cTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTgacgac                                                         >  3:152518/1‑149 (MQ=255)
                                                                                              gagaCCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGAcc                                                        <  8:301621/147‑1 (MQ=255)
                                                                                               agaCCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACctgc                                                     >  2:290557/1‑149 (MQ=255)
                                                                                                 aCCGGGGCGTCTGCTTAGCGGGTTAACCCGCACCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACctgctg                                                   >  2:351394/1‑149 (MQ=255)
                                                                                                   cGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGGATGACCTACACCCTTGACGACCTGCTGCa                                                 >  4:144585/1‑149 (MQ=255)
                                                                                                     gggCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGAcctgc                                                                                    <  4:312051/114‑1 (MQ=255)
                                                                                                     gggCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGAcctgc                                                                                    >  3:312051/1‑114 (MQ=255)
                                                                                                             gCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGcc                                       >  1:124441/1‑149 (MQ=255)
                                                                                                               cTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGccc                                      <  3:110157/149‑1 (MQ=255)
                                                                                                                     gggTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTgc                                <  7:202862/149‑1 (MQ=255)
                                                                                                                          aaCCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGAc                           >  2:268616/1‑149 (MQ=255)
                                                                                                                           aCCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGAcc                          >  1:175994/1‑149 (MQ=255)
                                                                                                                             ccGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACctgc                                                     >  5:363537/1‑120 (MQ=255)
                                                                                                                             ccGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACctgc                                                     <  6:363537/120‑1 (MQ=255)
                                                                                                                              cGCTCCCACAGGTACGGTGCAGATTTCAAACCTGCATCCATCCTGTAGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCg                        >  5:7749/1‑148 (MQ=21)
                                                                                                                              cGCTCCCACAGGTACGGTGCAGATTTCAAACCTGCATCCATCCTGTAGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCg                        <  6:7749/148‑1 (MQ=21)
                                                                                                                                      cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                <  2:124441/148‑1 (MQ=255)
                                                                                                                                      cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                <  1:198402/148‑1 (MQ=255)
                                                                                                                                      cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACg                 <  1:140188/147‑1 (MQ=255)
                                                                                                                                          atggatgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTgcgc                              <  2:322032/130‑1 (MQ=255)
                                                                                                                                          atggatgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTgcgc                              >  1:322032/1‑130 (MQ=255)
                                                                                                                                             ggtgCAGATTTCAAACCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGAc                           >  2:269937/3‑130 (MQ=255)
                                                                                                                                             ggtgCAGATTTCAAACCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGAc                           <  1:269937/128‑1 (MQ=255)
                                                                                                                                              atgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGATCCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGt        <  1:268616/149‑1 (MQ=255)
                                                                                                                                              atgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTg       >  1:48873/1‑149 (MQ=255)
                                                                                                                                                gCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTgg      >  1:1020/1‑148 (MQ=255)
                                                                                                                                                 cAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGa    <  1:351394/149‑1 (MQ=255)
                                                                                                                                                 acggTTTTAAATATGCACCGTACCTGTGGTTGGGGGCAAGCCCGCGAATAATCCACCATCCCAATCCCCGGGCCACCTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGt        <  1:290557/146‑1 (MQ=255)
                                                                                                                                                  aGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                >  6:94663/1‑136 (MQ=255)
                                                                                                                                                  aGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                <  5:94663/136‑1 (MQ=255)
                                                                                                                                                   ggTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGAcc  <  2:175994/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
ATCGAGATTCCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCACCGAAGGACCTGCCTGTATGACCTACACCCTTGACGACCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACC  >  NC_002947/1853060‑1853354

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: