Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,740,884	C→G	19.6%	intergenic (+75/+677)	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,884	0	C	G	19.6%	110.9 / 19.4	45	intergenic (+75/+677)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (18/18);  new base G (3/6);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAC  >  NC_002947/4740737‑4741033                                                                                                                                                    |                                                                                                                                                        gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAATGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                                                                       >  2:1306449/1‑150 (MQ=255)   tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                                                                     >  1:162635/1‑150 (MQ=255)           ccGCGCCAGCTGTACACCGGCGAGCAGCAGCCCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                             <  2:2133415/150‑1 (MQ=255)             gcgcCAGCTGTAAACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCCAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                           <  2:2012912/150‑1 (MQ=255)              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                                                          <  1:1533871/150‑1 (MQ=39)              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                                                          <  2:1213353/150‑1 (MQ=39)              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                                                          <  1:483664/150‑1 (MQ=255)                ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                                                                        <  2:162635/150‑1 (MQ=255)                  aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                                                                      <  1:1306449/150‑1 (MQ=255)                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                           <  2:1673314/144‑1 (MQ=255)                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                           >  1:1673314/1‑144 (MQ=39)                     tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCCAGGCCGCGAAGAGGCCGGca                                                                                                                                   <  1:1213425/150‑1 (MQ=255)                          acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGac                                                                                                                              >  1:1149849/1‑150 (MQ=255)                                  gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcg                                                                                                                                                 >  1:1839061/1‑123 (MQ=255)                                  gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcg                                                                                                                                                 <  2:1839061/123‑1 (MQ=255)                                     gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                           >  1:2278748/1‑126 (MQ=255)                                     gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                           <  2:2278748/126‑1 (MQ=255)                                                         tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTccgc                                                                                                  <  2:829605/145‑1 (MQ=255)                                                         tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTccgc                                                                                                  >  1:829605/1‑145 (MQ=255)                                                             aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACg                                                                                         <  2:1149849/150‑1 (MQ=255)                                                                         aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                             <  2:998262/150‑1 (MQ=255)                                                                              cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                                                        >  2:2187870/1‑150 (MQ=255)                                                                                                    tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCtt                                                                                 >  2:1511898/1‑119 (MQ=255)                                                                                                    tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCtt                                                                                 <  1:1511898/119‑1 (MQ=255)                                                                                                        cTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCAGATGTCGCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttg                                              >  1:1969961/1‑150 (MQ=255)                                                                                                                 ccTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                              >  1:980420/1‑109 (MQ=255)                                                                                                                 ccTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCa                                                                              <  2:980420/109‑1 (MQ=255)                                                                                                                 ccTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc                                     <  1:2187870/150‑1 (MQ=255)                                                                                                                                cTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTa                      >  2:1682557/1‑150 (MQ=255)                                                                                                                                      cGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                            <  2:2113080/90‑1 (MQ=25)                                                                                                                                      cGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGACCTTTCAgg                                                                            >  1:2113080/1‑90 (MQ=21)                                                                                                                                      cGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCt                >  1:1964460/1‑150 (MQ=255)                                                                                                                                       ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                                 >  1:940276/1‑116 (MQ=255)                                                                                                                                       ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                                 <  2:940237/116‑1 (MQ=255)                                                                                                                                       ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                                 <  2:940276/116‑1 (MQ=255)                                                                                                                                       ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                                 >  1:940237/1‑116 (MQ=255)                                                                                                                                          ctctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGCCA‑‑TGAGGGCCTATGGTGCTTTCTCCGCCCGCTCAGATAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt            >  1:414391/1‑150 (MQ=255)                                                                                                                                           tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                                                       <  1:689459/90‑1 (MQ=37)                                                                                                                                           tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTccgc                                                                                                  >  2:1180654/1‑63 (MQ=255)                                                                                                                                           tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTccgc                                                                                                  <  1:1180654/63‑1 (MQ=255)                                                                                                                                           tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCt                >  2:1397926/1‑145 (MQ=255)                                                                                                                                           tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCt                <  1:1397926/145‑1 (MQ=255)                                                                                                                                           tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc           >  2:2219231/1‑150 (MQ=255)                                                                                                                                            ctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac                                                     >  1:591473/1‑107 (MQ=38)                                                                                                                                                gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt      <  1:2219231/150‑1 (MQ=255)                                                                                                                                                    ccTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAc  >  1:243929/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                        GATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAC  >  NC_002947/4740737‑4741033

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