Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,382,478	C→A	intergenic (‑160/‑165)	frdA ← / → epmA	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit/Elongation Factor P Lys34 lysyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,382,478	0	C	A	100.0%	160.8 / NA	41	intergenic (‑160/‑165)	frdA/epmA	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit/Elongation Factor P Lys34 lysyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/17);  total (24/17)

AGATTGTTTTTATCCCACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGCAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAGACTTCGTGCCCTTGTCA  >  NC_000913/4382330‑4382624                                                                                                                                                     |                                                                                                                                                   aGATTGTTTTTATCCCACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGaaa                                                                                                                                                  <  1:529603/151‑1 (MQ=255)     tGTTTTTATCCCACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAgg                                                                                                                                              >  1:757665/1‑151 (MQ=255)         tttATCCCACAGACACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGAGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGACTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCt                                                                                                                                           >  2:871661/1‑150 (MQ=255)            aTCCCACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTaaa                                                                                                                                        <  2:652933/150‑1 (MQ=255)            aTCCCACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAAt                                                                                                                                       <  2:302190/151‑1 (MQ=255)            aTCCCACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAAt                                                                                                                                       <  1:366005/151‑1 (MQ=255)              cccACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAAtt                                                                                                                                      >  1:242256/1‑150 (MQ=255)              cccACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAAtt                                                                                                                                      <  1:740839/150‑1 (MQ=255)              cccACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTa                                                                                                                                     <  2:774146/151‑1 (MQ=255)                   aGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACt                                                                                                                                >  2:54983/1‑151 (MQ=255)                     ccACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAAATGGACGAGATCGATTTTTTTAGTTCGAAAGAGGTTTAATTTATCACTGa                                                                                                                              <  2:738282/151‑1 (MQ=255)                            cTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTATTGCGAAAGAGGGCTAAATTATCACTGAAgatgat                                                                                                                       >  2:117877/1‑151 (MQ=255)                                   gTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGAttaattta                                                                                                                <  2:488524/151‑1 (MQ=255)                                                aGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCAt                                                                                                   >  1:430802/1‑151 (MQ=255)                                                     cGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATc                                                                                                  >  2:489341/1‑147 (MQ=255)                                                         cTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACAATCAGATCGTg                                                                                          >  2:208344/1‑151 (MQ=255)                                                           gCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTg                                                                                          <  2:812441/149‑1 (MQ=255)                                                            cGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCtt                                                                                       <  2:189038/151‑1 (MQ=255)                                                             gAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCttt                                                                                      <  1:798609/151‑1 (MQ=255)                                                              aaCGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCtttt                                                                                     >  2:934215/1‑151 (MQ=255)                                                              aaCGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCtttt                                                                                     <  2:101817/151‑1 (MQ=255)                                                                 gCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGA‑TTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCttttttt                                                                                  >  1:755854/1‑150 (MQ=255)                                                                    aTTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGt                                                                               >  1:627831/1‑151 (MQ=255)                                                                       ccACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGAttaatttaa                                                                                                               <  2:598422/116‑1 (MQ=255)                                                                       ccACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGAttaattta                                                                                                                >  1:598481/1‑115 (MQ=255)                                                                          cTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGc                                                                         >  2:61176/1‑151 (MQ=255)                                                                                cTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCg                                                                   >  2:717125/1‑151 (MQ=255)                                                                                 tttATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGc                                                                  >  2:526163/1‑151 (MQ=255)                                                                                          acaacaGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGAtt                                                         >  2:107074/1‑151 (MQ=255)                                                                                                      tCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAAt                                             >  1:715493/1‑151 (MQ=255)                                                                                                       cTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAAtt                                            >  2:752729/1‑151 (MQ=255)                                                                                                         gATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAAtt                                            <  2:317387/149‑1 (MQ=255)                                                                                                                aaGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGcc                                   >  1:214864/1‑151 (MQ=255)                                                                                                                  gTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGt                                 <  2:807446/151‑1 (MQ=255)                                                                                                                              aCGAGATCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGt                     <  2:286187/151‑1 (MQ=255)                                                                                                                                 agaTCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAg                   <  1:489394/150‑1 (MQ=255)                                                                                                                                  gaTCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAGAc                 >  1:321544/1‑151 (MQ=255)                                                                                                                                    tCGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAGACtt               >  1:661757/1‑151 (MQ=255)                                                                                                                                     cGATTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAGACTTc              <  2:635763/151‑1 (MQ=255)                                                                                                                                       aTTTTTTTAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAGACTTCGt            >  2:602778/1‑151 (MQ=255)                                                                                                                                              tAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAGACTTCGTGCCCTTg     >  2:108794/1‑151 (MQ=255)                                                                                                                                              tAGTGCGAAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAGACTTCGTGCCCTTg     >  1:351785/1‑151 (MQ=255)                                                                                                                                                  gcgaAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCACACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCCCGTTGCGAGTAGACTTCGTGCCCTTGTCa  >  2:16363/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AGATTGTTTTTATCCCACAGCCACGTACTTCAGGGTAAGTACCTGAAAGTTACGGTCTGCGAACGCTATTCCACTGCTCCTTTATAGGTACAACAGTATAGTCTGATGGATAAGTCTGAAATTTGACGAGATCGATTTTTTTAGTGCGCAAGAGGGCTAAATTATCACTGAAGATGATTAATTTAATTACTAAACCATCAGATCGTGCTTTTTTTAGTCACTGCCCACCGCTGTTTGATTCCTGCGTAAAATTTGTCTCGCCGTCGCGTTGCGAGTAGACTTCGTGCCCTTGTCA  >  NC_000913/4382330‑4382624

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑