Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,685,542	(C)6→5	coding (1047/1404 nt)	fumC ←	fumarase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,685,537	0	C	.	100.0%	173.3 / NA	41	coding (1052/1404 nt)	fumC	fumarase C
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (19/22);  total (19/22)

GCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685400‑1685668                                                                                                                                          |                                                                                                                                    gCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                                   <  2:374666/140‑1 (MQ=255)  cTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcccc                                                                                                                                  >  2:497465/1‑140 (MQ=255)            ggTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATc                                                                                                                       >  2:657577/1‑140 (MQ=255)             gTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCg                                                                                                                      >  2:54509/1‑140 (MQ=255)              ttCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGc                                                                                                                     <  1:54509/140‑1 (MQ=255)                 aaTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAc                                                                                                                  >  2:668493/1‑140 (MQ=255)                     cccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcg                                                                                                              >  1:199299/1‑140 (MQ=255)                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         <  2:834910/140‑1 (MQ=255)                           gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGcc                                                                                                                    <  1:674517/128‑1 (MQ=255)                           gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGcc                                                                                                                    >  2:674517/1‑128 (MQ=255)                            cgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCa                                                                                                       <  2:307698/140‑1 (MQ=255)                             gcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                      <  1:974957/140‑1 (MQ=255)                                      gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAca                                                                                             >  1:913176/1‑140 (MQ=255)                                       ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                            <  2:650430/140‑1 (MQ=255)                                       ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                                                            <  2:920432/140‑1 (MQ=255)                                        tAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                                                           <  1:587892/140‑1 (MQ=255)                                                tCCATGCCATCTGCCAGCAAGCACACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCGCCTGACAGCAGAGCATg                                                                                   >  2:889129/1‑140 (MQ=255)                                                 ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCAATCa                                                                                                    <  2:513496/122‑1 (MQ=255)                                                 ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCAATCa                                                                                                    >  1:513496/1‑122 (MQ=255)                                                    tGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                               <  1:668493/140‑1 (MQ=255)                                                      ccATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAAt                                                                             <  2:813914/140‑1 (MQ=255)                                                             ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                                      <  1:497465/140‑1 (MQ=255)                                                             ccAGCAAGCACACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                                      <  1:889129/140‑1 (MQ=255)                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  2:928632/140‑1 (MQ=255)                                                                  aaGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACtgt                                                                 <  1:785265/140‑1 (MQ=255)                                                                     cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCTCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                              >  2:898583/1‑140 (MQ=255)                                                                       cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                                                            >  1:794953/1‑140 (MQ=255)                                                                                        tgtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCa                                           <  1:319056/140‑1 (MQ=255)                                                                                          tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATg                                         >  1:835433/1‑140 (MQ=255)                                                                                             aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCAGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGAtt                                      >  2:39452/1‑140 (MQ=255)                                                                                               cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGa                                    >  2:844420/1‑140 (MQ=255)                                                                                                aCCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAg                                   >  1:705197/1‑140 (MQ=255)                                                                                                    ttGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGc                               <  2:840341/140‑1 (MQ=255)                                                                                                         cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGc                          >  2:743453/1‑140 (MQ=255)                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       <  1:898583/140‑1 (MQ=255)                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       >  2:561151/1‑140 (MQ=255)                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       >  2:608064/1‑140 (MQ=255)                                                                                                                 gttcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc                  >  2:258123/1‑140 (MQ=255)                                                                                                                  ttcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTcc                 <  1:743453/140‑1 (MQ=255)                                                                                                                   tcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                <  1:844420/140‑1 (MQ=255)                                                                                                                             gTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt      <  1:657577/140‑1 (MQ=255)                                                                                                                             gTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt      <  1:258123/140‑1 (MQ=255)                                                                                                                                 ccGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCa  >  1:841471/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                    GCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685400‑1685668

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: ‑