Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,381,225	T→A	D365V (GAT→GTT)	frdA ←	fumarate reductase flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,381,225	0	T	A	100.0%	50.1 / NA	16	D365V (GAT→GTT)	frdA	fumarate reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (6/10);  total (6/10)

GCCAGACGGCCGAAGACCACCAGTTCCGCCAGGGAGTTAGAACCCAGACGGTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCC  >  NC_000913/4381099‑4381314                                                                                                                               |                                                                                          gCCAGACGGCCGAAGACCACCAGTTCCGCCAGGGAGTTAGAACCCAGACGGTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGttt                                                                                     <  2:80878/133‑1 (MQ=255)     gACGGCCGAATACCACCAGTTCCGCCAGGGAGTTAGAACCCAGACGGTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGAt                                                                                 <  2:513757/133‑1 (MQ=255)          ccGAAGACCACCAGTTCCGCCAGGGAGTTAGAACCCAGACGGTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGc                                                                            <  1:329850/133‑1 (MQ=255)            gAAGACCACCAGTTCCGCCAGGGAGTTAGAACCCAGACGGTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGccc                                                                          <  1:66080/133‑1 (MQ=255)                             ccAGGGAGATAGAACCCACACGCTTTGCACTGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGt                                                         <  1:131988/133‑1 (MQ=255)                                      tAGAACCCAGACGGTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGACCAAGGTGTAGTATGCGTTCGGACGTAc                                                >  2:131988/1‑133 (MQ=255)                                              aGACGGTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCg                                        <  2:31861/133‑1 (MQ=255)                                                   gTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCt                                   <  1:493269/133‑1 (MQ=255)                                                            gTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGa                          >  2:362481/1‑133 (MQ=255)                                                                     aaCAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCTCAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCa                 <  1:48110/133‑1 (MQ=255)                                                                            gAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACg          >  1:425847/1‑133 (MQ=255)                                                                            gAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTAAACCGGATCGACGCCAACGTACg          >  2:436018/1‑133 (MQ=255)                                                                            gAACATTCACCCACGGCGAACAGACCTTTAATGCGGGGTTCTCAGTTCTGAACGGATTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACCTACCGGAATCGGCTCTATAACCGGACCGACCTCAACGTACg          >  2:226167/1‑133 (MQ=255)                                                                                 ttCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTc     >  1:490120/1‑133 (MQ=255)                                                                                   cACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGc   <  1:279844/133‑1 (MQ=255)                                                                                   cACCCACGCCGAACACACTTTTAATGCGGGCTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGc   <  1:145998/133‑1 (MQ=255)                                                                                    aCCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGAACGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGcc  >  1:46722/1‑133 (MQ=255)                                                                                                                               |                                                                                          GCCAGACGGCCGAAGACCACCAGTTCCGCCAGGGAGTTAGAACCCAGACGGTTTGCACCGTGCAGACCAACAGAGGAACATTCACCCACGGCGAACAGACCTTTAATGCGGGTTTCACAGTTCTGATCGGTTTCGATACCGCCCATGGTGTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCC  >  NC_000913/4381099‑4381314

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