Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	765,135	T→C	13.5%	intergenic (+86/+18)	sucD → / ← mngR	succinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit/transcriptional repressor for the mannosyl‑D‑glycerate catabolic operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	765,135	0	T	C	13.5%	99.8 / 8.4	37	intergenic (+86/+18)	sucD/mngR	succinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit/transcriptional repressor for the mannosyl‑D‑glycerate catabolic operon
Reads supporting (aligned to +/- strand):  ref base T (18/14);  new base C (3/2);  total (21/16)
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.42e-01

ATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGTTCGCAGCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGCAGACGGCTCATCTCTTCGGTTGGCATCAGAGGAATGATTTCCTGTTCGCTACGATCGATAACCATTTTCTTCACTTCTTCGATAAAGTGATATTTCGAAT  >  NC_000913/764912‑765400                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                          aTGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGTTCGCAGCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTa                                                                                                                                                                                                                        >  4:79516/1‑275 (MQ=255)        aCGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGTTCGCAGCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATCAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACGTGGTTGCCCATCGTATGATGGCCTTTTTTGTTCTTATCGCGATTGTTTTCGCTGCGCTATCAGGGTAAATTtat                                                                                                                                                                                                                 >  4:337828/1‑275 (MQ=255)         cgcgGGTGCCACCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGTTCGCAGCCTGGCGGATATCGGTGAAGCACTGCAAACTGTTCTGAAATAAATATCTGTAATACGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTTATTATCGTGCTGATTTTCGCTGCGCTTTCAGGGTAAAttct                                                                                                                                                                                                                  >  1:55232/1‑271 (MQ=255)         ccgggggaccATCATGGCCGGTGGGAAAGGGACTGCGAATGAGAAATTTGCTGCTCTGGAAACTGCAGGCGTGAAGACCGTTCGCAGCGTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTtata                                                                                                                                                                                                                <  1:67325/267‑1 (MQ=255)                gcaatatttGCGGGGGGGAAAGGGGNGNGNGAGGAGAATTTGGCTGCTAGGGAAGCCGCAGGCGTGAAAACCGTTNGCAGCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCNCCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATc                                                                                                                                                                                                          <  3:54603/267‑1 (MQ=255)                                ggAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGTTCGCAGCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTGTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTg                                                                                                                                                                                          <  2:104360/274‑1 (MQ=255)                                                                          cAGGCGTGAAAACCGTTCGCAGCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATACAAACACCCGACCATCAACTAAATATCCAcgc                                                                                                                                                <  2:79517/274‑1 (MQ=255)                                                                                       cGTTCGCAGCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGTTATCAGGGTAAATTTATAGTCCTCGGTATTAAACGCGTTGCGGCTATATTCAAACACCCGACCACCAACTAAATAACCAAGCCAGCCttttta                                                                                                                                     >  2:527619/1‑271 (MQ=255)                                                                                               gCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCt                                                                                                                          <  2:543722/275‑1 (MQ=255)                                                                                                   ggCGGATATCGGCGAAGCACTGAAAACTGTTCTGAAATAAAAATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCGCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGt                                                                                                                      <  4:204906/275‑1 (MQ=255)                                                                                                                                   tGAAATAAATCTCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGTAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTatgat                                                                                                                                                                                                                                                                        >  7:374053/1‑97 (MQ=255)                                                                                                                                   tGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGTCTGTCGGTTTCGACATGGTTGGCCATCGTatgat                                                                                                                                                                                                                                                                        <  8:374053/97‑1 (MQ=255)                                                                                                                                        taaataTCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTatga                                                                                                                                                                                                                                                                         >  8:1059151/1‑91 (MQ=255)                                                                                                                                        taaataTCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTatga                                                                                                                                                                                                                                                                         <  7:1059151/91‑1 (MQ=255)                                                                                                                                                      aaGAAATAGCCCTCGCCGCTTCCCTCTACAGGAGTGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCATACGATGGc                                                                                                                                                                                                                                                                     <  6:145405/81‑1 (MQ=255)                                                                                                                                                      aaGAAATAGCCCTCGCCGCTTCCCTCTACAGGAGTGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCATACGATGGc                                                                                                                                                                                                                                                                     >  5:145405/1‑81 (MQ=255)                                                                                                                                                          aaTAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATCTTCAAACACCCGACCATCAACTAAATCTCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGAAGACGGCTAATCTTTTCGGTTGGCATCAGAGGAATGAtt                                                               >  1:103786/1‑275 (MQ=255)                                                                                                                                                            tAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTgcg                                                                                                                                                                                                                                   >  4:17319/1‑109 (MQ=255)                                                                                                                                                            tAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGCAGACGGCTCATCTCTTCGGTTGGCATCAGAGGAATGAttt                                                              >  2:321847/1‑274 (MQ=255)                                                                                                                                                                ccTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTAt                                                                                                                                                                                                                               <  2:17319/109‑1 (MQ=255)                                                                                                                                                                         ttCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTtgtg                                                                                                                                                                                                                                                           <  8:88230/72‑1 (MQ=255)                                                                                                                                                                         ttCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTtgtg                                                                                                                                                                                                                                                           >  7:88230/1‑72 (MQ=255)                                                                                                                                                                         ttCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGCAg                                                                                                    >  2:206124/1‑223 (MQ=255)                                                                                                                                                                             ctctACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGCAGACgg                                                                                                <  4:206123/223‑1 (MQ=255)                                                                                                                                                                                          gCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGCAGACGGCTCATCTCTTCGGTTGGCATCAGAGGAATGATTTCCTGTTCGCTACGATCGATAACCATTttctt                               >  2:317223/1‑275 (MQ=255)                                                                                                                                                                                                    gTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGAc                                                                                                                                                                    <  8:131546/132‑1 (MQ=255)                                                                                                                                                                                                    gTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGAc                                                                                                                                                                    >  7:131546/1‑132 (MQ=255)                                                                                                                                                                                                       ggTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTtat                                                                                                                                                                                                                 <  6:586759/84‑1 (MQ=255)                                                                                                                                                                                                       ggTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTtat                                                                                                                                                                                                                 >  5:586759/1‑84 (MQ=255)                                                                                                                                                                                                         tttCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTgcgc                                                                                                                                                                                                                                  >  8:10902/1‑65 (MQ=255)                                                                                                                                                                                                         tttCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTgcgc                                                                                                                                                                                                                                  <  7:10902/65‑1 (MQ=255)                                                                                                                                                                                                          ttCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGCAGACGCCTCATCTCTTCGGTTGGCATCCAAGGAATGATTTCCTGTTCGCTACGATCGATAACCATTTTCTTCACTTCTTCGATAAAg               >  2:319381/1‑275 (MQ=255)                                                                                                                                                                                                                        ccATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCCACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGCAGACGGCTCATCTCTTCGGTTGGCATCCGAGGAATGATTTTCTGTTCGCTTCGATCGAATACCATTTTCTTCCCTTCTTCACTAAAGTGCTTTTTCGAAt  >  4:294397/1‑274 (MQ=255)                                                                                                                                                                                                                        ccATCATACGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGc                                                                                                                                                                                      >  5:113980/1‑94 (MQ=255)                                                                                                                                                                                                                        ccATCATACGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGc                                                                                                                                                                                      <  6:113980/94‑1 (MQ=255)                                                                                                                                                                                                                        ccATCATACGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGaa                                                                                                                                >  5:908163/1‑148 (MQ=255)                                                                                                                                                                                                                         cATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCAcg                                                                                                                                                 <  7:703307/130‑1 (MQ=255)                                                                                                                                                                                                                         cATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCAcg                                                                                                                                                 >  8:703307/1‑130 (MQ=255)                                                                                                                                                                                                                             gTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCt                                                                                                                          >  4:551340/1‑149 (MQ=255)                                                                                                                                                                                                                              tATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAat                                                                                                                                                        >  7:34927/1‑118 (MQ=255)                                                                                                                                                                                                                              tATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAat                                                                                                                                                        <  8:34927/118‑1 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                          ATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGTTCGCAGCCTGGCGGATATCGGTGAAGCACTGAAAACTGTTCTGAAATAAATATCTGTAATAAGAAATAGCCCTCGCCGCTTCCCTCTACAGGAATGGCGAAGGGCTGTCGGTTTCGACATGGTTGGCCATCGTATGATGGCCTTTTTTGTGCTTATCGCGATGATTTTCGCTGCGCTATCAGGGTAAATTTATAGTCATCGGTATTAAAAGCGTTGCGGCTATATTCAAACACCCGACCATCAACTAAATATCCACGCGATACTTTTTCAAGAATCGGCTTTGTCTGGCTGATATTAAGCAGACGGCTCATCTCTTCGGTTGGCATCAGAGGAATGATTTCCTGTTCGCTACGATCGATAACCATTTTCTTCACTTCTTCGATAAAGTGATATTTCGAAT  >  NC_000913/764912‑765400

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑