Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,154,084 G→A 8.7% intergenic (‑26/+57) PP_3657 ← / ← PP_3658 p‑nitrobenzoate reductase NfnB/aromatic compound MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,154,0840GA 8.7% 157.4 / 4.8 57intergenic (‑26/+57)PP_3657/PP_3658p‑nitrobenzoate reductase NfnB/aromatic compound MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (21/31);  new base A (2/3);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

AGGTTGACGCCTGTTGCGCTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGCGGGATTGGTTGATCATGGTCACTGCAACAGCCC  >  NC_002947/4153937‑4154229
                                                                                                                                                   |                                                                                                                                                 
aGGTTGACGCCTGTTGCGCTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCgg                                                                                                                                                  >  7:325830/1‑149 (MQ=255)
 ggTTGACGCCTGTTGCGCTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGt                                                                                                                                                 <  6:347758/149‑1 (MQ=255)
  gTTGACGCCTGTTGCGCTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGt                                                                                                                                                 >  2:129431/1‑148 (MQ=255)
      aCGCCTGTTGCGCTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTc                                                                                                                                            <  2:225180/149‑1 (MQ=255)
        gCCTGTTGCGCTAAAGCGAGCAACGTTCAGCGTCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTcgc                                                                                                                                          >  6:317823/1‑149 (MQ=255)
        gCCTGTTGCGCTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGACCTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGCCAATTGCCCAGTCGAGCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTcgc                                                                                                                                          >  3:261429/1‑149 (MQ=255)
               gcgcTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGc                                                                                                                                    >  7:195603/1‑148 (MQ=255)
                cgcTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGa                                                                                                                                  <  1:331588/149‑1 (MQ=255)
                   tAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCaa                                                                                                                                      >  7:228615/1‑142 (MQ=255)
                   tAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCaa                                                                                                                                      <  8:228615/142‑1 (MQ=255)
                   tAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATc                                                                                                                                >  8:147842/1‑148 (MQ=255)
                      aGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGATACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAAt                                                                                                                             <  3:30379/148‑1 (MQ=255)
                       gCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATaa                                                                                                                           <  7:107422/149‑1 (MQ=255)
                        cgagataCGTTCAGCTTCGCTTCGATTTCCTCTAGCGGAAGTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGTTCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAg                                                                                                                          <  8:103203/143‑1 (MQ=255)
                         gatcAACGTTAAGCATCGCTTCGATTTCCTCTAGCTGAACTGGGGTTGGTAGGAACGACCTGATGCTCCTACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGc                                                                                                                         <  5:33823/146‑1 (MQ=255)
                         gAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGc                                                                                                                         <  6:247268/149‑1 (MQ=255)
                         gAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGc                                                                                                                         <  8:325830/149‑1 (MQ=255)
                               cGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCAGTGATTCGCGCAAGCGATCAATAAGCATACtt                                                                                                                   <  4:361021/149‑1 (MQ=255)
                                      cATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCt                                                                                                            >  8:169914/1‑149 (MQ=255)
                                        tCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGt                                                                                                          <  3:161667/149‑1 (MQ=255)
                                            ttCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCAGTGATTCGCGCAAGCGa                                                                                                                                  <  3:350957/121‑1 (MQ=255)
                                            ttCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCAGTGATTCGCGCAAGCGa                                                                                                                                  <  3:356624/121‑1 (MQ=255)
                                            ttCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCAGTGATTCGCGCAAGCGa                                                                                                                                  >  4:350957/1‑121 (MQ=255)
                                            ttCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCAGTGATTCGCGCAAGCGa                                                                                                                                  >  4:356624/1‑121 (MQ=255)
                                               gATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCa                                                                                                    <  8:195603/148‑1 (MQ=255)
                                                     ctctAGCGGATCTGGGGTTGGGAGGAACGACCTGATGCTCCGACGGGAGACAATTGCCCAGTCGATCTAAGGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTAc                                                                                             <  4:261429/149‑1 (MQ=255)
                                                     ctctAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTAc                                                                                             >  8:288951/1‑149 (MQ=255)
                                                           cGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACggtcag                                                                                       <  5:237735/149‑1 (MQ=255)
                                                                cTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGgtcaggtcag                                                                                  <  7:169914/149‑1 (MQ=255)
                                                                 tGGGGTTGGGAGGAACGACCTGATTCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATaa                                                                                                                           <  4:376697/107‑1 (MQ=255)
                                                                     gTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGt                                                                             <  7:147842/149‑1 (MQ=255)
                                                                        gggAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTgg                                                                           <  2:285020/148‑1 (MQ=255)
                                                                             gaACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCATGGGCAAcc                                                                     >  5:109864/1‑149 (MQ=255)
                                                                                   ccTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGcc                                                               >  2:122398/1‑149 (MQ=255)
                                                                                   ccTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGcc                                                               >  6:384778/1‑149 (MQ=255)
                                                                                       aTGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCcttc                                                            <  5:317823/148‑1 (MQ=255)
                                                                                          cTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCttcttc                                                         <  7:277537/148‑1 (MQ=255)
                                                                                          cTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCttcttc                                                         <  1:129431/148‑1 (MQ=255)
                                                                                           tCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCt                                                       >  3:196214/1‑149 (MQ=255)
                                                                                             cGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCa                                                     <  3:155355/149‑1 (MQ=255)
                                                                                             cGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCCGTCTCTCAGAATTACGGTCAGGTCAGc                                                                                 >  1:164840/1‑121 (MQ=255)
                                                                                             cGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCCGTCTCTCAGAATTACGGTCAGGTCAGc                                                                                 <  2:164840/121‑1 (MQ=255)
                                                                                                    agaCAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAg                                               <  7:288951/148‑1 (MQ=255)
                                                                                                        aaTTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAg                                          >  2:231502/1‑149 (MQ=255)
                                                                                                              ccAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGa                                     <  7:95284/148‑1 (MQ=255)
                                                                                                               cAGTCGATCAATCGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGc                                   >  4:109841/1‑149 (MQ=255)
                                                                                                               cAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGc                                   <  2:145283/149‑1 (MQ=255)
                                                                                                                        aatcGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGCGGGATTGGt                          <  3:109841/146‑1 (MQ=255)
                                                                                                                          aCGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACg                                                                                            <  6:271955/81‑1 (MQ=255)
                                                                                                                          aCGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACg                                                                                            >  5:271955/1‑81 (MQ=255)
                                                                                                                            gggCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGCgg                                 >  2:369248/1‑138 (MQ=255)
                                                                                                                            gggCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGCgg                                 <  1:369248/138‑1 (MQ=255)
                                                                                                                                      gTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGCGGGATTGGTTGATCATGGTCAc             <  5:384778/148‑1 (MQ=255)
                                                                                                                                       tcagAACCTACCGGAGATTCGCGCAAGCGATCAATAAGCATACTAGGCACCTGTCTCTCAGAATTACGGTCACGTCAGCTCGTGGGCAACCGTGGCCTACTTCCTCACAACATCCGAGTCCGAGCGGGATTGGTTGATCagg                  >  8:149125/4‑140 (MQ=255)
                                                                                                                                       gccGTACCTACCGTTGATTCGCGGAAGCGATCAATTAGCATACTTGGCACCTGTCGCTCAGAATTCCGGTCAGGTCAGCTCGTTGGCAACGGTGGCCTTCTTCCTGACTACCGCCGAGTCCGAGCTGGATTGGTTGATCATg                  <  7:149125/141‑1 (MQ=255)
                                                                                                                                         cGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGCGGGATTGGTTGATCATGGTCACTGc          >  2:102371/1‑148 (MQ=255)
                                                                                                                                           aaCCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGCGGGATTGGTTGATCATGGTCACTGCAAc       >  7:331010/1‑149 (MQ=255)
                                                                                                                                                aCCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGACAGCTCGTGGGCAACCGTGGCCTTCTTCCTCCCTACAGCCGAGTCAGAGCGGGATTGGTTGATCATGGTCACTGCAACAGccc  >  6:110789/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGTTGACGCCTGTTGCGCTAAAGCGAGCAACGTTCAGCATCGCTTCGATTTCCTCTAGCGGAACTGGGGTTGGGAGGAACGACCTGATGCTCCGACGTGAGACAATTGCCCAGTCGATCAAACGGGCTGTTTCGTCCGAACCTACCGGTGATTCGCGCAAGCGATCAATAAGCATACTTGGCACCTGTCTCTCAGAATTACGGTCAGGTCAGCTCGTGGGCAACCGTGGCCTTCTTCCTCACTACAGCCGAGTCCGAGCGGGATTGGTTGATCATGGTCACTGCAACAGCCC  >  NC_002947/4153937‑4154229

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