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%	216.7 / NA	51	coding (1052/1404 nt)	fumC	fumarase C
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (19/32);  total (19/32)

TCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685398‑1685668                                                                                                                                            |                                                                                                                                        tCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAgcgc                                                                                                                                         >  1:692193/1‑140 (MQ=255)  cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                        <  1:83830/140‑1 (MQ=255)  cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                        >  2:1147372/1‑140 (MQ=255)    cTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcccc                                                                                                                                      >  1:20819/1‑140 (MQ=255)                 tCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGcc                                                                                                                        <  1:259216/140‑1 (MQ=255)                   aaTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAc                                                                                                                      >  1:1228300/1‑140 (MQ=255)                    aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcg                                                                                                                     <  2:962611/140‑1 (MQ=255)                    aTACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcg                                                                                                                     <  1:625486/140‑1 (MQ=255)                        ccACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcgt                                                                                                                 <  2:692193/140‑1 (MQ=255)                         cACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGtt                                                                                                                <  1:921631/140‑1 (MQ=255)                          aCTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcg                                                                                                                  <  2:345556/137‑1 (MQ=255)                          aCTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCAcgtcg                                                                                                                  >  1:345556/1‑137 (MQ=255)                            tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTccc                                                                                                             >  2:1120352/1‑140 (MQ=255)                            tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTccc                                                                                                             >  2:747679/1‑140 (MQ=255)                             gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑CCCCATGTTGATCGCCACGTCGTTCCCCa                                                                                                           >  2:254656/1‑140 (MQ=255)                                    gtttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                                                     <  1:1287934/140‑1 (MQ=255)                                    gtttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                                                     <  1:1150443/140‑1 (MQ=255)                                      ttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑Ga                                                                                               >  1:1161476/1‑140 (MQ=255)                                        gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GAca                                                                                             >  2:209612/1‑140 (MQ=255)                                        gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GAca                                                                                             <  2:1352925/140‑1 (MQ=255)                                         ttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GAcag                                                                                            <  1:1333553/140‑1 (MQ=255)                                           aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GAcagca                                                                                          <  1:812246/140‑1 (MQ=255)                                           aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GAcagca                                                                                          <  1:1147372/140‑1 (MQ=255)                                               cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGc                                                                                      <  1:220483/140‑1 (MQ=255)                                                          aTCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGc                                                                           <  2:180089/140‑1 (MQ=255)                                                           tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCCCCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGcc                                                                          <  2:1228300/140‑1 (MQ=255)                                                           tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑CCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATgtct                                                                         <  1:254656/140‑4 (MQ=255)                                                               ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTcac                                                                      <  2:20819/140‑1 (MQ=255)                                                                cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  1:397984/140‑1 (MQ=255)                                                                cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  1:1150491/140‑1 (MQ=255)                                                                    aaGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACtgt                                                                 <  2:1161476/140‑1 (MQ=255)                                                                     aGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACtgtg                                                                >  1:985558/1‑140 (MQ=255)                                                                       cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                              <  1:209612/140‑1 (MQ=255)                                                                       cgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                                                              >  1:1201534/1‑140 (MQ=255)                                                                                tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc                                                      <  1:410172/139‑1 (MQ=255)                                                                                tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc                                                      >  2:410172/1‑139 (MQ=255)                                                                                        aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                             <  2:637080/140‑1 (MQ=255)                                                                                         ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGc                                            >  1:343/1‑140 (MQ=255)                                                                                               aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                             <  2:1160538/133‑1 (MQ=255)                                                                                               aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCgg                                             >  1:1160538/1‑133 (MQ=255)                                                                                                   ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGc                                  >  2:446449/1‑140 (MQ=255)                                                                                                       tGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                              <  1:1120352/140‑1 (MQ=255)                                                                                                           cGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGc                          >  1:247737/1‑140 (MQ=255)                                                                                                              aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGACAGCAGAGCATGGTTAATGCCTCACACCGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCgg                           >  1:184484/1‑140 (MQ=255)                                                                                                              aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       <  1:753294/140‑1 (MQ=255)                                                                                                                gACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                              <  2:624470/131‑1 (MQ=255)                                                                                                                gACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGcc                              >  1:624470/1‑131 (MQ=255)                                                                                                                   gttcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc                  <  2:665893/140‑1 (MQ=255)                                                                                                                     tcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                <  1:592865/140‑1 (MQ=255)                                                                                                                        gttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc                  >  1:1310425/1‑135 (MQ=255)                                                                                                                        gttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc                  <  2:1310425/135‑1 (MQ=255)                                                                                                                              aGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTgaga       >  2:1246411/1‑140 (MQ=255)                                                                                                                               gTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt      <  2:761413/140‑1 (MQ=255)                                                                                                                               gTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAt      <  2:985558/140‑1 (MQ=255)                                                                                                                                   ccGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCa  >  2:387236/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                        TCGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCT‑‑‑‑GACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685398‑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: ‑