Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,740,844	C→G	22.5%	intergenic (+81/+672)	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,844	0	C	G	22.5%	68.7 / 19.0	40	intergenic (+81/+672)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (17/14);  new base G (3/6);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.01e-01

GGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGG  >  NC_002947/4740695‑4740992                                                                                                                                                       |                                                                                                                                                     ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGc                                                                                                                                                      <  7:48944/149‑1 (MQ=255)   ccGTCCGCGCCAGCTGTACACCGGCTAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGccc                                                                                                                                                    <  5:268913/149‑1 (MQ=255)        cgcgCCAGCTGTACAGCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAACTCAGCTCCAGAAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGaa                                                                                                                                               <  2:330529/149‑1 (MQ=255)        cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGC‑AAGGCCGCGaa                                                                                                                                               <  1:271945/149‑1 (MQ=255)        cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGaa                                                                                                                                               <  7:347066/149‑1 (MQ=255)         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGC‑AAGGCCGCGAag                                                                                                                                              <  1:371339/149‑1 (MQ=255)          cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGC‑AAGGCCGCGAaga                                                                                                                                             <  6:116072/149‑1 (MQ=38)          cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAaga                                                                                                                                             <  4:46668/149‑1 (MQ=255)          agcCATCTGTACACCGGCGAGCAGAAGCTCGACATCGTGGCCCTGAAAGACCGAAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCTCGGGC‑AAGGCCGCGAaga                                                                                                                                             <  4:335/148‑1 (MQ=255)               gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCg                                                                                                                                        <  2:68832/149‑1 (MQ=255)                   tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGc                                                                                                                                      <  4:310219/147‑1 (MQ=255)                        cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGC‑AAGGCCGCGAaga                                                                                                                                             <  4:230922/135‑1 (MQ=38)                        cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGC‑AAGGCCGCGAaga                                                                                                                                             >  3:230922/1‑135 (MQ=38)                                 gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAgg                                                                                                                                           >  2:40043/1‑128 (MQ=255)                                 gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAgg                                                                                                                                           <  1:40043/128‑1 (MQ=255)                                           aTCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtttt                                                                                                             >  8:2878/1‑148 (MQ=255)                                                   ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgccc                                                                                                    >  2:257417/1‑149 (MQ=255)                                                           aCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACGGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                                                                             >  2:232089/1‑148 (MQ=255)                                                             cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                                                                          >  6:340861/1‑149 (MQ=255)                                                                        tGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑ATGCCCGCGAAGCGGCCGGCACCGACACCACAAAACCACAGGTTTTCTCCGCCCGCTCACCCAGCCCGTTGAggg                                                                               >  5:218662/1‑149 (MQ=255)                                                                            gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCC‑TTGGCCCGGAAGAGGGCGGCACAGACACCACAAAACCGCACTTTTTATACGCCCCCGCACGGCGCCCGTTAAgggagat                                                                           >  6:17862/1‑145 (MQ=255)                                                                               gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGGC‑GTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTATCCGCCCGCGAACGCAGCCCGTTAAGGGTATTGa                                                                         >  4:212738/1‑148 (MQ=255)                                                                                 cGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                                                                      <  7:2878/149‑1 (MQ=255)                                                                                   aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACggt                                                                    <  1:257417/149‑2 (MQ=255)                                                                                      gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTT‑TTTGTCGTGTGCCGGCCTCTTCGCGGCG‑TTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                                <  3:128804/134‑1 (MQ=255)                                                                                                              ccttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCCCGAAGAGGCCGGCACACACACCACAAAACCTCAGTTTTTATCCGCCCGCGCACGCAGCCCTTTCAGGGGAATGTCAGGCGCATCAACCACGCACGTGTTGGcccc                                          >  4:216418/3‑147 (MQ=255)                                                                                                                ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGGATTGAACGGCGCATCCACCACGAACTTGTTGGcccc                                          >  1:247971/1‑146 (MQ=255)                                                                                                                            cctgTGCCGCCCTCTTCGCGGCC‑TGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGc                            <  3:212738/146‑1 (MQ=255)                                                                                                                                 gCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCCCAAAACCTCAGGTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTacca                       >  4:329960/1‑149 (MQ=255)                                                                                                                                    ggCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCACCGGCACGCTACCACCt                    >  3:17370/1‑149 (MQ=255)                                                                                                                                       ctctTCGCGGGC‑AAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                                 <  5:341333/85‑1 (MQ=37)                                                                                                                                       ctctTCGCGGGC‑AAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                                 >  6:341333/1‑85 (MQ=37)                                                                                                                                        tctTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt                >  7:99267/1‑149 (MQ=255)                                                                                                                                        tctTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt                >  7:93123/1‑149 (MQ=255)                                                                                                                                        tctTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTCTTGGCCACCATCGCCGGCACGCTACCACCTGGCt                >  6:348873/1‑149 (MQ=255)                                                                                                                                        tctTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACATCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCGCCCACGAACTTGTTGGCCACCATCGCCGGCACGCTCCCCCCTGGCt                >  2:374619/1‑149 (MQ=255)                                                                                                                                            cgcgGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg             <  2:247971/148‑1 (MQ=255)                                                                                                                                             gcgGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt            <  8:99267/148‑1 (MQ=255)                                                                                                                                                    gtGCCCGCGAATCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAcc     <  8:352836/148‑1 (MQ=255)                                                                                                                                                      gCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCCCGAACTTGTTGGCCACCAACGCCGGCACGCTACCACCTGGCTCGGTGTGAACCTg   >  7:248836/1‑149 (MQ=255)                                                                                                                                                      gCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg             >  6:28061/1‑139 (MQ=255)                                                                                                                                                      gCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg             <  5:28061/139‑1 (MQ=255)                                                                                                                                                       cccGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTgg  <  5:340861/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     GGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCTGG  >  NC_002947/4740695‑4740992

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