Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,740,811	(T)7→6	intergenic (+48/+705)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,740,805	0	T	.	93.5%	119.6 / 2.1	31	intergenic (+42/+711)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (12/17);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740658‑4740947                                                                                                                                                    |                                                                                                                                                gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGAAATAGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                                 >  3:47214/1‑148 (MQ=255)     tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                              >  5:126973/1‑147 (MQ=255)       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                           >  1:205140/1‑147 (MQ=255)       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                           >  6:12601/1‑147 (MQ=255)        aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACCCCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                           <  8:87329/147‑2 (MQ=255)        aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                            <  1:198720/146‑1 (MQ=255)        aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCACAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                           <  4:139971/147‑2 (MQ=255)          agaaaTGCTCTCCAGCCCGTACAAGATCGGCCGGCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                        <  7:92140/148‑5 (MQ=255)                                aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCTTTGCCCTGAAAGACCGCAAGTAAGCTGAAGTAGCCGAACGAAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                                 >  6:42946/1‑148 (MQ=255)                                       gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGc                                                                                                         <  6:125622/148‑1 (MQ=255)                                              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGaa                                                                                                  <  1:156101/148‑1 (MQ=255)                                              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGaa                                                                                                  <  2:45144/148‑1 (MQ=255)                                               cgcCAGCGGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAag                                                                                                 <  8:74013/148‑1 (MQ=255)                                                tccAGCTGTACACCGGCGAGCAGAAGCGCGAAATCGTTTCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAaga                                                                                                <  8:128942/147‑1 (MQ=255)                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGG‑CAAGGCCGCGAaga                                                                                                <  5:82770/148‑1 (MQ=38)                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAaga                                                                                                <  7:14403/148‑1 (MQ=255)                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAaga                                                                                                <  2:46186/148‑1 (MQ=255)                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGC‑CTTGCCCGCGAaga                                                                                                <  3:15621/148‑1 (MQ=255)                                                   aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGc                                                                                             <  2:62519/148‑1 (MQ=255)                                                      tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGCCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCATAGCGGC‑CCTGCCCCCGACGCGGCCgg                                                                                          >  5:114338/1‑148 (MQ=255)                                                         acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGc                                                                                         >  7:41031/1‑146 (MQ=255)                                                                           gcgACATGGTGGCCCTGAACGACCCAAATTACGCTGCAGTAGCCGAACGCAATACGGGTGCCCTCGGGCAGCTTTTTTTTTTTGTCCTGG‑GCCGGCCTCTTCGCGGC‑CTTGCCCGCGAATAGGCCGGCACAGACACCACAAAACCTCa                                                                     <  8:30034/148‑1 (MQ=255)                                                                                      gCCCTGAAAGACGGCATGTAAGGTGCAGTAGCCTAACGCAAGAAGGGTGCCGTGGGGCAG‑CTTTTTTTATTTGTCTTGTGCCGGCCTCTGCGCGGC‑TTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                           <  8:21874/147‑1 (MQ=255)                                                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGC‑CTTGCCCCCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGGTTTTATccgccc                                                       >  8:84809/1‑147 (MQ=255)                                                                                                  cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCa                                              >  6:661/1‑148 (MQ=255)                                                                                                  cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCAGCCAAGAGGCCGGCCACGACACCACAAAACCTCAGGTTTTCTACGCACCcgcacgca                                              >  6:7769/1‑148 (MQ=255)                                                                                                        tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGCAGCGGCCGGCCCAGCCACCCCAAAACCTCACTTTTTATCCCCCCGCGCACCCAGCCCtt                                        >  1:152334/1‑148 (MQ=255)                                                                                                           gcgccAGTCGCCCACCGCACGACGGCGGCCCGCGGGCAGCC‑TTTTTTATTTGTCGTGGGCCGGCCGCTTCGCGGC‑CTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                     <  8:119189/144‑1 (MQ=255)                                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCGGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                     <  7:120876/148‑1 (MQ=255)                                                                                                             tccAGTAGCCGAACGAAAGACGGCTGCCCGCGGGCAGCC‑TTTTTAATTTGCCCGGGGCCGGCCTCTGCGCGGC‑CTGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                   <  5:12601/146‑1 (MQ=255)                                                                                                               cAGTAGCCGCACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CGTGCCCGCGACGCGGCCGGCCCAGACCCCACAAACCCGCAGGTTTTCTCCGCCCGCTCCCGCAGCCCTTTGAgggg                                 >  7:9562/1‑147 (MQ=255)                                                                                                               cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CATGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTATACGCCCGCTCACGCAGCCCGTTCAGGGt                                 >  8:159540/1‑148 (MQ=255)                                                                                                               cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCCCGGC‑CTTGCCCGCGAATAGGCCCGCACAGACACCACAAAACCGCAGTTTTTCTACGCCCGCTCACGCAGCCCTTTCAGGGt                                 >  2:116238/1‑148 (MQ=255)                                                                                                                 gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                               <  7:24613/148‑1 (MQ=255)                                                                                                                           taaaaaagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTACGCCCCCTCACCCAGCCCGTTCAGGGTATTGAACGgcg                      >  2:161141/6‑147 (MQ=255)                                                                                                                            cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCa                    <  2:5804/148‑1 (MQ=255)                                                                                                                              aaaagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAt                   >  8:152185/3‑147 (MQ=255)                                                                                                                                 agGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac               <  7:53518/148‑1 (MQ=255)                                                                                                                                          tCGGGCCGCC‑TTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGC‑GTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt      <  5:107265/148‑1 (MQ=255)                                                                                                                                                gctgttttttttttttTGCGGGGCCCGCCTCGGCGGGGCGCTGGCCCGCTTCGCGCCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg  <  6:126973/144‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740658‑4740947

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑