Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,044,933	A→G	12.9%	S159G (AGC→GGC)	pstC →	phosphate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,044,933	0	A	G	12.9%	73.9 / 2.9	31	S159G (AGC→GGC)	pstC	phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (2/2);  total (18/13)
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.90e-01

CCAGGATTTCCTTCCGGCCAAGCTTGCCGCGGTGTTTCGCACCATCGTCGAGTTGCTCGCCGCCATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGATGCGCTGACCAGCGTACCTATGAAAAC  >  NC_002947/3044791‑3045075                                                                                                                                               |                                                                                                                                               ccAGGATTTCCTTCCGGCCAAGCTTGCCGCGGTGTTTCGCACCATCGTCGAGTTGCTCGCCGCCATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAgc                                                                                                                                          >  5:224895/1‑149 (MQ=255)                    aaGCTTGCCGCGGTGTTTCGCACCATCGTCGAGTTGCTCGCCGCCATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGttttt                                                                                                                       <  6:224895/149‑1 (MQ=255)                                      cGCACCATCGTCGAGTTGCTCGCCGCCATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTg                                                                                                      <  3:35119/148‑1 (MQ=255)                                                     ttGCTCGCCGCCATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGctgcc                                                                                      >  1:173148/1‑149 (MQ=255)                                                         tCGCCGCCATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTgcg                                                                                  >  2:266258/1‑149 (MQ=255)                                                         tCGCCGCCATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCAAATCGTACGCAGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATAACGTTTTTCGGAACCACCCTGAGATGCCCAGGCCTGCTGCCTgcg                                                                                  >  1:76920/1‑149 (MQ=255)                                                               ccATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAg                                                                                                        <  6:93617/121‑1 (MQ=255)                                                               ccATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAg                                                                                                        >  5:93617/1‑121 (MQ=255)                                                                 aTCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGt                                                                          >  1:248543/1‑149 (MQ=255)                                                                     ccAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGTCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTg                                                                      >  1:142355/1‑149 (MQ=255)                                                                            ggtgTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTgg                                                                     <  4:99744/143‑1 (MQ=255)                                                                            ggtgTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTgg                                                                     >  3:99744/1‑143 (MQ=255)                                                                               gtACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTgg                                                                     <  1:45354/140‑1 (MQ=255)                                                                               gtACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTgg                                                                     >  2:45354/1‑140 (MQ=255)                                                                                  cGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTgcgg                                                                                 >  6:154124/1‑125 (MQ=255)                                                                                  cGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTgcgg                                                                                 <  5:154124/125‑1 (MQ=255)                                                                                                           aTCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTcc                                  >  5:53005/1‑147 (MQ=255)                                                                                                             cccGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGGTTCGCAg                               >  1:93541/1‑148 (MQ=255)                                                                                                              ccGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGa                             <  1:266258/149‑1 (MQ=255)                                                                                                                 gCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGa                             <  4:19842/146‑1 (MQ=255)                                                                                                                 gCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGa                             >  3:19842/1‑146 (MQ=255)                                                                                                                        gTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGATGCGCTGAc                    >  3:91786/1‑148 (MQ=255)                                                                                                                         tCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCgcggcg                                        >  7:219312/1‑127 (MQ=255)                                                                                                                         tCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCgcggcg                                        <  8:219312/127‑1 (MQ=255)                                                                                                                                 ccGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTACCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGATGCGCTGACCAGCGTACCt          <  8:102673/149‑1 (MQ=255)                                                                                                                                  cGCCTGGCTGAACGGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGATGCGCTGACCAGCGTACCTa         >  2:236844/1‑149 (MQ=255)                                                                                                                                  cGCCTGGCTGAACGGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGATGCGCTGACCAGCGTACCTa         >  2:232886/1‑149 (MQ=255)                                                                                                                                         cTGAACGGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGATGCGCTGACCAGCGTACCTATGAAAAc  <  1:236844/149‑1 (MQ=255)                                                                                                                                         cTGAACGGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGATGCGCTGACCAGCGTACCTATGAAAAc  <  1:232886/149‑1 (MQ=255)                                                                                                                                         cTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGccc                                                                                                  >  5:99210/1‑53 (MQ=255)                                                                                                                                         cTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGccc                                                                                                  <  6:99210/53‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               CCAGGATTTCCTTCCGGCCAAGCTTGCCGCGGTGTTTCGCACCATCGTCGAGTTGCTCGCCGCCATCCCCAGCGTGGTGTACGGCCTGTGGGGCATCTATGTGGTCATCCCGGCGATTCGTCCGCTGACCGCCTGGCTGAACAGCGAGCTGGGCTGGATCCCGTTTTTCGGCACCTCCCTGAGTGGCCCAGGCCTGCTGCCTGCGGCGCTGGTGCTGGCGATCATGATCCTGCCGACCATCGCGGCGGTTTCCCAGGATGCGCTGACCAGCGTACCTATGAAAAC  >  NC_002947/3044791‑3045075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑