Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,688,066	C→G	12.5%	P348A (CCG→GCG)	PP_4148 →	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,688,066	0	C	G	12.5%	80.3 / 3.8	32	P348A (CCG→GCG)	PP_4148	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (15/13);  new base G (2/2);  total (17/15)
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.23e-01

GCCTGATGAGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGAA  >  NC_002947/4687941‑4688213                                                                                                                              |                                                                                                                                                    gcctgATGAGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCGTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCgg                                                                                                                              <  1:104284/149‑1 (MQ=255) gcctgATGAGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCgg                                                                                                                              <  5:22341/149‑1 (MQ=255)     gatgaGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAg                                                                                                                           >  4:178278/1‑148 (MQ=255)      atgaGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCa                                                                                                                         <  3:232654/149‑1 (MQ=255)              cGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACAt                                                                                                                  <  3:178278/148‑1 (MQ=255)                     gCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTg                                                                                                           >  8:141590/1‑148 (MQ=255)                          cAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTccccc                                                                                                     >  6:322495/1‑149 (MQ=255)                             ccGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTc                                                                                                  <  7:141590/149‑1 (MQ=255)                                             gtCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTc                                                                                   <  5:322495/148‑1 (MQ=255)                                                    gCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGACAGCCGGGAGCACTGACATGGCCTTGTCCACCCTCAACCGCCGGCGCTTCGAGCGtt                                                                            >  3:92126/1‑148 (MQ=255)                                                    gCACCCGGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGttt                                                                           <  7:245329/149‑1 (MQ=255)                                                                 cTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCACCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAAacgc                                                              >  4:225139/1‑149 (MQ=255)                                                                                  tggtggTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCa                                                                         <  1:161083/121‑1 (MQ=255)                                                                                  tggtggTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCa                                                                         >  2:161083/1‑121 (MQ=255)                                                                                          aaaCTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCtgg                                                     >  8:205065/1‑133 (MQ=255)                                                                                          aaaCTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCtgg                                                     <  7:205065/133‑1 (MQ=255)                                                                                            aCTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTc                                   >  6:177396/1‑149 (MQ=255)                                                                                                  cGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGAt                              >  4:36226/1‑148 (MQ=255)                                                                                                    gCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGAAATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCt                           >  6:162362/1‑149 (MQ=255)                                                                                                      gacctgacctACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGt                         >  1:11040/1‑149 (MQ=255)                                                                                                        cctgacctACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCAGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTACGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTc                       >  7:70863/1‑149 (MQ=255)                                                                                                         ctgacctACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTTGTCGCTGTGGCTGTTCCTGATTCTGTTCg                      <  5:162362/149‑1 (MQ=255)                                                                                                                    cTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCtgt                                           >  4:77565/1‑117 (MQ=255)                                                                                                                    cTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCtgt                                           <  3:77565/117‑1 (MQ=255)                                                                                                                    cTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTg           >  4:141349/1‑149 (MQ=255)                                                                                                                    cTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTg           >  4:136357/1‑149 (MQ=255)                                                                                                                      aagTCGATGCGCGGATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAAc                                                                 >  7:28700/3‑93 (MQ=255)                                                                                                                      aagTCGATGCGCGGATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAAc                                                                 <  8:28700/91‑1 (MQ=255)                                                                                                                       agTCGATGCGCGGATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCtggtg                                                   <  2:247476/105‑1 (MQ=255)                                                                                                                       agTCGATGCGCGGATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCtggtg                                                   >  1:247476/2‑106 (MQ=255)                                                                                                                             tCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGaa  <  3:141349/149‑1 (MQ=255)                                                                                                                             tCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGaa  <  3:136357/149‑1 (MQ=255)                                                                                                                              |                                                                                                                                                    GCCTGATGAGTTTCGAAGCGGCCATCAACCGCGACTACCCGGTGGTCTTCGGCACCCTGTTCATCTTCACCCTGCTAGGGCTGGTGGTGAAACTGATCGGCGACCTGACCTACACCTTGGTCGATCCGCGCATCGACTTTGCCAGCCGGGAGCACTGACATGGCCTTGTCCCCCCTCAACCGCCGGCGCTTCGAGCGTTTCAAGGCCAACCGCCGTGGCTGGTGGTCGCTGTGGCTGTTCCTGATTCTGTTCGGGCTGAGCCTGGGCGCCGAA  >  NC_002947/4687941‑4688213

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