Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,346,977	T→A	K117I (AAA→ATA)	fumB ←	fumarase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,346,977	0	T	A	100.0%	104.2 / NA	30	K117I (AAA→ATA)	fumB	fumarase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (13/17);  total (13/17)

ATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTT  >  NC_000913/4346901‑4347054                                                                             |                                                                              aTCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtattt                                                                            >  2:798160/1‑80 (MQ=255) aTCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtattt                                                                            <  2:798110/80‑1 (MQ=255)  tCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtatttt                                                                           <  2:798111/80‑1 (MQ=255)     tCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTAcc                                                                        >  2:798157/1‑80 (MQ=255)      cGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACcga                                                                      >  1:798168/1‑81 (MQ=255)             ggTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACgatgat                                                               <  2:798114/81‑1 (MQ=255)             ggTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACgatgat                                                               <  2:798113/81‑1 (MQ=255)             ggTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACgatgat                                                               <  2:798112/81‑1 (MQ=255)              gtgtTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATc                                                              <  2:798115/81‑1 (MQ=255)               tgtTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATcg                                                             <  1:798116/81‑1 (MQ=255)                gtTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATcgc                                                            >  2:798165/1‑81 (MQ=255)                 ttATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATcgc                                                            <  1:798117/80‑1 (MQ=255)                    taGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGt                                                         <  1:798118/80‑1 (MQ=255)                        cGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCgg                                                    <  1:798119/81‑1 (MQ=255)                           cTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTAt                                                 >  2:798159/1‑81 (MQ=255)                            tttCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTTCCGGTATc                                                >  2:798175/1‑81 (MQ=255)                            tttCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATc                                                >  1:798162/1‑81 (MQ=255)                              tCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCt                                              >  2:798178/1‑81 (MQ=255)                                gACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTgg                                            >  2:798163/1‑81 (MQ=255)                                    gCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAgg                                        <  2:798120/81‑1 (MQ=255)                                      gTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTc                                      <  1:798121/81‑1 (MQ=255)                                          cttcATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCgg                                    >  1:798172/1‑79 (MQ=255)                                          cttcATCACCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGgc                                   <  1:798122/80‑1 (MQ=255)                                              atcaCCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCAc                              >  2:798170/1‑81 (MQ=255)                                                caCCGCCGCCGGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGc                            >  1:798167/1‑81 (MQ=255)                                                   cgccgccgGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCtt                         <  1:798123/81‑1 (MQ=255)                                                    gccgccgGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCttt                        <  1:798124/81‑1 (MQ=255)                                                      cgccgGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTgg                      >  1:798148/1‑81 (MQ=255)                                                       gccgGTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                     >  2:798176/1‑81 (MQ=255)                                                          ggTCCACACGCGCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                  <  2:798125/81‑1 (MQ=255)                                                                  cgcgCTGGCCTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCgg          <  1:798126/81‑1 (MQ=255)                                                                          ccTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCtt  <  1:798127/81‑1 (MQ=255)                                                                          ccTATTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCtt  <  1:798128/81‑1 (MQ=255)                                                                             |                                                                              ATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTT  >  NC_000913/4346901‑4347054

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