Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,602,068	T→A	11.1%	intergenic (+81/‑86)	PP_1404 → / → PP_5469	hypothetical protein/SAM‑dependent methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,602,068	0	T	A	11.1%	117.7 / 5.3	45	intergenic (+81/‑86)	PP_1404/PP_5469	hypothetical protein/SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base A (3/2);  total (23/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 = 9.78e-01

CCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTTCTCGA  >  NC_002947/1601924‑1602209                                                                                                                                                 |                                                                                                                                              ccAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCg                                                                                                                                           >  7:113711/1‑149 (MQ=255)  cAGTACGTGCCGCCAACCCATGCACGCAAGGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACCGGCGACAGGGTGTCTGGCACCGGCTTTGCCgg                                                                                                                                          >  5:25257/1‑149 (MQ=255)     cACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGACACCGGCAAAGCCGgtgc                                                                                                                                       <  6:204542/149‑2 (MQ=255)         tGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                   <  6:52196/149‑1 (MQ=255)         tGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTcgc                                                                                                                                   >  3:119977/1‑149 (MQ=255)            cgccAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTCGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCggg                                                                                                                                <  5:328396/149‑1 (MQ=255)               cAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCggg                                                                                                                                >  2:12628/1‑146 (MQ=255)               cAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCggg                                                                                                                                <  1:12628/146‑1 (MQ=255)                   ccATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGccc                                                                                                                         >  3:307943/1‑149 (MQ=255)                             aaCGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAg                                                                                                               <  3:249555/149‑1 (MQ=255)                                      cccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTg                                                                                                       <  6:321895/148‑1 (MQ=255)                                       ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTggg                                                                                                     >  6:275465/1‑149 (MQ=255)                                       ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTggg                                                                                                     >  6:270289/1‑149 (MQ=255)                                       ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCAAAGCCGgtgc                                                                                                                                       <  8:108908/115‑2 (MQ=255)                                       ccAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCAAAGCCGgtgc                                                                                                                                       >  7:108908/1‑114 (MQ=255)                                        cAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                                                                                    <  6:181005/149‑1 (MQ=255)                                         aaCCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTgggg                                                                                                    >  2:227611/1‑148 (MQ=255)                                           ccGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCAc                                                                                                 <  4:307943/149‑1 (MQ=255)                                                 ggTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTa                                                                                             <  4:119977/147‑1 (MQ=255)                                                     aTTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAg                                                                                       <  6:205536/149‑1 (MQ=255)                                                         tcGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGt                                                                                   >  8:292316/1‑149 (MQ=255)                                                             cTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCg                                                                               >  8:159487/1‑149 (MQ=255)                                                                aaGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCtt                                                                            <  8:246836/149‑1 (MQ=255)                                                                aaGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTACCGCtt                                                                            >  8:167193/1‑149 (MQ=255)                                                                    aGCGTCCAGCCTTGCTGGATTTCCTGTGGACACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTccc                                                                         >  1:27022/1‑148 (MQ=255)                                                                      cGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAg                                                                      <  7:349892/149‑1 (MQ=255)                                                                       gTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAgg                                                                     <  8:216487/149‑1 (MQ=255)                                                                        tCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGc                                                                    >  6:14380/1‑149 (MQ=255)                                                                           aGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCGTTGTAAACAGGGGACAGGGTGTCTGGAACCGGCTTTGCCGGTGTTCGCGGGAAAGCCCGATCCCAAAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGcccc                                                                 <  6:167465/149‑1 (MQ=38)                                                                             ccTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCt                                                               >  6:66341/1‑149 (MQ=255)                                                                               ttGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCttt                                                             >  2:177209/1‑149 (MQ=255)                                                                                  cTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGGCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCttttt                                                           <  3:212999/148‑1 (MQ=255)                                                                                  cTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCttttt                                                           <  7:159487/148‑1 (MQ=255)                                                                                     gATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATg                                                       >  4:80162/1‑149 (MQ=255)                                                                                     gATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATg                                                       >  1:26721/1‑149 (MQ=255)                                                                                             gtgGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAg                                               >  4:186388/1‑149 (MQ=255)                                                                                                ggCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTAc                                            <  8:113711/149‑1 (MQ=255)                                                                                                    caGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCt                                        <  6:311636/149‑1 (MQ=255)                                                                                                                  aaaCAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGcc                          >  1:23152/1‑149 (MQ=255)                                                                                                                           aCAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATgcgcg                 <  1:227611/149‑1 (MQ=255)                                                                                                                                     gaacACCGGCAAAGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGttt       >  4:261559/4‑149 (MQ=255)                                                                                                                                      aacACCGGCAAAGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTGCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGACGGCTGAGCGCGCTCAGTTtc      >  8:121138/3‑149 (MQ=38)                                                                                                                                       gCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTtc      <  3:186388/148‑1 (MQ=255)                                                                                                                                        cACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATACGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTtct     <  6:25257/148‑1 (MQ=255)                                                                                                                                          ccGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTTCTCGa  >  4:162343/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                              CCAGCACGTGCCGCCAACCCATGCACGCAACGAGTTGCCCAACCGCTTGGTGATTCTCGACTGAAGGAGCGTCCAGCCTTGCTGGATTTCCTGTGGGCACAGGCTTGCCTTGGAAACAGGCGACAGGGTGTCTGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGTGCTTGGGGCACGTTAGAATAGGCGTTCCGCTTCCCAAGGCCCCCTTTTTCATGTCCGCCAGTACTTCTGCCTCTTGCGCGCCGGATGCGCGCGCTCAGTTTCTCGA  >  NC_002947/1601924‑1602209

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