Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,667,259	A→C	9.1%	intergenic (+37/+34)	nhaP → / ← PP_4975	NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,667,259	0	A	C	9.1%	118.9 / 4.6	44	intergenic (+37/+34)	nhaP/PP_4975	NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base A (21/19);  new base C (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

GAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTT  >  NC_002947/5667112‑5667397                                                                                                                                                    |                                                                                                                                           gAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCATGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCCGGAGaa                                                                                                                                           >  3:347306/1‑149 (MQ=255)    gAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGAAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCATTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAAccc                                                                                                                                        >  5:113496/1‑149 (MQ=255)      aCGGGACCTGCTGCTGTCGATCACCTACATCGTCGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCcgc                                                                                                                                      <  5:250447/149‑1 (MQ=255)         ggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCGGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTTCACCCGCGaa                                                                                                                                   <  7:326556/149‑1 (MQ=255)              tgctgctgTGGTTGATCTTCTTGGTGGTGTTGTGGTCGATTGTGGTGCAGGGGTTGAGCTTTGGGGGGGTGGTTCGCTAGGTCAGCGACCAGCCTTGAGATTGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGc                                                                                                                              <  4:347306/149‑1 (MQ=25)                ctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTTCACCCGCGaatgggc                                                                                                                              <  8:106339/147‑6 (MQ=255)                ctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGcc                                                                                                                             <  5:105200/148‑1 (MQ=255)                   ctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGca                                                                                                                         <  2:297200/149‑1 (MQ=255)                      tCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGcac                                                                                                                        <  2:78413/147‑1 (MQ=255)                      tCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAg                                                                                                                      <  2:185564/149‑1 (MQ=255)                        gATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCGGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAgg                                                                                                                     <  4:150498/148‑1 (MQ=255)                                      tggtgTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTAtt                                                                                                      <  5:2006/149‑1 (MQ=255)                                           ttgtCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCTAAGAGTCCGGCACAGGCAACCAGGGTCTATTCGACt                                                                                                 >  1:129973/1‑149 (MQ=255)                                                gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCg                                                                                             >  6:218179/1‑148 (MQ=255)                                                gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCg                                                                                             >  1:333946/1‑148 (MQ=255)                                                  cgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGAGCCCAGCCTTGAGAATGCCGTGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTCTGCGCATGCCGaa                                                                                           >  8:288590/1‑148 (MQ=255)                                                        tGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGcc                                                                                              >  6:361780/1‑139 (MQ=255)                                                        tGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGcc                                                                                              <  5:361780/139‑1 (MQ=255)                                                          gggCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCggg                                                                                     <  6:113496/144‑1 (MQ=255)                                                             cAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGACCCAGCGAAGAGGCCGGCACTGGCACCCAGGGTATATTCGCCTGCCGAATCCGGGAACTgg                                                                               >  8:256703/1‑149 (MQ=255)                                                              accgcttGACCATTGGGCGGCTGCTGCCCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGCGTCAACCCGCCAAGACGCCGGCACACGCAACCAGCGTCTATACGACTCCCGAATcc                                                                                        >  2:85353/6‑139 (MQ=255)                                                              aGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATcc                                                                                        <  1:85353/139‑1 (MQ=255)                                                                   ttGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACt                                                                                                 <  7:229468/125‑1 (MQ=255)                                                                   ttGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACt                                                                                                 >  8:229468/1‑125 (MQ=255)                                                                     gAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                             <  7:315106/111‑1 (MQ=255)                                                                     gAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                             >  8:315106/1‑111 (MQ=255)                                                                       gCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACg                                                                      <  2:333946/148‑1 (MQ=255)                                                                        cATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACAAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGt                                                                     >  5:24106/1‑148 (MQ=255)                                                                          ttGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                      <  1:344041/129‑1 (MQ=255)                                                                          ttGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                      >  2:344041/1‑129 (MQ=255)                                                                                     gtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTg                                                       >  5:205860/1‑149 (MQ=255)                                                                                                     cccAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCg                                       <  5:218179/149‑1 (MQ=255)                                                                                                               agaATGCCGGGGCCGCTTTTCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCAAAGGAAACAAGGGTCTATTAGAATGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATc                             >  6:20976/1‑149 (MQ=255)                                                                                                                 cgcgaatgggCCGCAAAGCGGCCCATTCGCGGGTGACCCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTCATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATcc                            >  3:181002/8‑148 (MQ=255)                                                                                                                     ccGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGcc                       >  5:306976/1‑149 (MQ=255)                                                                                                                     aatgggCCGCAAAGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGcc                       >  4:17833/4‑149 (MQ=255)                                                                                                                              cTTTGCGGCACATTCGCGGGAGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATTGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAc              <  7:288590/149‑1 (MQ=255)                                                                                                                                 tGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACg             <  6:306976/147‑1 (MQ=255)                                                                                                                                   cGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAaca         >  3:345793/1‑149 (MQ=255)                                                                                                                                    ggCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAacac        >  7:125751/1‑149 (MQ=255)                                                                                                                                       ctcttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCa      >  6:299354/4‑148 (MQ=255)                                                                                                                                        cATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAaca         <  2:280459/144‑1 (MQ=255)                                                                                                                                        cATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAaca         >  1:280459/1‑144 (MQ=255)                                                                                                                                          ttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGc                                                                                               >  6:349686/1‑56 (MQ=18)                                                                                                                                          ttCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCtt  >  4:106918/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                           GAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTT  >  NC_002947/5667112‑5667397

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