Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,741,277	C→T	intergenic (+468/+284)	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,741,277	0	C	T	79.6%	105.6 / 17.8	54	intergenic (+468/+284)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base T (23/20);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑T‑GC‑ATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCA  >  NC_002947/4741155‑4741403                                                                                                                             |                                                                                                                                    gTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGaatccag                                                                                                       <  2:888679/150‑6 (MQ=255)                 cctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGaca                                                                                      >  1:362082/1‑150 (MQ=255)                      ttGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcc                                                                                 >  2:1113983/1‑150 (MQ=255)                            gTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCAATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGc                                                                           >  1:598644/1‑150 (MQ=255)                             tCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑T‑GC‑ATTCGTGa                                                                                                             <  1:53243/116‑1 (MQ=255)                             tCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCa                                                                          >  1:1631800/1‑150 (MQ=255)                               gTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAgc                                                                        >  1:832052/1‑150 (MQ=255)                                        agcaCGGAATCACGCGGCGTCACCGGCCATGGCTTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGTCACAGGCCCCGGacac                                                                                     >  2:417113/1‑128 (MQ=255)                                        agcaCGGAATCACGCGGCGTCACCGGCCATGGCTTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGTCACAGGCCCCGGacac                                                                                     <  1:417113/128‑1 (MQ=255)                                          caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                 >  1:1071345/1‑146 (MQ=255)                                          caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                 <  2:1071345/146‑1 (MQ=255)                                           aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                            <  2:1617207/150‑1 (MQ=255)                                           aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                            <  2:1631800/150‑1 (MQ=255)                                             ggAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GCAAGTGTTGAT‑CCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGc                                                          <  2:857323/150‑1 (MQ=255)                                             ggAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGc                                                          >  2:1358990/1‑150 (MQ=255)                                             ggAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGACACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGc                                                          >  2:353452/1‑150 (MQ=255)                                               aaTCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCa                                                        <  2:832052/150‑1 (MQ=255)                                                 tCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTAT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCATACGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                      <  1:33042/150‑1 (MQ=255)                                                   aCGCGGCGTCACCGGCAATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                    <  2:598644/150‑1 (MQ=255)                                                     gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                  >  1:760656/1‑150 (MQ=255)                                                      cggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGTGgcgc                                                 <  1:298752/150‑1 (MQ=255)                                                           tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                            >  2:675384/1‑150 (MQ=255)                                                            cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                           >  1:470922/1‑150 (MQ=255)                                                              ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACgtcgtc                                                                               <  2:274648/112‑1 (MQ=255)                                                              ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcgtc                                                                               >  1:274648/1‑109 (MQ=255)                                                               cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                        >  2:1010641/1‑150 (MQ=255)                                                               cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                        >  1:1130094/1‑150 (MQ=255)                                                               cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                        >  1:850569/1‑150 (MQ=255)                                                                 gCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                      <  2:362082/150‑1 (MQ=255)                                                                 gCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTCCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                      <  2:1164829/150‑1 (MQ=255)                                                                         gTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcg                              >  1:1485922/1‑150 (MQ=255)                                                                          tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                             <  1:1113983/150‑1 (MQ=255)                                                                            ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                   <  2:1054634/110‑1 (MQ=255)                                                                            ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                   >  1:1054634/1‑110 (MQ=255)                                                                              gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGc                                                                           <  2:500991/100‑1 (MQ=255)                                                                              gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa                         <  2:470922/150‑1 (MQ=255)                                                                              gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCGTCCTGCAGCGCCTc                                                                   >  1:1741397/1‑108 (MQ=255)                                                                              gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGc                                                                           >  1:500991/1‑100 (MQ=255)                                                                              gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                   <  2:1741397/108‑1 (MQ=255)                                                                               aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                        >  2:1736021/1‑150 (MQ=255)                                                                               aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGTCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                       <  2:1065878/135‑1 (MQ=255)                                                                               aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGACTT‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                       >  1:1065878/1‑135 (MQ=255)                                                                                aCTGGGTATAGGTCCAGCTCTGGACACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACAGCCTCCTGCAACGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACCCCGCGCTAgg                       >  2:237715/1‑150 (MQ=255)                                                                                  tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                     <  1:12208/150‑1 (MQ=255)                                                                                     gTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                  <  2:1130094/150‑1 (MQ=255)                                                                                         aGGTACAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACt              >  1:716267/1‑150 (MQ=255)                                                                                           gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACGCCGCGATAGGCCTTGTACTTg            >  1:414198/1‑150 (MQ=255)                                                                                           gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg            >  2:1203289/1‑150 (MQ=255)                                                                                              cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctc                                                                               >  2:1018169/1‑80 (MQ=255)                                                                                              cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctc                                                                               <  1:1018169/80‑1 (MQ=255)                                                                                               aGCTCTGGTCACCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                  <  1:690896/140‑1 (MQ=255)                                                                                               aGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCTGCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                  >  2:690896/1‑140 (MQ=255)                                                                                                   ctGGTCACCTTCGTGCTTGAG‑CA‑GCTTC‑TGCGAC‑TGGT‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGc    <  1:387968/150‑1 (MQ=255)                                                                                                         aCCTTCGTGCTTGAGCAG‑CTTCT‑GCGACT‑T‑GC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCa  <  1:839902/146‑1 (MQ=255)                                                                                                                             |                                                                                                                                    GTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CCTCC‑GCCAAC‑T‑GC‑ATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCA  >  NC_002947/4741155‑4741403

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