Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,381,680	A→G	G213G (GGT→GGC)	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,680	0	A	G	100.0%	111.4 / NA	33	G213G (GGT→GGC)	frdA	fumarate reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/14);  total (19/14)

TTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCG  >  NC_000913/4381544‑4381812                                                                                                                                          |                                                                                                                                      ttttGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCg                                                                                                                                    <  2:489564/140‑1 (MQ=255) ttttGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCg                                                                                                                                    >  1:549662/1‑140 (MQ=255)  tttGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGt                                                                                                                                   <  1:114478/140‑1 (MQ=255)  tttGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGt                                                                                                                                   <  1:723340/140‑1 (MQ=255)         aCCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGt                                                                                                                                   <  2:92598/133‑1 (MQ=255)         aCCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGt                                                                                                                                   >  1:92598/1‑133 (MQ=255)                     ccTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGt                                                                                                                            <  2:341873/128‑1 (MQ=255)                     ccTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGt                                                                                                                            >  1:341873/1‑128 (MQ=255)                     ccTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATgccgcc                                                                                                               >  2:114338/1‑140 (MQ=255)                     ccTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATgccgcc                                                                                                               >  1:286762/1‑140 (MQ=255)                       ttCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGt                                                                                                             <  2:286762/140‑1 (MQ=255)                       ttCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGt                                                                                                             <  1:108969/140‑1 (MQ=255)                        tCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATgccgc                                                                                                                <  1:144537/136‑1 (MQ=255)                        tCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATgccgc                                                                                                                >  2:144537/1‑136 (MQ=255)                        tCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGtt                                                                                                            >  1:132229/1‑140 (MQ=255)                        tCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGt                                                                                                             <  1:15989/139‑1 (MQ=255)                        tCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGt                                                                                                             >  2:15989/1‑139 (MQ=255)                          aCCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTgg                                                                                                          >  2:74307/1‑140 (MQ=255)                              cgGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTTATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGtgtt                                                                                                      <  2:411350/140‑1 (MQ=255)                                gCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGtgttgt                                                                                                    <  1:114338/140‑1 (MQ=255)                                gCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGtgttgt                                                                                                    <  2:569181/140‑1 (MQ=255)                                gCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGtgttgt                                                                                                    <  2:549662/140‑1 (MQ=255)                                              tcaGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGa                                                                                      <  1:74307/140‑1 (MQ=255)                                               caGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGAc                                                                                     >  1:406182/1‑140 (MQ=255)                                                     aCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCAc                                                                               >  1:332183/1‑140 (MQ=255)                                                     aCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCAc                                                                               >  2:473947/1‑140 (MQ=255)                                                            ccTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGt                                                                        >  1:278946/1‑140 (MQ=255)                                                                   gccgggTGGGGTGATAC‑GAACCAATTCCCATGTCACGCGGCGAACCGCTGTGTCTTGCCCCCATCCCCAT‑CCCCCCCCGGTAAAGGGTGCCGCTGTTGGTGTTTAAGCAAAAAACGCGCCCCGCACCGCCAGTAGCCATa                                                                 <  2:200570/136‑1 (MQ=255)                                                                                  cctaaaCGATTC‑CCATGTCACGCGGCGGACCGCGGTGGTTTAGCCCCATCCCCATCCGGTCCCGGGT‑AACGAGGCCGCTGTTGGTGTTGAAACGAAAACCGCCCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                                  <  2:152453/136‑1 (MQ=255)                                                                                  aCTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                                  >  1:213774/1‑140 (MQ=255)                                                                                   cTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCAc                                                 >  1:407306/1‑140 (MQ=255)                                                                                      aaCGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGa                                              >  2:124764/1‑140 (MQ=255)                                                                                         ctattcccATGCCCCGCGGGGGACCGCGGTGGTTTAGCCCCATCCCCAGCCCGCCCCGGGA‑AACGATGCCGCCGTTGGTGTTGTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATc                                            <  2:743637/134‑1 (MQ=255)                                                                                           aTT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                         >  1:125095/1‑140 (MQ=255)                                                                                             t‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCAc                                       >  1:743637/1‑140 (MQ=255)                                                                                                       cGCAGCGAAACGCGGGGGTTTGGCCCCACCCCCAGCCCGCCCCGGGA‑AACGAGCCCGCCGTTGGTGTTGAAACAAAAACCGCGCCCCGCCCCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGcc                              <  2:132229/140‑1 (MQ=255)                                                                                                         ggggggACCGCGGTGGCTTAGCCCCATCCCCATCCCGCCCCCGGA‑AACGATGCCGCCGTTGGTTTTGTAACAAAAACCGCGCCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCtt                            <  2:332183/136‑1 (MQ=255)                                                                                                         cAGCGGACCGCGGTGGCTTACCCCCATCCCCAGGCCGCCCCCGGA‑AACGATCCCCCCGTTGGTGTTGAAACAAAAACCGCCCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCtt                            <  2:278946/140‑1 (MQ=255)                                                                                                          aGCGGAACGCCGTGGCTTAGCGCCATACCCATGCCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTc                           >  2:253617/1‑140 (MQ=255)                                                                                                                                  acccccatgccGCCCCCGGA‑AACGATGCCGCGTTTGGTTTTGAAACAAAAACCGCCCCCCCCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGcc   <  2:429146/137‑1 (MQ=255)                                                                                                                                   tacctatgccGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCcg  >  1:80385/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                      TTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATT‑CCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGT‑AACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCTACCAGGCCG  >  NC_000913/4381544‑4381812

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