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	97.9%	144.9 / NA	49	intergenic (+56/+697)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (30/17);  minor base A (1/0);  total (31/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_002947	4,740,820	0	G	T	100.0%	170.7 / NA	49	intergenic (+57/+696)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (31/18);  total (31/18)

CAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGC  >  NC_002947/4740679‑4740958                                                                                                                                             ||                                                                                                                                            cAGCCCGTACAATAACTTCCCTCCCCTCCAACTCTACACCCGACAGCAGAAGCGAGACATCTTAACCATGAAATACAACAATTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGAACTCGGGAAGCC‑TTATTTATTTCTCATGTGcccg                                                                                                                                      >  4:99175/1‑146 (MQ=255)  aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                       <  3:56211/146‑1 (MQ=255)   gCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGcc                                                                                                                                    >  1:282516/1‑148 (MQ=255)        tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                                               >  3:189617/1‑148 (MQ=255)                cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                                      >  7:91973/1‑149 (MQ=255)                   ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                                   <  2:195121/149‑1 (MQ=255)                           gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                            <  6:271427/148‑1 (MQ=255)                           gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                            <  2:261060/148‑1 (MQ=255)                           gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGAC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                              >  1:261060/1‑146 (MQ=255)                                cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCg                                                                                                       >  6:29574/1‑148 (MQ=255)                                cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTATTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                                       >  4:35119/1‑148 (MQ=255)                                    acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAATAGGCCGGc                                                                                                     >  5:385439/1‑146 (MQ=255)                                     cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                                 >  3:293917/1‑117 (MQ=255)                                     cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                                                    >  6:309645/1‑130 (MQ=255)                                      acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                                                <  4:293917/117‑1 (MQ=255)                                      acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                                   <  5:309645/130‑1 (MQ=255)                                                  aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAcc                                                                                    >  1:307651/1‑149 (MQ=255)                                                     ggcgaattCGTTGCCCTGAACGCCCGAAAGTACGCTGCAGTAGCCTAACGCAAGAAGGCGGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑c                                                                                 <  3:99175/141‑1 (MQ=255)                                                      gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAg                                                                               >  2:290121/1‑149 (MQ=255)                                                           aTCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCT‑CAGgtt                                                                            >  3:171329/1‑147 (MQ=255)                                                            tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATACCGGGCATGCCCGCGCAGAGGCCGGCACACCCACCACAAAACCT‑CAGatttt                                                                          >  4:25451/1‑148 (MQ=255)                                                            tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGCCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTAATCGCGGCCTTGCCCGAGAAGACGCCGGCACAAACCCCACAAAACCA‑CAGttta                                                                           >  6:170556/1‑146 (MQ=255)                                                                gcggttGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGCAGAGGCCGGGACAGACACCCCAAAAACG‑CAGTTTTTCTcc                                                                      >  6:23756/6‑148 (MQ=255)                                                                  cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCTGGC‑GTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCGCGCGATGCGGCGGGCCCAGACGCCGCAAAACCG‑GAGtttttgtaagc                                                                    >  6:107853/1‑142 (MQ=255)                                                                   ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACATCGCGGCCATGCCCGCCCAGAGGCCGGCCCCGACCCCCCAAAACCG‑CAGTTTTTATcccccc                                                                  >  2:373562/1‑149 (MQ=255)                                                                      gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGCGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTATccgcccg                                                                 >  5:122178/1‑147 (MQ=255)                                                                       aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACCCCCCAAAACCT‑CAGTTTTTCTCCGCCCGCt                                                               >  3:260036/1‑148 (MQ=255)                                                                             cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAg                                                        >  2:57953/1‑149 (MQ=255)                                                                              gCAATTATGCTGCATTAGCCGATGTAAAGATGGTTGCCCTCGGGCAG‑CTTTTTTTATTTGTCTTGTGCGGGCTTTTTCGGGGCTTTGCCCGCTAATTGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGc                                                       <  2:100371/149‑1 (MQ=255)                                                                                   tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAACAGGCCGGCCCAGAGACCACAAAACCT‑CACTTTTTCTCCGCCCGCTCACGCCGCCCgtt                                                  >  7:100886/1‑149 (MQ=255)                                                                                    attatGCATTAGCCTAACGCACCTAGGGTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGACTCTTCGCGCCCTTGCCCGCGAAGAGGCCTGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCttt                                                  <  6:283126/144‑1 (MQ=255)                                                                                      ggtGCATTAGCCGAACGCAAGAAGGCTGCCCTCGGGCTGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAATAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                                 <  6:122178/145‑1 (MQ=255)                                                                                      gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                               <  1:57953/149‑1 (MQ=255)                                                                                      gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCGTGCCCCGGCAGCGGCCGGCCCAGCCACCACAAAACCA‑CAGGTTTTATGCGCCCCCGCACCCAGcccggtca                                                >  5:22932/1‑143 (MQ=255)                                                                                         gCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATAGCGGCCGTGCCCCCGAAAAGGCCCGCCCAGACACCCCCAAACCA‑CAGTTTTTTTACGCCCGCTCACGCCGCCCgttgacggg                                            >  8:141899/1‑143 (MQ=255)                                                                                          cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                                        >  4:287509/1‑73 (MQ=37)                                                                                           aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                                       <  3:287509/73‑1 (MQ=37)                                                                                           aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                          >  2:60732/1‑149 (MQ=255)                                                                                             tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                                             >  6:147644/1‑81 (MQ=37)                                                                                              aGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                            <  5:147644/81‑1 (MQ=37)                                                                                                  gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGGGAAGCGGGCCGCCCCGACCCCACCACACCT‑CGGGTTTTCTCCGCCCGCGACCGCCGCGCTTGAAGGGGATTGATCgg                                   >  8:262326/1‑149 (MQ=255)                                                                                                  gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                                    <  1:290121/149‑1 (MQ=255)                                                                                                   aaCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                  <  1:60732/149‑1 (MQ=255)                                                                                                        aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCG‑CAGTGTTTGTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCCTc                             >  7:159011/1‑149 (MQ=255)                                                                                                           acgccgccccGCGGGCAGCC‑TTTTTTATTGGCCCTGGCCCGCCCTCTGCGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac                          <  7:358975/142‑1 (MQ=255)                                                                                                                    cTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt                 <  2:307651/149‑1 (MQ=255)                                                                                                                                ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                               >  6:174519/1‑108 (MQ=255)                                                                                                                                 tttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                              <  5:174519/108‑1 (MQ=255)                                                                                                                                    ttATTTGTCCTGTGCCGGCCTATTAGCGGCCTTGCCCGCCAAGAGGCCGGCCCACAAACCACAAAACCT‑AAGTTTTTATACGCCCCCAAACGCACCCCTTTCACGGTATTGAACGGCGCATACACCACGAACTTGTTGGCCACCATAGc  >  2:141834/1‑149 (MQ=255)                                                                                                                                             ||                                                                                                                                            CAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCT‑CAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGC  >  NC_002947/4740679‑4740958

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