Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,178,383	C→T	R61H (CGT→CAT)	fbaB ←	fructose‑bisphosphate aldolase class I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,178,383	0	C	T	100.0%	176.0 / NA	44	R61H (CGT→CAT)	fbaB	fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/29);  total (15/29)

GCCACACAGTTACAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGACGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTCAGAAGGAATTGTCATACAACGGTGCTGTAAAAGGTTGTCG  >  NC_000913/2178236‑2178529                                                                                                                                                    |                                                                                                                                                   gCCACACAGTTACAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGcc                                                                                                                                                 <  1:453449/151‑1 (MQ=255)  ccACACAGTTACAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGccc                                                                                                                                                >  2:598219/1‑151 (MQ=255)   cacacaGTTACAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCg                                                                                                                                               <  1:246459/151‑1 (MQ=255)       caGTTACAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGgtgt                                                                                                                                           >  1:290361/1‑151 (MQ=255)        aGTTACAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGtgtt                                                                                                                                          <  1:620612/151‑1 (MQ=255)          ttACAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGtgttgt                                                                                                                                        >  1:413345/1‑151 (MQ=255)            aCAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTa                                                                                                                                       <  2:136361/150‑1 (MQ=255)                 ccGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGt                                                                                                                                 >  1:219853/1‑151 (MQ=255)                  cGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTc                                                                                                                                <  2:230305/151‑1 (MQ=255)                      tCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTg                                                                                                                              <  1:231836/149‑1 (MQ=255)                       cGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCat                                                                                                                           >  1:250526/1‑151 (MQ=255)                         aTCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCatat                                                                                                                         <  2:690449/151‑1 (MQ=255)                           cGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATAtt                                                                                                                        <  2:23835/150‑1 (MQ=255)                            gCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTAc                                                                                                                      <  2:576179/151‑1 (MQ=255)                               aGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTa                                                                                                                   <  1:378982/151‑1 (MQ=255)                                gTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTaa                                                                                                                  <  1:257508/151‑1 (MQ=255)                                  tCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAac                                                                                                                 >  1:289323/1‑150 (MQ=255)                                   caacaaTGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAacac                                                                                                               >  2:332173/1‑151 (MQ=255)                                     acaaTGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCg                                                                                                             <  1:542324/151‑1 (MQ=255)                                       aaTGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCt                                                                                                           >  2:221123/1‑151 (MQ=255)                                        aTGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTg                                                                                                          <  2:289311/151‑1 (MQ=255)                                              ttCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTggcggg                                                                                                     <  2:75083/150‑1 (MQ=255)                                                  ggTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGatt                                                                                                <  2:744383/151‑1 (MQ=255)                                                     cAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGattatt                                                                                             >  1:249734/1‑151 (MQ=255)                                                        aGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTc                                                                                          <  2:713541/151‑1 (MQ=255)                                                                 ggTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTAc                                                                                 >  1:33086/1‑151 (MQ=255)                                                                  gTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAGAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTAcg                                                                                >  1:8792/1‑151 (MQ=255)                                                                   ttAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGattatt                                                                                             >  2:407525/1‑137 (MQ=255)                                                                   ttAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGattatt                                                                                             <  1:407546/137‑1 (MQ=255)                                                                     agcagcAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTAcgcg                                                                              <  1:273736/150‑1 (MQ=255)                                                                        agcaAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCt                                                                          <  1:533833/151‑1 (MQ=255)                                                                                      cGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTcc                                                            <  2:351735/151‑1 (MQ=255)                                                                                           gagTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGga                                                       <  2:344609/151‑1 (MQ=255)                                                                                             gTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGgaga                                                     <  2:758682/151‑1 (MQ=255)                                                                                             gTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGgaga                                                     <  1:344335/151‑1 (MQ=255)                                                                                                   aaCGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGc                                               <  1:517069/151‑1 (MQ=255)                                                                                                          tGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTc                                          >  2:99512/1‑149 (MQ=255)                                                                                                          tGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTc                                          >  1:57305/1‑149 (MQ=255)                                                                                                           ggTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTCa                                         <  1:226850/149‑1 (MQ=255)                                                                                                                         aTAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACGCCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTCAGAAGGAATTGTCATAc                         <  2:76333/151‑1 (MQ=255)                                                                                                                         aTAGAAAGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTCAGAAGGAATTGTCATAc                         <  1:99513/151‑1 (MQ=255)                                                                                                                             aaaGATATCCTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTCAGAAGGAATTGTCATACAACg                     >  2:6461/1‑151 (MQ=255)                                                                                                                                     ccTGTGCCAGCCAGATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTCAGAAGGAATTGTCATACAACGGTGCTGTa             <  2:531511/151‑1 (MQ=255)                                                                                                                                                cagATGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTCAGAAGGAATTGTCATACAACGGTGCTGTAAAAGGTTGTCg  <  1:441690/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCCACACAGTTACAGCCCGCTTCGATCGCCAGTTCAACAATGTTTTTCGGGTCAAAGTAGAGCGGGTTAGCAGCAAATGAAGCTCCGGCAGAGTGCTCAACGCCCTGGTCAACCGGCAGAATAGAAAGATATCCTGTGCCAGCCAGACGCCCGGTGTTGTACAACGTCTGCATATTACGTAACACCGCTGGCGGGCGATTATTGTCAATCATTACGCGGTCTACGTAGTCATGTCCGGGGAGATAAAGCTGGTCAGAAGGAATTGTCATACAACGGTGCTGTAAAAGGTTGTCG  >  NC_000913/2178236‑2178529

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