Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,687,729	G→A	R217C (CGC→TGC)	fumA ←	fumarase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,687,729	0	G	A	100.0%	184.7 / NA	53	R217C (CGC→TGC)	fumA	fumarase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/22);  total (31/22)

TTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGG  >  NC_000913/1687601‑1687863                                                                                                                                 |                                                                                                                                       ttCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCaa                                                                                                                             >  1:746294/1‑140 (MQ=255)   ccGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAAcc                                                                                                                           >  1:121219/1‑140 (MQ=255)   ccGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGCGCACATCTTCTCAAcc                                                                                                                           >  1:852043/1‑140 (MQ=255)     gTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAg                                                                                                                         >  1:124055/1‑140 (MQ=255)     gTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAg                                                                                                                         >  2:529088/1‑140 (MQ=255)         gCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTaa                                                                                                                     >  2:27973/1‑140 (MQ=255)             ttCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAAtttt                                                                                                                 >  1:618076/1‑140 (MQ=255)                   atAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGtt                                                                                                           >  1:59878/1‑140 (MQ=255)                       tATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTcc                                                                                                        <  2:77468/139‑1 (MQ=255)                       tATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTcc                                                                                                        >  1:77468/1‑139 (MQ=255)                         tttCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCg                                                                                                     >  1:76480/1‑140 (MQ=255)                           tCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGc                                                                                                   >  1:201193/1‑140 (MQ=255)                            cgcgGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGc                                                                                                   >  2:510064/1‑139 (MQ=255)                            cgcgGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGc                                                                                                   <  1:510064/139‑1 (MQ=255)                            cgcgGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCg                                                                                                  >  1:30207/1‑140 (MQ=255)                             gcgGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGt                                                                                                 >  2:750033/1‑140 (MQ=255)                              cgGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTc                                                                                                >  2:462228/1‑140 (MQ=255)                              cgGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTc                                                                                                >  2:88131/1‑140 (MQ=255)                                       gTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCt                                                                                       <  2:746294/140‑1 (MQ=255)                                            aCCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGt                                                                                  <  2:30207/140‑1 (MQ=255)                                             ccGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt                                                                                 >  1:229457/1‑140 (MQ=255)                                                ttttAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTcc                                                                              <  1:462228/140‑1 (MQ=255)                                                ttttAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTcc                                                                              <  2:59878/140‑1 (MQ=255)                                                  ttAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTg                                                                            >  2:728119/1‑140 (MQ=255)                                                  ttAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTg                                                                            >  2:175885/1‑140 (MQ=255)                                                   tAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                                           <  1:728119/140‑1 (MQ=255)                                                   tAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                                           <  2:76480/140‑1 (MQ=255)                                                    aaGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGAt                                                                          <  2:719319/140‑1 (MQ=255)                                                     aGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATa                                                                         <  1:88131/140‑1 (MQ=255)                                                          cGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGAt                                                                    >  1:649375/1‑140 (MQ=255)                                                                       gTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGc                                                       >  2:557942/1‑140 (MQ=255)                                                                       gTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGc                                                       >  1:432689/1‑140 (MQ=255)                                                                                 tAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccacca                                             >  1:154223/1‑140 (MQ=255)                                                                                 tAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccacca                                             <  2:432689/140‑1 (MQ=255)                                                                                   acgaacgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccacc                                           <  2:124055/140‑1 (MQ=255)                                                                                        cgCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGCTTGCCGAACCACCACCTTTgg                                      <  1:529088/140‑1 (MQ=255)                                                                                                gATACGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACag                              <  1:27973/140‑1 (MQ=255)                                                                                                   aCGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTgg                                      <  2:592487/129‑1 (MQ=255)                                                                                                   aCGGAGGACAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTgg                                      >  1:592487/1‑129 (MQ=255)                                                                                                     ggaggaCAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACag                              >  2:568356/1‑135 (MQ=255)                                                                                                     ggaggaCAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACag                              <  1:568356/135‑1 (MQ=255)                                                                                                       aggaCAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAtt                       >  1:410436/1‑140 (MQ=255)                                                                                                         gaCAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                     >  1:491573/1‑140 (MQ=255)                                                                                                           cAGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTa                   >  1:737900/1‑140 (MQ=255)                                                                                                            aGGCCGCCGTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTAt                  >  2:181006/1‑140 (MQ=255)                                                                                                                  ccgTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgtcg            <  2:649375/140‑1 (MQ=255)                                                                                                                   cgTACCCAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGc           >  1:903924/1‑140 (MQ=255)                                                                                                                        ccAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCa      <  2:737900/140‑1 (MQ=255)                                                                                                                         cAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa     <  2:903924/140‑1 (MQ=255)                                                                                                                         cAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa     <  2:491573/140‑1 (MQ=255)                                                                                                                         cAGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa     <  2:618076/140‑1 (MQ=255)                                                                                                                          aGCGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAAc    <  2:532631/140‑1 (MQ=255)                                                                                                                            cGTGCACATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACgg  >  2:302660/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                       TTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCACCAATAACGAACGCAATATGATACGGAGGACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGG  >  NC_000913/1687601‑1687863

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