Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	503,835	A→T	13.2%	intergenic (+35/‑116)	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_002947	503,835	0	A	T	13.2%	141.8 / 13.7	53	intergenic (+35/‑116)	PP_0414/rpe	polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (21/25);  new base T (3/4);  total (24/29)
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.13e-01

CCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCA  >  NC_002947/503687‑503978                                                                                                                                                     |                                                                                                                                                 ccGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTaaa                                                                                                                                                 >  1:246895/1‑150 (MQ=255)     ccATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTCGCGAGCAGCTGCGGACTCACCCGCAGGCTTGCTGACGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCg                                                                                                                                             <  2:235163/150‑1 (MQ=255)       aTTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑Ga                                                                                                                                          >  1:99172/1‑150 (MQ=255)                  ggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGcg                                                                                                                               >  1:285967/1‑149 (MQ=255)                   gcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGtac                                                                                                                              <  1:671533/150‑1 (MQ=38)                     gACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAg                                                                                                                            <  2:246895/150‑1 (MQ=255)                      aCCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGCTTGCTTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAgg                                                                                                                           <  1:351694/150‑1 (MQ=255)                           gCTGATTGCGTTCTCGGTGTTGTTTCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCTGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCg                                                                                                                      <  1:116444/150‑1 (MQ=255)                                  gCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGTGCGTAGCGAGAAGCTAAAAACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGGCCAGGTTGCGTTGGTCg                                                                                                               >  1:506790/1‑150 (MQ=255)                                    gTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGac                                                                                                             <  2:285967/150‑1 (MQ=255)                                     ttCCCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGaca                                                                                                            <  2:611737/150‑1 (MQ=255)                                      tctcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                <  1:725277/145‑1 (MQ=255)                                      tctcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                >  2:725277/1‑145 (MQ=255)                                        tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGtaca                                                                                                                             <  2:596000/130‑1 (MQ=37)                                        tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGtaca                                                                                                                             >  1:596000/1‑130 (MQ=37)                                        tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAagag                                                                                                                                      <  2:119841/121‑1 (MQ=255)                                        tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAagag                                                                                                                                      >  1:119841/1‑121 (MQ=255)                                         cGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCaa                                                                                                        <  2:709125/150‑1 (MQ=255)                                          ggTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                       <  1:336403/150‑1 (MQ=255)                                          ggTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                       >  2:336403/1‑150 (MQ=255)                                           gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                      <  2:516796/150‑1 (MQ=255)                                             gttgttGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTAAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGa                                                                                                    <  2:398498/150‑1 (MQ=255)                                                gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                 <  2:99172/150‑1 (MQ=255)                                                 ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAcc                                                                                                >  1:396195/1‑150 (MQ=255)                                                      gctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCa                                                                                           >  1:585551/1‑150 (MQ=255)                                                      gctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCa                                                                                           >  2:282752/1‑150 (MQ=255)                                                            ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTc                                                                                     >  2:299883/1‑150 (MQ=255)                                                              tGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCAt                                                                                   <  1:299883/150‑1 (MQ=255)                                                                     gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                                 >  2:491371/1‑113 (MQ=255)                                                                     gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                                 <  1:491371/113‑1 (MQ=255)                                                                      ttGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTcc                                                                           <  2:585551/150‑1 (MQ=255)                                                                          gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                       <  1:185139/150‑1 (MQ=255)                                                                          gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                       <  2:396195/150‑1 (MQ=255)                                                                                           cTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                                           >  2:659359/1‑145 (MQ=255)                                                                                           cTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATgcgc                                                           <  1:659359/145‑1 (MQ=255)                                                                                              aGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCCATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGt                                                   <  1:282752/150‑1 (MQ=255)                                                                                                gACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGa                                                 >  1:132538/1‑150 (MQ=255)                                                                                                  cTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGGGTTGGt                                                                                                                 >  1:170260/1‑84 (MQ=255)                                                                                                  cTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGt                                                                                                                 <  2:170260/84‑1 (MQ=255)                                                                                                   tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAtt                                              >  2:416781/1‑150 (MQ=255)                                                                                                           cGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACa                                      >  2:472885/1‑150 (MQ=255)                                                                                                                     atgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCAt                            <  1:472885/150‑1 (MQ=255)                                                                                                                       acgcaaccTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                                           >  1:273255/7‑133 (MQ=21)                                                                                                                       acgcaaccTGTGCCGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                                           <  2:273255/127‑1 (MQ=21)                                                                                                                        gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGtt                                                                                                                         <  1:175764/54‑1 (MQ=255)                                                                                                                        gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGtt                                                                                                                         >  2:175764/1‑54 (MQ=255)                                                                                                                              cTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTa                   >  1:211518/1‑150 (MQ=255)                                                                                                                                    cGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGCTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAt             >  1:581559/1‑150 (MQ=255)                                                                                                                                      gCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                               <  2:423113/114‑1 (MQ=37)                                                                                                                                      gCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                               >  1:423113/1‑114 (MQ=37)                                                                                                                                           ttCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                      <  1:686375/134‑1 (MQ=255)                                                                                                                                           ttCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                      >  2:686375/1‑134 (MQ=255)                                                                                                                                               cgGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCa  >  2:47280/1‑150 (MQ=255)                                                                                                                                               cgGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCa  <  2:132538/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCA  >  NC_002947/503687‑503978

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