Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,534,950	A→G	7.8%	intergenic (+99/‑84)	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,950	0	A	G	7.8%	152.5 / 2.3	51	intergenic (+99/‑84)	secA/argJ	protein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/20);  new base G (2/2);  total (29/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.87e-01

GTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGC  >  NC_002947/1534803‑1535098                                                                                                                                                    |                                                                                                                                                     gtgtGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACGCAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATc                                                                                                                                                    >  2:807491/1‑150 (MQ=255) gggtgGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGGGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTATCAACACAAAGCCGCTCCCACAGGAACGGCGCTAGTCCATGAAGGGGCGCGacc                                                                                                                                                    >  2:1160338/3‑148 (MQ=255)     ggTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTtg                                                                                                                                                <  1:940920/150‑1 (MQ=255)            gCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGggc                                                                                                                                         >  1:1129279/1‑150 (MQ=255)                gaagTTCAAGCATTGCCACGTTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGAGCGGcc                                                                                                                                     <  1:1199164/150‑1 (MQ=18)                gaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTGAAGTATGAACGGGCGCGATCCTTGTGGGGGCGGcc                                                                                                                                     <  1:1001634/150‑1 (MQ=39)                            ttGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGAt                                                                                                                         <  2:1129279/150‑1 (MQ=255)                                cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATgggg                                                                                                                     >  2:732251/1‑150 (MQ=255)                                         aTCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCg                                                                                                            >  2:482281/1‑150 (MQ=255)                                               tGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCt                                                                                                      >  2:593835/1‑150 (MQ=255)                                               tGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGagcg                                                                                                                                        >  2:758901/1‑116 (MQ=14)                                               tGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGGACGAGCGCGGTACCTGTGGGagcg                                                                                                                                        <  1:758901/116‑1 (MQ=2)                                                  tCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCtttt                                                                                                   >  2:519344/1‑150 (MQ=255)                                                   cTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCttttt                                                                                                  <  2:551074/150‑1 (MQ=255)                                                    tGNACCTTACGTTGTACAAGCCTGGGGNCGCTTTGCGGCCNTTTCGCNANACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                                                 >  1:747697/1‑150 (MQ=255)                                                          ttACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCAGCTTCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcatc                                                                                           >  1:85941/1‑150 (MQ=255)                                                               ttGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGa                                                                                                                                                      <  1:199258/86‑1 (MQ=255)                                                               ttGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGa                                                                                                                                                      >  2:199258/1‑86 (MQ=255)                                                                  tACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATg                                                                                                                        >  2:1087286/1‑113 (MQ=255)                                                                  tACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATg                                                                                                                        <  1:1087286/113‑1 (MQ=255)                                                                  tACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                                                   >  1:189235/1‑150 (MQ=255)                                                                   aCAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGttttt                                                                                  >  2:796697/1‑150 (MQ=255)                                                                       gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGcc                                                                                                                                     >  1:1117586/1‑95 (MQ=255)                                                                       gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGcc                                                                                                                                     <  2:1117586/95‑1 (MQ=255)                                                                          tGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGa                                                                           <  1:796697/150‑1 (MQ=255)                                                                               tCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGTACCCGCCCCTCTTTCGACATGCTCGCGTTTTGCAGgcgcgca                                                                      >  1:1199898/1‑149 (MQ=255)                                                                                             tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                                                         <  1:1032252/149‑1 (MQ=255)                                                                                             tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                                                         >  2:1032252/1‑149 (MQ=255)                                                                                                 gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtctt                                                    >  1:1186730/1‑150 (MQ=255)                                                                                                   aaCACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTTTTTTTTTTTGTTTTTTCTAGGAGCGCTCTAATGTCTGttggtctttg                                                  >  2:732701/1‑148 (MQ=255)                                                                                                     aacaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttggtc                                                <  2:85941/149‑1 (MQ=255)                                                                                                       caAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                                              <  1:164817/150‑1 (MQ=255)                                                                                                          ggCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTg                                           <  2:35884/150‑1 (MQ=255)                                                                                                                  ccgcaGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCt                                   <  2:196660/147‑1 (MQ=255)                                                                                                                   cacaGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTg                                  >  1:217801/1‑150 (MQ=255)                                                                                                                    acaGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGc                                 >  1:731847/1‑150 (MQ=255)                                                                                                                    acaGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTGTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGATGTTGGTATTAGTCCTTTGCCCACCTTGc                                 >  2:561170/1‑150 (MQ=255)                                                                                                                     caGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCa                                <  2:731847/150‑1 (MQ=255)                                                                                                                       ggTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAcc                              <  2:217801/150‑1 (MQ=255)                                                                                                                          atcgcgCCCGTTCATGGACGAGCGCGGTACCTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGGGCTCTAATGGCTGCTGGGCTGGGTCCGTTGCCCACCCTGCACCCgg                           >  2:1184241/3‑150 (MQ=25)                                                                                                                             gcgcTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGctc                                                                     >  1:962119/1‑105 (MQ=255)                                                                                                                             gcgcTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGctc                                                                     <  2:962119/105‑1 (MQ=255)                                                                                                                                   gTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGt                                                           <  1:354927/109‑1 (MQ=255)                                                                                                                                   gTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGt                                                           >  2:354927/1‑109 (MQ=255)                                                                                                                                     ccATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGtttt                <  1:519344/150‑1 (MQ=255)                                                                                                                                      cATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGa                                                                                               >  1:1153555/1‑70 (MQ=255)                                                                                                                                      cATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGa                                                                                               <  2:1153555/70‑1 (MQ=255)                                                                                                                                      cATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCt                                                             >  2:592624/1‑104 (MQ=255)                                                                                                                                      cATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCt                                                             <  1:592624/104‑1 (MQ=255)                                                                                                                                      cATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTg               >  2:1172594/1‑150 (MQ=255)                                                                                                                                       aTGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGa              >  2:373241/1‑150 (MQ=255)                                                                                                                                                 gcgATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATc    >  1:353915/1‑150 (MQ=255)                                                                                                                                                  cgATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCg   <  1:373241/150‑1 (MQ=255)                                                                                                                                                   gATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGc  >  2:304016/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGAACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCCGGTTCCGGGTTTTGAACTCGGTATCGC  >  NC_002947/1534803‑1535098

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑