Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,430,404	Δ1 bp	10.2%	intergenic (‑86/+75)	PP_4769 ← / ← PP_4770	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,430,404	0	T	.	10.2%	141.4 / 11.7	49	intergenic (‑86/+75)	PP_4769/PP_4770	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/23);  new base . (3/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

TGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCC  >  NC_002947/5430257‑5430551                                                                                                                                                     |                                                                                                                                                    tGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTa                                                                                                                                                    <  3:4198/149‑1 (MQ=255)      cTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCTCTCCCACAGGTAT‑CGCGCGTACAccc                                                                                                                                               <  5:115876/149‑1 (MQ=255)              cGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGc                                                                                                                                       <  2:80471/149‑1 (MQ=255)              cGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGc                                                                                                                                       <  2:77073/149‑1 (MQ=255)                gggtagtaGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCgg                                                                                                                                     <  6:205490/147‑1 (MQ=255)                                gCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                                     >  3:32824/1‑149 (MQ=255)                                 cacaTGGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCGGCTCCCACAGGTAT‑CGCGCGTACACCGGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                                       <  2:62909/146‑1 (MQ=255)                                 cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                                       <  2:58523/146‑1 (MQ=255)                                                cgcgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCg                                                                                                     >  1:99546/1‑149 (MQ=255)                                                 gcgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCgg                                                                                                    <  2:99546/149‑1 (MQ=255)                                                   gATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCCCAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAg                                                                                                             >  4:261504/1‑138 (MQ=255)                                                   gATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAg                                                                                                             <  3:261504/138‑1 (MQ=255)                                                       tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                                                <  5:126831/147‑1 (MQ=255)                                                       tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                                                >  6:126831/1‑147 (MQ=255)                                                         tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGaaa                                                                                                 <  3:131941/144‑1 (MQ=255)                                                                 cgagctcTCCTTTTGGGCCTGGAGGGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCCCAGGTCT‑CGCGCGGACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATg                                                                                    <  4:103367/143‑1 (MQ=38)                                                                  ggAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGc                                                                                   <  4:32824/149‑1 (MQ=255)                                                                       cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGCCACCCCCGCAACCCCACGTAT‑AGCCCGGACACCCCTGGGCCCCGGCGGTCCCGCCACCGGGCGCAAACCGCCCCCCGCAATGCCCGCAGAATTcccgcc                                                                              >  6:19223/1‑145 (MQ=255)                                                                        tCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAcc                                                                             >  2:170708/1‑149 (MQ=255)                                                                        tCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                              >  8:113565/1‑148 (MQ=255)                                                                                 gCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACg                                                                    >  7:245446/1‑149 (MQ=255)                                                                                  ccTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACg                                                                    >  7:46635/1‑148 (MQ=255)                                                                                  ccTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCCCGCCCCCAGGTAT‑AGCGCGTACACCCGTGGGAGCGGGGGTGCCCGCCCCGGGGCGCAAAGCGCCCCCGGCAATGCCAGCAGTATGCCa                                                                                 >  8:158295/1‑135 (MQ=255)                                                                                   cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCGGTATGCCAGACCTGTCAGACGt                                                                   >  5:135267/1‑148 (MQ=255)                                                                                   cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGt                                                                   >  2:4456/1‑148 (MQ=255)                                                                                   cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGt                                                                   >  2:200836/1‑148 (MQ=255)                                                                                         tGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtc                                                            <  6:135267/149‑1 (MQ=255)                                                                                           ccGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcg                                                           <  1:200836/148‑1 (MQ=255)                                                                                            cGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                         <  1:170708/149‑1 (MQ=255)                                                                                            cGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                         <  1:4456/149‑1 (MQ=255)                                                                                                          gcAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGc                                                                                              <  5:54595/98‑1 (MQ=255)                                                                                                          gcAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGc                                                                                              >  6:54595/1‑98 (MQ=255)                                                                                                                cTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGt                                     >  2:20308/1‑149 (MQ=255)                                                                                                                 ttCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGt                                     >  7:242355/1‑148 (MQ=255)                                                                                                                   cgcgGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTc                                    <  4:779/147‑1 (MQ=255)                                                                                                                        acacGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTg                             <  8:46635/149‑1 (MQ=255)                                                                                                                         cacGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGc                            >  1:31265/1‑149 (MQ=255)                                                                                                                            gCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAGCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCa                          >  5:231435/1‑148 (MQ=255)                                                                                                                                gCTCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                                     >  7:109348/1‑117 (MQ=21)                                                                                                                                gCTCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                                     <  8:109348/117‑1 (MQ=21)                                                                                                                                  tCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAaccac                   >  2:112533/1‑149 (MQ=37)                                                                                                                                  tCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCt                                                                            <  7:138819/92‑1 (MQ=18)                                                                                                                                  tCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCt                                                                            >  8:138819/1‑92 (MQ=18)                                                                                                                                      acaGGTAT‑CGCGCGTACACCCGGGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGGCATCGCTGCCAAAACCACCGTc               <  5:19223/149‑1 (MQ=255)                                                                                                                                        aGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCgg             >  2:191325/1‑149 (MQ=255)                                                                                                                                        aGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCgg             <  8:245446/149‑1 (MQ=255)                                                                                                                                                    gTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                                 >  3:255809/1‑86 (MQ=255)                                                                                                                                                    gTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                                 <  4:255809/86‑1 (MQ=255)                                                                                                                                                    gTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGccc  <  8:242355/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCC  >  NC_002947/5430257‑5430551

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