Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,685,614	C→T	V325V (GTG→GTA)	fumC ←	fumarase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,685,614	0	C	T	100.0%	196.8 / NA	56	V325V (GTG→GTA)	fumC	fumarase C
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/20);  total (36/20)

AATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCTTTCAAC  >  NC_000913/1685484‑1685749                                                                                                                                   |                                                                                                                                            aaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCg                                                                                                                                   <  2:63909/140‑1 (MQ=255)  aTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCgg                                                                                                                                  >  1:752670/1‑140 (MQ=255)    tgtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCa                                                                                                                                >  1:171877/1‑140 (MQ=255)      tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTccc                                                                                                                                     >  2:322377/1‑133 (MQ=255)      tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTccc                                                                                                                                     <  1:322377/133‑1 (MQ=255)      tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATg                                                                                                                              >  1:1180525/1‑140 (MQ=255)      tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATg                                                                                                                              >  2:427488/1‑140 (MQ=255)           cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCa                                                                                                                                <  1:953583/133‑1 (MQ=255)           cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCa                                                                                                                                >  2:953583/1‑133 (MQ=255)           cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGC‑ACCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGa                                                                                                                         >  1:1198023/1‑140 (MQ=255)                 tGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGa                                                                                                                             >  2:147934/1‑130 (MQ=255)                 tGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGa                                                                                                                             <  1:147934/130‑1 (MQ=255)                   gACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCg                                                                                                                 >  2:1201369/1‑140 (MQ=255)                        aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                                                                                                            >  1:255672/1‑140 (MQ=255)                        aaGACGTTCAGTTCAAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                                                                                                            >  2:238113/1‑140 (MQ=255)                             gttcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc                                                                                                       <  2:255672/140‑1 (MQ=255)                             gttcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTc                                                                                                       >  1:1112008/1‑140 (MQ=255)                               tcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                                                                                                     <  1:1201369/140‑1 (MQ=255)                               tcagttcaAAGTTACCGGAAGCG‑CCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                                                                                                     <  1:393539/140‑1 (MQ=255)                                                   agcgcccccATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAAt                                                                              >  1:1203640/5‑140 (MQ=255)                                                                  aTCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGcc                                                               >  1:1079522/1‑140 (MQ=255)                                                                  aTCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGcc                                                               <  1:958338/140‑1 (MQ=255)                                                                      ccACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAgcgc                                                                    <  2:916475/131‑1 (MQ=255)                                                                      ccACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAgcgc                                                                    >  1:916475/1‑131 (MQ=255)                                                                      ccACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCagag                                                           >  2:639346/1‑140 (MQ=255)                                                                      ccACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCagag                                                           >  2:96462/1‑140 (MQ=255)                                                                        aCGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGgc                                                         <  1:96462/140‑1 (MQ=255)                                                                            cgtTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGgccagc                                                     >  1:384069/1‑140 (MQ=255)                                                                                   cATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGa                                              >  1:631806/1‑140 (MQ=255)                                                                                    aTCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGgccagcca                                                   >  2:1098499/1‑134 (MQ=255)                                                                                    aTCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGgccagcca                                                   <  1:1098499/134‑1 (MQ=255)                                                                                    aTCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGAc                                             <  2:1036371/140‑1 (MQ=255)                                                                                    aTCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGAc                                             <  2:526804/140‑1 (MQ=255)                                                                                      cACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGTCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGAcat                                           >  2:851008/1‑140 (MQ=255)                                                                                      cACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGa                                              >  2:630619/1‑137 (MQ=255)                                                                                      cACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGa                                              <  1:630619/137‑1 (MQ=255)                                                                                        ccTGACAGCAGAGCATGGTTAATGCCTCACGCTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGAcatca                                         <  1:677842/140‑1 (MQ=255)                                                                                               gcagAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGa                                  <  2:631806/140‑1 (MQ=255)                                                                                                cagAGCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGAt                                 >  1:562114/1‑140 (MQ=255)                                                                                                   agCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGa                                  >  1:894061/1‑140 (MQ=255)                                                                                                   agCATGGTTAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGAtttt                              >  2:857223/1‑140 (MQ=255)                                                                                                          ttAATGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGt                       >  1:737335/1‑140 (MQ=255)                                                                                                             aTGCCTCACACTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGAc                    <  2:737335/140‑1 (MQ=255)                                                                                                                   cacaCTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGcc              <  2:490301/140‑1 (MQ=255)                                                                                                                    acacTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTtcat                           <  2:855454/126‑1 (MQ=255)                                                                                                                    acacTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTtcat                           >  1:855454/1‑126 (MQ=255)                                                                                                                    acacTGTGTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCATTGACGCAGCCa             >  1:1092975/1‑140 (MQ=255)                                                                                                                        tgtgTTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAAccc         >  1:424155/1‑140 (MQ=255)                                                                                                                          tgtTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGcc              <  1:111500/133‑1 (MQ=255)                                                                                                                          tgtTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGcc              >  2:111500/1‑133 (MQ=255)                                                                                                                          tgtTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCtt       >  1:1190642/1‑140 (MQ=255)                                                                                                                          tgtTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCtt       >  1:345056/1‑140 (MQ=255)                                                                                                                          tgtTGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCtt       >  2:952053/1‑140 (MQ=255)                                                                                                                            ttGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCTTTc     >  2:658872/1‑140 (MQ=255)                                                                                                                             tGGGTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCTTTCa    >  2:596218/1‑140 (MQ=255)                                                                                                                               ggTTTACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCTTTCAAc  >  1:337099/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                            AATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGG‑‑‑‑GATTGAGATTTCACCAATTCCGCAGCGCGGGCCAGAGGCCAGCCAGCGGACATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCTTTCAAC  >  NC_000913/1685484‑1685749

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