Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,740,886	T→A	17.7%	intergenic (+77/+675)	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,886	0	T	A	17.7%	115.7 / 11.9	46	intergenic (+77/+675)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base A (3/5);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCT  >  NC_002947/4740737‑4741035                                                                                                                                                      |                                                                                                                                                        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                                                                                                                                             >  1:1673314/1‑144 (MQ=39)                   gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                             <  2:1673314/144‑1 (MQ=255)                     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                                                                                                    >  1:829605/1‑145 (MQ=255)                                                         tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTccgc                                                                                                    <  2:829605/145‑1 (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                                                                                   <  1:1511898/119‑1 (MQ=255)                                                                                                    tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCtt                                                                                   >  2:1511898/1‑119 (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                                                   <  2:940237/116‑1 (MQ=255)                                                                                                                                       ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                                   <  2:940276/116‑1 (MQ=255)                                                                                                                                       ggCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                                   >  1:940276/1‑116 (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)                                                                                                                                                     cTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCAcc   >  2:844821/1‑150 (MQ=255)                                                                                                                                                      ttGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCt  >  1:495118/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                        GATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTT‑‑TTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGCACCT  >  NC_002947/4740737‑4741035

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