Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,534,919	G→A	8.5%	intergenic (+110/‑73)	secA → / → argJ	protein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,534,919	0	G	A	8.5%	129.8 / 2.9	47	intergenic (+110/‑73)	secA/argJ	protein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (21/22);  new base A (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

CAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGCTTCTGCCGGT  >  NC_002947/1534773‑1535066                                                                                                                                                   |                                                                                                                                                     cAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggc                                                                                                                                                    >  5:53106/1‑149 (MQ=255)      agTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGcc                                                                                                                                                <  6:53106/148‑1 (MQ=255)       gTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCtt                                                                                                                                              <  7:221345/149‑1 (MQ=255)        ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTtg                                                                                                                                             >  5:77243/1‑149 (MQ=255)          cAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTtgtg                                                                                                                                           >  3:207025/1‑149 (MQ=255)          cAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTtgtg                                                                                                                                           >  1:123117/1‑149 (MQ=255)                        gTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGCGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑c                                                                                                                            <  5:231275/149‑1 (MQ=255)                                     aTCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGCTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCtt                                                                                                               >  8:92811/1‑149 (MQ=255)                                          aaCCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTc                                                                                                           <  4:29492/148‑1 (MQ=255)                                               tACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTgggg                                                                                                                                                      >  1:92939/1‑101 (MQ=255)                                               tACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTgggg                                                                                                                                                      <  2:92939/101‑1 (MQ=255)                                                      tACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGACCTAGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGttt                                                                                              >  8:170087/1‑149 (MQ=255)                                                         aaGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGttttt                                                                                            >  2:144084/1‑148 (MQ=255)                                                            ccTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGAc                                                                                                        <  4:186642/133‑1 (MQ=255)                                                             cTGGGGTCGGTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                                                       <  5:169752/149‑1 (MQ=255)                                                               ggggTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTg                                                                                                                              <  1:14447/109‑1 (MQ=255)                                                               ggggTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTg                                                                                                                              >  2:14447/1‑109 (MQ=255)                                                                gggTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGa                                                                                     <  4:207025/148‑1 (MQ=255)                                                                   tCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAgcgc                                                                                 >  2:91694/1‑149 (MQ=255)                                                                     gCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGctc                                                                               <  3:191960/149‑1 (MQ=255)                                                                       tttGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGctct                                                                              >  4:100028/1‑148 (MQ=255)                                                                           cGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCGTTTTTCGACATCATCGTTTTTTATAGGAGCGAACTAATgg                                                                         >  6:228239/1‑149 (MQ=255)                                                                           cGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATg                                                                          <  1:91694/148‑1 (MQ=255)                                                                            ggCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGc                                                                        <  3:85408/149‑1 (MQ=255)                                                                            ggCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGc                                                                        >  2:34254/1‑149 (MQ=255)                                                                                 tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                                                                   >  8:28462/1‑149 (MQ=255)                                                                                 tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                                                                   >  2:208750/1‑149 (MQ=255)                                                                                    cGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGCGCCGGCCTTGTGTCGCGATGGCGTG‑CCAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                                                                >  5:65984/1‑149 (MQ=255)                                                                                            aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                                        >  8:83168/1‑149 (MQ=255)                                                                                              ggCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                                        <  6:248026/147‑1 (MQ=255)                                                                                              ggCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCttt                                                      >  4:52844/1‑149 (MQ=18)                                                                                                 cGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGcc                                                   <  1:34254/149‑1 (MQ=255)                                                                                                 cGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTTGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGcc                                                   >  4:67407/1‑149 (MQ=18)                                                                                                  gCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGcc                                                   >  8:114197/1‑148 (MQ=255)                                                                                                     cccACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAcc                                               >  6:184744/1‑149 (MQ=255)                                                                                                           ggTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAc                                         >  3:44606/1‑149 (MQ=255)                                                                                                            gTACCGCGCTCGTCCATGAACGGGTGCGATGCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAcc                                        <  6:77243/149‑1 (MQ=255)                                                                                                             tACCGCGCTCGTGCATGATCGGGGGCGATCGTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAccc                                       <  6:65984/149‑1 (MQ=255)                                                                                                                cgcgCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCATTTGCGCACCCTGCACCCGGt                                    >  6:182503/1‑149 (MQ=255)                                                                                                                       gTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCgg                               <  7:83168/147‑1 (MQ=255)                                                                                                                       gTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCgg                               <  1:208750/147‑1 (MQ=255)                                                                                                                       gTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTGGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGt                             >  8:99411/1‑149 (MQ=255)                                                                                                                         ccATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGttt                           <  7:92811/149‑1 (MQ=255)                                                                                                                          cATGGACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGttt                           <  3:52844/148‑1 (MQ=38)                                                                                                                              aaCGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAAc                      <  1:208306/149‑1 (MQ=255)                                                                                                                                cGGGCGCGATCCTTGTGGGGGCGGCCTTGTTTCGCGTTGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTc                    <  5:228239/149‑1 (MQ=255)                                                                                                                                 gggCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGTCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGGTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTc                    >  1:188410/1‑149 (MQ=255)                                                                                                                                   gcgcGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGt                 <  4:44606/149‑1 (MQ=255)                                                                                                                                       ggtacctgtgGGAGCGGCCTTGTGTTGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCg             <  3:67407/143‑1 (MQ=255)                                                                                                                                           cTTGTGGGGGCGACCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCTGTATCGCtt          <  7:170087/148‑1 (MQ=255)                                                                                                                                                  gggCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGCTTCTGCCGGt  <  7:28462/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTG‑CGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGCTTCTGCCGGT  >  NC_002947/1534773‑1535066

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