Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,686,252	G→T	R113R (CGG→AGG)	fumC ←	fumarase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,686,252	0	G	T	100.0%	170.7 / NA	48	R113R (CGG→AGG)	fumC	fumarase C
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (21/27);  total (21/27)

CTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCT  >  NC_000913/1686115‑1686386                                                                                                                                          |                                                                                                                                       cttgcGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGt                                                                                                                                      >  1:639241/1‑140 (MQ=255)    gcgcAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATTGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAg                                                                                                                                   >  2:158068/1‑140 (MQ=255)             cagcagCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACtt                                                                                                                          <  2:1131947/140‑1 (MQ=255)             cagcagCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACtt                                                                                                                          <  1:756492/140‑1 (MQ=255)                 agcGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCgtt                                                                                                                      <  2:881845/140‑1 (MQ=255)                    gccgcAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCgttcat                                                                                                                   <  2:706118/140‑1 (MQ=255)                    gccgcAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCgttcat                                                                                                                   <  2:810230/140‑1 (MQ=255)                    gccgcAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCgttcat                                                                                                                   <  2:938052/140‑1 (MQ=255)                      cgcAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCgttcatgt                                                                                                                 <  2:243138/140‑1 (MQ=255)                        cAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCgttcatgttc                                                                                                               <  2:811774/140‑1 (MQ=255)                        cAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCgttcatgttc                                                                                                               <  2:1019922/140‑1 (MQ=255)                        cAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCgtccatgttc                                                                                                               >  1:954561/1‑140 (MQ=255)                                        cGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGc                                                                                               <  1:937314/140‑1 (MQ=255)                                              gACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGc                                                                                         <  1:912101/140‑1 (MQ=255)                                                cATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCg                                                                                       <  1:905551/140‑1 (MQ=255)                                                  tCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGt                                                                                     <  1:608534/140‑1 (MQ=255)                                                   cGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTc                                                                                    >  1:996722/1‑140 (MQ=255)                                                    gTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCt                                                                                   <  1:634445/140‑1 (MQ=255)                                                     ttGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTg                                                                                  <  1:158068/140‑1 (MQ=255)                                                          aCTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGa                                                                             <  2:949822/140‑1 (MQ=255)                                                           cTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGAt                                                                            >  1:32922/1‑140 (MQ=255)                                                              tGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGc                                                                         >  1:995668/1‑140 (MQ=255)                                                                gCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCa                                                                       >  2:1035982/1‑140 (MQ=255)                                                                 cTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAg                                                                      <  1:1035982/140‑1 (MQ=255)                                                                                     aGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCt                                                  <  2:728876/140‑1 (MQ=255)                                                                                       ggTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGt                                                <  2:996722/140‑1 (MQ=255)                                                                                       ggTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGt                                                <  2:995668/140‑1 (MQ=255)                                                                                        gTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTc                                               >  2:927617/1‑140 (MQ=255)                                                                                           aaCTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCAt                                                     <  2:1143343/131‑1 (MQ=255)                                                                                           aaCTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCAt                                                     >  1:1143343/1‑131 (MQ=255)                                                                                              tttACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCa                                         >  1:737688/1‑140 (MQ=255)                                                                                                tACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCt                                                  <  2:1031901/129‑1 (MQ=255)                                                                                                tACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCt                                                  >  1:1031901/1‑129 (MQ=255)                                                                                                 aCGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTa                                      >  2:1061555/1‑140 (MQ=255)                                                                                                     tCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTc                                  >  1:767980/1‑140 (MQ=255)                                                                                                              cgcACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCcgccgc                         <  2:44221/140‑1 (MQ=255)                                                                                                              cgcACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCcgccgc                         >  2:534296/1‑140 (MQ=255)                                                                                                              cgcACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCcgccgc                         >  2:297390/1‑140 (MQ=255)                                                                                                               gcACACCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCgccgcc                        >  2:857079/1‑140 (MQ=255)                                                                                                                   acCGCCGAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGAc                    >  1:259481/1‑140 (MQ=255)                                                                                                                     cgccgAGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGa                  >  2:754808/1‑140 (MQ=255)                                                                                                                          aGTAATTCACTGGCCCTGTTAGCCAGCACTTCGTCCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGgc             >  2:554346/1‑140 (MQ=255)                                                                                                                           gTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGgcg            <  2:737688/140‑1 (MQ=255)                                                                                                                           gTAATTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGgcg            <  2:259481/140‑1 (MQ=255)                                                                                                                              aTTCACTGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGcgctc         <  2:525826/140‑1 (MQ=255)                                                                                                                                   cTGGCCCTGTTAGCCAGCACTTCGTCCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCGTCCGCCGCCTGACGAATGGCGCTCGCttt    <  2:954561/140‑1 (MQ=255)                                                                                                                                    tGGCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTtc   >  1:79103/1‑140 (MQ=255)                                                                                                                                     ggCCCTGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGTTCTGCCAGATAGCCAGCGGGAATTCGTCGCCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTtct  >  1:554506/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       CTTGCGCAGCGCCAGCAGCGCCGCAACGTGCATCGCCGTCGGAAAGACATCGTTGGAACTTTGGCTTTTGTTCACGTCGTCGTTAGGGTGAACTTTACGTTCCATCCCGCGCACACCGCCGAGTAATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCT  >  NC_000913/1686115‑1686386

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑