Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,718,027	G→A	10.5%	P68P (CCG→CCA)	PP_2377 →	acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,718,027	0	G	A	10.5%	97.4 / 2.1	38	P68P (CCG→CCA)	PP_2377	acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (18/16);  new base A (2/2);  total (20/18)
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.76e-01

CACCATTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGTTCCTGTCGCTGCTGTT  >  NC_002947/2717878‑2718173                                                                                                                                                      |                                                                                                                                                   caccaTTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAAATACTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTGATCTACCTGTCGGCCCGCGACAAAACCGTCAAGCcg                                                                                                                                                    >  2:465180/1‑150 (MQ=255)  accaTTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCcgc                                                                                                                                                   >  1:711590/1‑150 (MQ=255)  accaTTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTAGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCcgc                                                                                                                                                   >  1:205237/1‑150 (MQ=255)   ccaTTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCcgcg                                                                                                                                                  >  2:87373/1‑150 (MQ=255)   ccaTTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCATCAGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCcgcg                                                                                                                                                  >  1:365955/1‑150 (MQ=255)    caTTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCgcgc                                                                                                                                                 >  1:834796/1‑150 (MQ=255)               gATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGc                                                                                                                                       >  2:388475/1‑149 (MQ=255)               gATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTCGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGc                                                                                                                                       <  1:388475/149‑1 (MQ=255)                  cttcttCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCaa                                                                                                                                   <  1:87373/150‑1 (MQ=255)                   ttcttcGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAAt                                                                                                                                  >  1:860520/1‑150 (MQ=255)                      ttcGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCACGCCAGTTCCTGCTCAATCgg                                                                                                                               >  2:1335766/1‑150 (MQ=255)                      ttcGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCACGCCAGTTCCTGCTCAATCgg                                                                                                                               >  2:1335815/1‑150 (MQ=255)                            ttCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTg                                                                                                                         >  1:788035/1‑150 (MQ=255)                                      ccGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTcc                                                                                                               >  2:683663/1‑150 (MQ=255)                                            ccATCGGCCAACTGCTCACCGCCCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTTATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGCTCCTGCCt                                                                                                         >  1:127692/1‑150 (MQ=255)                                                             accgaccgTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGCCGCTTCTCATCGTCCCTGCCTACTGGTTCTACACAGcc                                                                                        >  2:691423/1‑150 (MQ=255)                                                               cgaccgTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCt                                                                                      <  2:711590/150‑1 (MQ=255)                                                                     tGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGc                                                                                <  2:579330/150‑1 (MQ=255)                                                                         gcCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCata                                                                            <  1:683663/150‑1 (MQ=255)                                                                               ggCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTa                                                                      <  2:937756/150‑1 (MQ=255)                                                                                gCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTAc                                                                     >  1:1093087/1‑150 (MQ=255)                                                                                gCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTAc                                                                     >  2:106283/1‑150 (MQ=255)                                                                                 cGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACt                                                                    <  2:1093087/150‑1 (MQ=255)                                                                                 cGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACt                                                                    >  1:1008047/1‑150 (MQ=255)                                                                                      aCATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCa                                                                              <  1:1080013/135‑1 (MQ=255)                                                                                      aCATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCa                                                                              >  2:1080013/1‑135 (MQ=255)                                                                                           ttcttcGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGc                                                          >  2:894647/1‑150 (MQ=255)                                                                                               tcGTCATCAGTGGGCCGATCATCTCCCTGTCTACCCGCGACAAAACCGTCGACCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAAt                                                      <  2:1240499/150‑1 (MQ=255)                                                                                                    atcaGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCACGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATAGTGCTACTGGCATTCAGCCAATGGATg                                                 <  1:1335766/150‑1 (MQ=255)                                                                                                    atcaGTGGGCTGGTAATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCACGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCTTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATCGCTACTGGCATTCAGCCAATGGATg                                                 <  1:1335815/150‑1 (MQ=255)                                                                                                           ggCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATc                                          >  2:1298908/1‑150 (MQ=255)                                                                                                             cTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAg                                        <  2:834796/150‑1 (MQ=255)                                                                                                              tGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGt                                       >  1:844197/1‑150 (MQ=255)                                                                                                                 tCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGt                                    <  1:1298908/150‑1 (MQ=255)                                                                                                                           gTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTggcg                          <  2:1008047/150‑1 (MQ=255)                                                                                                                                 ccGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCt                    >  2:243444/1‑150 (MQ=255)                                                                                                                                   gcgACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGt                  >  1:967814/1‑150 (MQ=255)                                                                                                                                   gcgACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGt                  <  1:894647/150‑1 (MQ=255)                                                                                                                                      aCAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACAt                                           <  2:1201997/122‑1 (MQ=255)                                                                                                                                      aCAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACAt                                           >  1:1201997/1‑122 (MQ=255)                                                                                                                                       cAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGTTCCt              <  1:691423/150‑1 (MQ=255)                                                                                                                                          aaCCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGTTCCTGTc           >  1:1089706/1‑150 (MQ=255)                                                                                                                                                  aGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTTGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGTTCCTGTCGCTGCTGt   <  1:446785/150‑1 (MQ=255)                                                                                                                                                   gCCGCGCCAGTTACTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGTTCCTGTCGCTGCTGtt  >  2:1417742/1‑150 (MQ=255)                                                                                                                                                      gcgcCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGTTc                >  2:975519/1‑133 (MQ=255)                                                                                                                                                      gcgcCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGTTc                <  1:975519/133‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CACCATTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCATCGGCCAACTGCTCACCGACCGTGGCGCCGTAGGCGTCGACATCTTCTTCGTCATCAGTGGGCTGGTCATCTACCTGTCGACCCGCGACAAAACCGTCGAGCCGCGCCAGTTCCTGCTCAATCGGGCGCTGCGCATCGTCCCTGCCTACTGGTTCTACACAGCCCTGATGGCCATATTGCTACTGGCATTCAGCCAATGGATGCCACATCAGTCGTTTGACTGGCGGCACCTGTTCCTGTCGCTGCTGTT  >  NC_002947/2717878‑2718173

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