Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	788,586	T→G	10.5%	intergenic (+36/‑177)	ettA → / → gdhA	ADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	788,586	0	T	G	10.5%	123.1 / 8.0	48	intergenic (+36/‑177)	ettA/gdhA	ADP/ATP ratio sensor/translation inhibitor/NADP‑dependent glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (19/24);  new base G (2/3);  total (21/27)
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.33e-01

TGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCT  >  NC_002947/788441‑788729                                                                                                                                                  |                                                                                                                                                 tggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGc                                                                                                                                               <  4:279629/149‑1 (MQ=255) tggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGc                                                                                                                                               <  1:259273/149‑1 (MQ=255) tggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGc                                                                                                                                               <  7:132457/149‑1 (MQ=255)  ggtgTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGcc                                                                                                                                              >  8:109710/1‑149 (MQ=255)     gtTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCa                                                                                                                                           >  4:65654/1‑149 (MQ=255)        cttcGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttt                                                                                                                                        >  5:70733/1‑149 (MQ=255)           cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCAttttt                                                                                                                                      <  8:43527/148‑1 (MQ=255)           cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGAGGCACAAAAATGGGGCCACTTGGGGCCCCAttctt                                                                                                                                      >  8:162978/1‑148 (MQ=39)           cGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGAGGCACAAAAATGGGGCCACTTGGGGCCCCAttctt                                                                                                                                      <  7:162978/148‑1 (MQ=39)            gAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAAAATGGGGCCACTTGGGGCCCCATTCTtg                                                                                                                                     <  6:90254/148‑1 (MQ=39)                 tAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTg                                                                                                                               <  2:175767/149‑1 (MQ=255)                   aCTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTgg                                                                                                                              <  2:266227/148‑1 (MQ=255)                     tACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCg                                                                                                                            >  1:153856/1‑148 (MQ=255)                        acCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGAcc                                                                                                                        <  2:20509/149‑1 (MQ=255)                             gTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCtat                                                                                                                   <  6:70733/149‑1 (MQ=255)                             gTACGAAGCCGACCGCAAGAAGCGTCTGGGGCATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGTACCGCtat                                                                                                                   <  8:245044/149‑1 (MQ=255)                                    gCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATggc                                                                                                            <  7:109710/149‑1 (MQ=255)                                               gaagCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGgt                                                                                                >  4:272350/1‑149 (MQ=255)                                                         ggTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGta                                                                                      <  2:28196/149‑1 (MQ=255)                                                           tGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGtata                                                                                    <  1:14865/149‑1 (MQ=255)                                                             aTGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATaca                                                                                  <  4:291529/149‑1 (MQ=255)                                                               gctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGta                                                                                      >  8:9449/1‑143 (MQ=255)                                                               gctgctgcCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGta                                                                                      <  7:9449/143‑1 (MQ=255)                                                                        cAGGCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCACAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGTCGAGGCT‑AAGGTGTATTTTGTATACACTTACATaaaa                                                                       <  2:192755/149‑1 (MQ=255)                                                                             gCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCAc                                                                  >  2:141559/1‑149 (MQ=255)                                                                                ccGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCaa                                                               <  6:161778/149‑1 (MQ=255)                                                                                  gtgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGta                                                                                      <  4:252844/124‑1 (MQ=255)                                                                                  gtgtACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGta                                                                                      >  3:252844/1‑124 (MQ=255)                                                                                     tACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCAGACCGATATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAaataat                                                          >  6:98919/1‑149 (MQ=255)                                                                                        ggCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAAttt                                                        >  6:295593/1‑148 (MQ=255)                                                                                         gCACAAGTAAGTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATATATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCa                                                      <  2:102282/149‑1 (MQ=255)                                                                                             aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCa                                                      >  8:12022/1‑145 (MQ=255)                                                                                             aaGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCa                                                      <  7:12022/145‑1 (MQ=255)                                                                                              aGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCaaaaa                                                  >  5:63057/1‑148 (MQ=255)                                                                                                  aCTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGa                                             <  5:98919/149‑1 (MQ=255)                                                                                                    tgacccAGTAAGCCGGTTTACTCATGCACAAGAATGGGGCCCCAAGTTCCCCCATTTTTGAGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTAACATAAAACGCACCAAATAATAACAAAAAAACTAc                                            >  6:253391/4‑148 (MQ=255)                                                                                                           ttattccGGTTTTCTTATGAAAAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCT‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTAAAAAAAACGACATTTCGcc                                    <  2:153856/144‑1 (MQ=255)                                                                                                                      ttCTGATGCACAAGAATGGTGCCCCAATTGGAACCATTTTTGTGCCTGGCGTACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGt                          <  3:65654/148‑1 (MQ=255)                                                                                                                               acaAGAATGGGGCCCCAAGTGGCCCCATTTTTATGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGctt                >  1:193739/1‑149 (MQ=255)                                                                                                                               acaAAAATGGGGCCACTTGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCg                            <  4:130654/137‑1 (MQ=39)                                                                                                                               acaAAAATGGGGCCACTTGGGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCg                            >  3:130654/1‑137 (MQ=39)                                                                                                                                caAGAATGGGGCCCCAAGTGGCCCCATTTTTGTTCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGcttc               >  6:171309/1‑149 (MQ=255)                                                                                                                                 aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACTACATTTCGCCCTATTCCGGTGCGATTGcttct              >  6:62646/1‑149 (MQ=255)                                                                                                                                 aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCg                       >  6:266069/1‑140 (MQ=255)                                                                                                                                 aaGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCg                       <  5:266069/140‑1 (MQ=255)                                                                                                                                    aaTGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCGTGAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTg            >  1:88074/1‑149 (MQ=255)                                                                                                                                        gggCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACa       >  7:226676/1‑149 (MQ=255)                                                                                                                                             ccAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCt  <  3:272350/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TGGTGTTCTTCGAAGGTAACTACACCGAGTACGAAGCCGACCGCAAGAAGCGTCTGGGTGATGCTGCTGCCCAGCCGCACCGTGTACGGCACAAGAAACTGGCCCAGTAAGCCGGTTTTCTGATGCACAAGAATGGGGCCCCAAGTGGCCCCATTTTTGTGCCTGGCGGACCGCTATAAATGGCGAGGCG‑AAGGTGTATTTTGTATACACTTACATAAAACGCACCAAATAATTTCAAAAAAACGACATTTCGCCCTATTCCGGTGCGATTGCTTCTTGGTACATTCCT  >  NC_002947/788441‑788729

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