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.8%	195.2 / 6.0	48	intergenic (+42/+711)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base . (26/19);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 8.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGT  >  NC_002947/4740658‑4740945                                                                                                                                                    |                                                                                                                                               gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                                >  2:139585/1‑148 (MQ=255)        aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                         >  8:43206/1‑146 (MQ=255)                     ccAGCCCGTACAAGATCGGCCGTACGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCGGAAAGACCGCAAGTAAGCTGCAGTAGCCGACCGCAAGAAGGCTGCCCTCGAGCCGCC‑TTTGTTATTTGTCATGTGcc                                                                                                                            >  6:315882/1‑147 (MQ=255)                      cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                          >  7:93281/1‑148 (MQ=255)                       agctagTACATGATCGGCCGTCCGCGCCAGCTGTACACCGTCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCTCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGATGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                          <  5:315882/142‑1 (MQ=255)                       aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttt                                                                                                                                              <  2:418942/128‑1 (MQ=255)                       aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttt                                                                                                                                              >  1:418942/1‑128 (MQ=255)                            gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCc‑‑tct                                                                                                                   >  6:148389/1‑147 (MQ=255)                             tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                        <  2:52262/128‑4 (MQ=255)                             tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                        >  1:52262/1‑125 (MQ=255)                               cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTcgc                                                                                                               <  8:93281/148‑1 (MQ=255)                                aaGATCGGTCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTcg                                                                                                                <  6:179056/146‑1 (MQ=255)                                aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTcgc                                                                                                               <  3:106750/147‑1 (MQ=255)                                   aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGcc                                                                                                           >  1:410518/1‑148 (MQ=255)                                            ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCcgcg                                                                                                  <  7:707/148‑1 (MQ=255)                                             cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCcgcg                                                                                                  <  2:410518/147‑1 (MQ=255)                                               cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTAGCGGCCTTGCCcgc                                                                                                   >  6:171722/1‑144 (MQ=255)                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGGCAAGGCCGCGAag                                                                                               <  5:148389/147‑1 (MQ=255)                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGGCAAGGCCGCGAag                                                                                               <  6:25825/147‑1 (MQ=255)                                                 ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAagag                                                                                             >  7:175445/1‑148 (MQ=255)                                                     cTGGACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTACGCTGCACTAGCCGCACGCAAGACGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                          <  3:84617/147‑1 (MQ=255)                                                      tGTACACGGGCGAGCATACGCGCGACATCGTTGCCCTGAAAGCCCGCAAGTACGCTGCAGTAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                         <  6:348862/147‑1 (MQ=255)                                                        tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATCTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                         <  1:306179/146‑1 (MQ=255)                                                         acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTAGCGGCCTTTCCCGCGAAGAGGCCGGca                                                                                      >  2:132402/1‑147 (MQ=255)                                                            ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCATGCCCGCGAAGAGGCCGGccca                                                                                    >  5:141197/1‑146 (MQ=255)                                                                     ataaGCGCGACATCGTTGCCCTCAACCCCCGAAATTAAGCTGCTGTAGCCGAACCCAAGAAGGCTGCCCTGGGGCAG‑CTTTTTTTTTTTGTCGTGTGCCGGCC‑‑TCTTCGCGGCGTGGCCCGCTAAGAGGCCGGCACAGACACCACaaa                                                                         <  2:147423/146‑1 (MQ=255)                                                                        aGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCCCCGCCACCCCcaaa                                                                        >  3:387471/1‑146 (MQ=255)                                                                           gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGCACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCCGCC‑‑ACTTCCCGGCCATGCCCGCGAAGACGCCGCCCCAGACACCACAAAAccac                                                                    >  3:48183/1‑145 (MQ=255)                                                                           gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAcc                                                                      >  6:209947/1‑145 (MQ=255)                                                                           gcgACATCGTTCCCCTTAAAGCCCGCATGTAAGCTGCATTAGCCGAACGCATGATGGCTGCCCTGGGGATGCGTTTTTTT‑TTTGTCCTGTGCGGGCG‑‑TCTTCGCGGCGTTGCCCGGTATTAGGCCGGCACAGACACCACAAAACCTc                                                                    <  6:141197/147‑1 (MQ=255)                                                                                aTCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGc                                                                                                                         <  4:171337/91‑1 (MQ=255)                                                                                aTCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGc                                                                                                                         >  3:171337/1‑91 (MQ=255)                                                                                aTCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTAGCGGCCTTGCCCGCGAAGAGGGCGGCACAGACACCACAAAACCTCAGttt                                                               >  4:223621/1‑147 (MQ=255)                                                                                         cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCATAGCGGCCGTGCCCGCGCAGCGGGCGGCACAGACACCCCAAAACCACAGTTTTTGTacgccc                                                     >  1:415437/1‑148 (MQ=255)                                                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                    >  7:206380/1‑148 (MQ=255)                                                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGAGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                    >  1:106788/1‑148 (MQ=255)                                                                                             aaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTc                                                 >  7:272067/1‑148 (MQ=255)                                                                                                 ccGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCTTATTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAc                                               <  1:132402/148‑1 (MQ=255)                                                                                                  cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCGTGCCCGCGGAGAGGCCGGCACCGACACCACAAAAACTCAGTTTTTCTCCGCCCGCGCACGc                                             >  4:235083/1‑147 (MQ=255)                                                                                                      aGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TATTCGCGGCCTTGCCCGCCAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAccc                                         >  2:348235/1‑147 (MQ=255)                                                                                                          aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGCCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                    <  1:139585/148‑1 (MQ=255)                                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                  <  7:43206/148‑1 (MQ=255)                                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                  <  7:81128/148‑1 (MQ=255)                                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑ACATAGCGGCCTTGCCCGCGAAGACGCCGGCACAGACACCACAAAACCACAGGTTTTATCCGCCCCCTCACGCACCCCattaa                                   >  3:149676/1‑145 (MQ=255)                                                                                                             tGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                  <  5:209947/147‑1 (MQ=255)                                                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGc                                                                                                            <  2:418417/70‑1 (MQ=37)                                                                                                                aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGc                                                                                                            >  1:418417/1‑70 (MQ=37)                                                                                                                      cGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TATTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGCCACCACAAACCCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGaa                         >  6:95901/1‑147 (MQ=255)                                                                                                                       gAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGCCACCACAAAACCTGCGTTTTTCTCCGCCCGCGAACGCAGCCCTTTCAGGGTTTTGAACg                       >  2:57146/1‑148 (MQ=255)                                                                                                                        aaCGCACGAAGGCTGCCCGCGGGCAGCC‑TTTTTTATTTCTCCTGTGCCGGCC‑‑GCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACgg                      <  8:110519/148‑1 (MQ=255)                                                                                                                          cGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                     >  6:195357/1‑147 (MQ=255)                                                                                                                              agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCCGCCGTGCCCGCGAAGCGGCCGGCACAGACCCCCCAAAACCACCGGTTTTCTCCGCCCGCTCACGCAGCCCTTTGAGGGTATTTAACGGCGCAt                 >  5:361903/1‑147 (MQ=255)                                                                                                                                  ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacc            <  8:272067/148‑1 (MQ=255)                                                                                                                                             ggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTGACGCAGCCCTTTCAGGGTATTGAACGGCGCCTCCACCCCGAACttgt  >  5:92585/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                               GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCC‑‑TCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGT  >  NC_002947/4740658‑4740945

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