Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,430,392	A→G	9.7%	intergenic (‑74/+87)	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,392	0	A	G	9.7%	139.0 / 3.5	51	intergenic (‑74/+87)	PP_4769/PP_4770	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/25);  new base G (3/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATGGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGG  >  NC_002947/5430244‑5430540                                                                                                                                                     |                                                                                                                                                      gAAATGGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACCCCCCCTCCcaca                                                                                                                                                       >  3:779/1‑149 (MQ=255)           ccGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑cgcgc                                                                                                                                            <  1:15833/149‑1 (MQ=255)           ccGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCCATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑cgcgc                                                                                                                                            <  1:219700/149‑1 (MQ=255)              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)                                   tagAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑cgcgcg                                                                                                                                           >  3:109845/1‑126 (MQ=255)                                   tagAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑cgcgcg                                                                                                                                           <  4:109845/126‑1 (MQ=255)                                             gCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                          >  3:32824/1‑149 (MQ=255)                                              cacaTGGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCGGCTCCCACAGGTAT‑CGCGCGTACACCGGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                            <  2:62909/146‑1 (MQ=255)                                              cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑c                                                                                                                                                <  3:4377/110‑1 (MQ=255)                                              cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑c                                                                                                                                                >  4:4377/1‑110 (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)                                                                        tCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAgg                                                                                                                                                     <  8:83136/80‑1 (MQ=255)                                                                        tCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAgg                                                                                                                                                     >  7:83136/1‑80 (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:200836/1‑148 (MQ=255)                                                                                                cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGt                                                        >  2:4456/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                                                                                   >  6:54595/1‑98 (MQ=255)                                                                                                                       gcAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGc                                                                                   <  5:54595/98‑1 (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)                                                                                                                                                     |                                                                                                                                                      GAAATGGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGG  >  NC_002947/5430244‑5430540

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