Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	503,852	C→T	8.9%	intergenic (+52/‑99)	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,852	0	C	T	8.9%	151.7 / 5.1	56	intergenic (+52/‑99)	PP_0414/rpe	polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (31/20);  new base T (3/2);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

GACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCG  >  NC_002947/503707‑503996                                                                                                                                                   |                                                                                                                                                 ggcgtggCTGATTGGGTTCTCGGTGTTGTTGCTGCTGACCCGGGGGTGGTTGCGTGGGCGTAGCGAGAGGCGGGGGACTCCACCGCCGGCTTGATGACGGCTTGCCTGTACCGGCCCCTTCCCCGCTAAAGGCG‑GAGGAGGCCGCCACAg                                                                                                                                              <  2:235149/144‑1 (MQ=255)  aCCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTTGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAgg                                                                                                                                             >  2:427323/1‑150 (MQ=255)       gCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCg                                                                                                                                        <  1:577535/150‑1 (MQ=255)         tGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGtt                                                                                                                                      >  1:1583856/1‑150 (MQ=255)                                        ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTc                                                                                                       >  2:608571/1‑150 (MQ=255)                                        ccTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTc                                                                                                       >  2:1345560/1‑150 (MQ=255)                                                 gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGACTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GGAGAGGCCGGCACAGGTTGCGGTGGTCGACACTAATTGATACCAATCAGCCTTCATTGCtttt                                                                                              >  1:1230548/1‑149 (MQ=255)                                                    gCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCt                                                                                           <  1:1345560/150‑1 (MQ=255)                                                      gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                                     >  1:1573692/1‑122 (MQ=255)                                                      gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                                     <  2:1573692/122‑1 (MQ=255)                                                      gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                                         >  2:1080265/1‑150 (MQ=255)                                                      gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                                         >  2:1535427/1‑150 (MQ=255)                                                          gCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTa                                                                                     >  2:1357130/1‑150 (MQ=255)                                                                       cTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAc                                                                        >  1:680545/1‑150 (MQ=255)                                                                               tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAtt                                                                >  1:1609098/1‑150 (MQ=255)                                                                                cAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTc                                                                                                       >  1:661438/1‑110 (MQ=255)                                                                                cAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTc                                                                                                       <  2:661438/110‑1 (MQ=255)                                                                                cAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCAGCAAAGGTTGCTTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                                                               <  1:608571/150‑1 (MQ=255)                                                                                  aCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                                                             <  2:680545/150‑1 (MQ=255)                                                                                  aCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                                                             <  1:1535427/150‑1 (MQ=255)                                                                                     gCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCaa                                                          <  1:427323/150‑1 (MQ=255)                                                                                                  tgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                                         >  1:95980/1‑106 (MQ=255)                                                                                                  tgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                                         <  2:95980/106‑1 (MQ=255)                                                                                                  tgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgc                                                   <  2:113252/144‑1 (MQ=255)                                                                                                  tgGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgc                                                   >  1:113252/1‑144 (MQ=255)                                                                                                    gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATc                                                                                                    >  1:395437/1‑93 (MQ=255)                                                                                                    gCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATc                                                                                                    <  2:395437/93‑1 (MQ=255)                                                                                                      ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                                         <  2:336576/70‑1 (MQ=255)                                                                                                      ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                                         <  2:336666/70‑1 (MQ=255)                                                                                                      ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                                         >  1:336666/1‑70 (MQ=255)                                                                                                      ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                                         >  1:336576/1‑70 (MQ=255)                                                                                                      ttGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGc                                         >  2:174775/1‑150 (MQ=255)                                                                                                             tgccGGCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCTCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGc                                  >  2:926769/3‑150 (MQ=37)                                                                                                                cGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                                       >  1:942490/1‑126 (MQ=255)                                                                                                                cGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                                       <  2:942490/126‑1 (MQ=255)                                                                                                                  gCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTg                             >  1:617090/1‑150 (MQ=255)                                                                                                                  gCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                                       <  2:1347160/124‑1 (MQ=37)                                                                                                                  gCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                                       >  1:1347160/1‑124 (MQ=37)                                                                                                                  gCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAAt                                                                                                               >  2:634081/1‑68 (MQ=18)                                                                                                                  gCCTCTTC‑CGGCTTTAGCCGCGAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAAt                                                                                                               <  1:634081/68‑1 (MQ=18)                                                                                                                     tctTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGt                                                     >  1:1161283/1‑123 (MQ=255)                                                                                                                     tctTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGt                                                     <  2:1161283/123‑1 (MQ=255)                                                                                                                     tctTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccc                          >  2:1452978/1‑150 (MQ=255)                                                                                                                     tctTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccc                          >  1:266981/1‑150 (MQ=255)                                                                                                                      ctTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGcccc                         >  2:426109/1‑150 (MQ=255)                                                                                                                       ttCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccccc                        >  2:377805/1‑150 (MQ=255)                                                                                                                       ttCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccccc                        >  1:144634/1‑150 (MQ=255)                                                                                                                            ggCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCAt                   >  1:1557209/1‑150 (MQ=255)                                                                                                                                  aGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATtctctc             >  2:339539/1‑150 (MQ=255)                                                                                                                                          agagGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCTCTCTGCCGAt      <  1:377805/150‑1 (MQ=255)                                                                                                                                             ggCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGt                                                                                       <  1:786534/66‑1 (MQ=255)                                                                                                                                             ggCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGt                                                                                       >  2:786534/1‑66 (MQ=255)                                                                                                                                              gCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGt                                                                                       >  1:976583/1‑65 (MQ=255)                                                                                                                                              gCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGt                                                                                       <  2:976583/65‑1 (MQ=255)                                                                                                                                              gCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTcg  <  2:144634/150‑1 (MQ=255)                                                                                                                                              gCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTcg  >  2:893227/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTCTCTCTGCCGATTTCG  >  NC_002947/503707‑503996

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