Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,816 A→G 8.9% intergenic (+16/‑135) PP_0414 → / → rpe polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,8160AG 8.9% 146.4 / 3.0 54intergenic (+16/‑135)PP_0414/rpepolyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base G (2/3);  total (24/30)
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

CGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTA  >  NC_002947/503672‑503961
                                                                                                                                                |                                                                                                                                                  
cgcgAGAACCTGAGCCCGTTCCTTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGc                                                                                                                                               <  2:722122/150‑1 (MQ=255)
cgcgAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGACTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGc                                                                                                                                               <  1:475389/150‑1 (MQ=255)
         cTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCCCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTcgcg                                                                                                                                      >  1:611737/1‑150 (MQ=255)
              cccGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTcgc                                                                                                                                       >  2:350418/1‑144 (MQ=255)
              cccGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTcgc                                                                                                                                       <  1:350418/144‑1 (MQ=255)
               ccGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTaaa                                                                                                                                >  1:246895/1‑150 (MQ=255)
                   ccATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTCGCGAGCAGCTGCGGACTCACCCGCAGGCTTGCTGACGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCg                                                                                                                            <  2:235163/150‑1 (MQ=255)
                     aTTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGa                                                                                                                         >  1:99172/1‑150 (MQ=255)
                                ggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGcg                                                                                                              >  1:285967/1‑149 (MQ=255)
                                 gcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGtac                                                                                                             <  1:671533/150‑1 (MQ=38)
                                   gACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAg                                                                                                           <  2:246895/150‑1 (MQ=255)
                                    aCCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGCTTGCTTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAgg                                                                                                          <  1:351694/150‑1 (MQ=255)
                                         gCTGATTGCGTTCTCGGTGTTGTTTCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCTGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCg                                                                                                     <  1:116444/150‑1 (MQ=255)
                                                gCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGTGCGTAGCGAGAAGCTAAAAACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGGCCAGGTTGCGTTGGTCg                                                                                              >  1:506790/1‑150 (MQ=255)
                                                  gTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGac                                                                                            <  2:285967/150‑1 (MQ=255)
                                                   ttCCCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGaca                                                                                           <  2:611737/150‑1 (MQ=255)
                                                    tctcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                               <  1:725277/145‑1 (MQ=255)
                                                    tctcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                               >  2:725277/1‑145 (MQ=255)
                                                      tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGtaca                                                                                                            <  2:596000/130‑1 (MQ=37)
                                                      tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGtaca                                                                                                            >  1:596000/1‑130 (MQ=37)
                                                      tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAagag                                                                                                                     <  2:119841/121‑1 (MQ=255)
                                                      tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAagag                                                                                                                     >  1:119841/1‑121 (MQ=255)
                                                       cGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCaa                                                                                       <  2:709125/150‑1 (MQ=255)
                                                        ggTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                      <  1:336403/150‑1 (MQ=255)
                                                        ggTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                      >  2:336403/1‑150 (MQ=255)
                                                         gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                     <  2:516796/150‑1 (MQ=255)
                                                           gttgttGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTAAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGa                                                                                   <  2:398498/150‑1 (MQ=255)
                                                              gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                <  2:99172/150‑1 (MQ=255)
                                                               ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAcc                                                                               >  1:396195/1‑150 (MQ=255)
                                                                   tgctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCg                                                                                                                                                 <  1:551982/81‑1 (MQ=255)
                                                                   tgctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCg                                                                                                                                                 >  2:551982/1‑81 (MQ=255)
                                                                    gctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCa                                                                          >  1:585551/1‑150 (MQ=255)
                                                                    gctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCa                                                                          >  2:282752/1‑150 (MQ=255)
                                                                          ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTc                                                                    >  2:299883/1‑150 (MQ=255)
                                                                            tGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCAt                                                                  <  1:299883/150‑1 (MQ=255)
                                                                                   gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                <  1:491371/113‑1 (MQ=255)
                                                                                   gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                >  2:491371/1‑113 (MQ=255)
                                                                                    ttGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTcc                                                          <  2:585551/150‑1 (MQ=255)
                                                                                        gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                      <  1:185139/150‑1 (MQ=255)
                                                                                        gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                      <  2:396195/150‑1 (MQ=255)
                                                                                                         cTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                          <  1:659359/145‑1 (MQ=255)
                                                                                                         cTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                          >  2:659359/1‑145 (MQ=255)
                                                                                                            aGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCCATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGt                                  <  1:282752/150‑1 (MQ=255)
                                                                                                              gACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGa                                >  1:132538/1‑150 (MQ=255)
                                                                                                                cTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGGGTTGGt                                                                                                >  1:170260/1‑84 (MQ=255)
                                                                                                                cTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                <  2:170260/84‑1 (MQ=255)
                                                                                                                 tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAtt                             >  2:416781/1‑150 (MQ=255)
                                                                                                                         cGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACa                     >  2:472885/1‑150 (MQ=255)
                                                                                                                                   atgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCAt           <  1:472885/150‑1 (MQ=255)
                                                                                                                                     acgcaaccTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                          >  1:273255/7‑133 (MQ=21)
                                                                                                                                     acgcaaccTGTGCCGGCCTCTTCCGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                          <  2:273255/127‑1 (MQ=21)
                                                                                                                                      gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGtt                                                                                                        <  1:175764/54‑1 (MQ=255)
                                                                                                                                      gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGtt                                                                                                        >  2:175764/1‑54 (MQ=255)
                                                                                                                                            cTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTa  >  1:211518/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCGGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTA  >  NC_002947/503672‑503961

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