Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,740,819	2 bp→CT	intergenic (+56/+696)	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,819	0	T	C	95.7%	137.4 / ‑3.2	46	intergenic (+56/+697)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (22/22);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_002947	4,740,820	0	G	T	97.8%	148.3 / ‑5.8	46	intergenic (+57/+696)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (22/23);  total (22/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 = 9.85e-01

CGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCC  >  NC_002947/4740684‑4740950                                                                                                                                           ||                                                                                                                                   cGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                                       >  5:281698/1‑149 (MQ=255)    aCAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAACTGGGGGGGGACGAAGGAGGCAAGAAGGCTGCCATCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                                       >  1:140138/1‑149 (MQ=255)     cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                                   <  8:112002/149‑1 (MQ=255)     cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                                   <  6:317096/149‑1 (MQ=255)        gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCCCTTGGCgggc                                                                                                                >  5:293997/1‑147 (MQ=255)              ccGGCCGCGCCAGCTGTACTCCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                          <  6:323069/149‑1 (MQ=255)                 tCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                        >  1:138777/1‑148 (MQ=255)                   cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                     <  3:76119/149‑1 (MQ=255)                   cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                      <  6:114010/148‑1 (MQ=255)                   cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                      >  2:390938/1‑148 (MQ=255)                     cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                   <  5:47333/149‑1 (MQ=38)                      gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                  <  3:61422/149‑1 (MQ=255)                       ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAgg                                                                                                 <  6:281698/149‑1 (MQ=38)                         aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCTCTTAGCGGCCTTGCGCGCGTAGCGGcc                                                                                               >  4:415425/1‑149 (MQ=255)                          gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                              <  6:166698/149‑1 (MQ=255)                          gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                              <  1:390938/149‑1 (MQ=255)                          gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                              <  5:243306/149‑1 (MQ=255)                          gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                              <  2:138777/149‑1 (MQ=255)                          gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCCAAGAGGCCg                                                                                              >  2:117637/1‑149 (MQ=255)                                 acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                  <  3:320016/138‑1 (MQ=38)                                 acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                  >  4:320016/1‑138 (MQ=38)                                                        cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGGCCGCGAAGAGGCCGGCACAGCCACCCCCAAACCGCAGATATtc                                                                 >  4:92756/1‑148 (MQ=255)                                                         gTTGCCCTGAAAGACCGCATGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                                <  1:117637/148‑1 (MQ=255)                                                              ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCGCAGTTTTTATacgccc                                                          >  4:230782/1‑149 (MQ=255)                                                               cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCAGAGTTTTTATacgcccc                                                         >  7:351385/1‑148 (MQ=255)                                                                tcaaaGACCGCAAGTAAGCTTCAGTACCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCAGGCCTCATCGAGGACTTGAACGCGAATAGGCAGGCCCAGACACCAACAAACCTCAGTATATCAacgccca                                                         >  6:151225/3‑147 (MQ=255)                                                                          cAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                              >  8:109988/1‑149 (MQ=255)                                                                          cAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAAAGGCCGGCACAGACACCACAAAAACTCAGTTTTTCTACGCCCGCTCACGCAGcc                                              >  1:138168/1‑149 (MQ=255)                                                                                 gCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTGGCCGCGGATAGGCCGGCACAGAACCCACAAAACCGAAGGTTTTCTCCGCCCGCGCACCCAGCCCGTTCAg                                       >  7:363249/1‑149 (MQ=255)                                                                                  cTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGTGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                      <  7:109988/149‑1 (MQ=255)                                                                                  cTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                      <  7:336897/149‑1 (MQ=255)                                                                                   tGCAGTAGCCGA‑‑‑AAGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                     >  2:119203/1‑149 (MQ=255)                                                                                           ccGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCGCAGACCCCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTTAAc                             >  3:37898/1‑149 (MQ=255)                                                                                                      aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTGTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                     >  1:343986/1‑149 (MQ=255)                                                                                                        gatggCTGCCGTCGGGCAG‑CTTTTTTTATTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca                   <  2:138168/146‑1 (MQ=255)                                                                                                         aagccggCCCTCGGGCAGCC‑TTTTTTATTTGGCCTGTGCCGGCCGCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac                  <  4:37898/143‑1 (MQ=255)                                                                                                         aagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac                  >  5:106519/1‑149 (MQ=255)                                                                                                           ggcggcCCTCGGGCAGCC‑TTTTTAATTGGCCCTGTGCCGCCCTCTGCGCGCCGTGCCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca                <  3:415425/149‑1 (MQ=255)                                                                                                           ggCTGCCGTCGGCCAGCTTTTTTTTTTTGCCCTTG‑GCCGCCCTCTTGCCCGCCTTCCCCGCTAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca                <  4:219462/149‑1 (MQ=255)                                                                                                               gCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCAAGaa            >  1:329460/1‑149 (MQ=255)                                                                                                                cccTCGGGCAGCTTTTTTTTTTTGTCCTGG‑GCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc           <  3:145983/149‑1 (MQ=255)                                                                                                                     gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                 <  7:246489/122‑1 (MQ=39)                                                                                                                     gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                 >  8:246489/1‑122 (MQ=39)                                                                                                                     gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgt       <  2:140138/148‑1 (MQ=255)                                                                                                                             cttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcc  <  3:78284/145‑1 (MQ=255)                                                                                                                             cttttttATTTGTCCTGTGCCGGCCACTTAGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACACAGTTTTTATACGCCCCCGCACGCCGCCCGTTAAGGGTATTTAACGGCGCATACACCACCAACATGTTGGcc  >  4:78284/2‑146 (MQ=255)                                                                                                                                           ||                                                                                                                                   CGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGA‑‑‑ACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCC  >  NC_002947/4740684‑4740950

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