Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,455,282	T→G	D110A (GAC→GCC)	fbp ←	fructose‑1,6‑bisphosphatase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	4,455,282	0	T	G	100.0%	162.4 / NA	47	D110A (GAC→GCC)	fbp	fructose‑1,6‑bisphosphatase I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/26);  total (21/26)

GTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGC  >  NC_000913_3_bme_pgi/4455137‑4455425                                                                                                                                                  |                                                                                                                                                gTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCAt                                                                                                                                             <  1:301161/150‑1 (MQ=255)   aaCCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTCCCGTTACCGGCGTGCCAACACGCGTACCGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa                                                                                                                                           <  1:299565/150‑1 (MQ=255)                cTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTAttt                                                                                                                              >  3:8079/1‑150 (MQ=255)                   tttgttACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGc                                                                                                                           <  4:353185/150‑1 (MQ=255)                   tttgttACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGc                                                                                                                           <  1:411763/150‑1 (MQ=255)                       ttACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgtgt                                                                                                                       >  2:55575/1‑150 (MQ=255)                        tACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGtt                                                                                                                      <  1:175154/150‑1 (MQ=255)                         aCCAGGCTGGAGGAATGCCTGGTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTc                                                                                                                     <  1:334486/150‑1 (MQ=255)                                  gaggaAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTc                                                                                                            >  1:550153/1‑150 (MQ=255)                                        aTCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAAAAAATGGTACCGACAGAGAAGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGGCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgac                                                                                                      <  1:417425/150‑1 (MQ=255)                                          cttcttCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacga                                                                                                    >  3:303464/1‑150 (MQ=255)                                                   ttACCGGCGCGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGTCGGAGACGTTAACCTCGATGTTGGACGCGCCATCCAGGGGGGCCCTCAGCACCACGTATTTTGCGTGCTCACGGCCTACCACGACGCCACTCTCAt                                                                                           <  4:72667/150‑1 (MQ=255)                                                          cGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtcttcttc                                                                                    <  1:27544/150‑1 (MQ=255)                                                            tGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATcttcttctt                                                                                  <  1:103484/150‑1 (MQ=255)                                                              ccAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCa                                                                                >  1:451180/1‑150 (MQ=255)                                                              ccAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCa                                                                                >  1:76272/1‑150 (MQ=255)                                                                  cAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAg                                                                                                        <  4:40748/122‑1 (MQ=255)                                                                  cAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAg                                                                                                        >  3:40748/1‑122 (MQ=255)                                                                   aGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGc                                                                           >  3:6095/1‑150 (MQ=255)                                                                     gCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCaa                                                                         >  1:553802/1‑150 (MQ=255)                                                                      cGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAAt                                                                        <  2:537217/150‑1 (MQ=255)                                                                        tAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGc                                                                      >  4:322250/1‑150 (MQ=255)                                                                        tAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGc                                                                      <  4:272580/150‑1 (MQ=255)                                                                          aCGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGccc                                                                    <  2:131606/150‑1 (MQ=255)                                                                               gcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCGATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAAc                                                               <  4:161433/150‑1 (MQ=255)                                                                                  gTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGat                                                            <  2:137045/150‑1 (MQ=255)                                                                                       tGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATcgc                                                       <  3:259821/150‑1 (MQ=255)                                                                                              aaTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCttt                                                >  2:236075/1‑150 (MQ=255)                                                                                                 ggTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAg                                             <  2:76272/150‑1 (MQ=255)                                                                                                 ggTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGACTTCCAAGAAGAAAATCTCATATTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAg                                             >  1:239347/1‑150 (MQ=255)                                                                                                  gTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGt                                            >  2:563631/1‑150 (MQ=255)                                                                                                      cGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCgg                                        >  1:360670/1‑150 (MQ=255)                                                                                                      cGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCgg                                        <  3:125535/150‑1 (MQ=255)                                                                                                       gACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGc                                       >  1:464695/1‑150 (MQ=255)                                                                                                       gACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGc                                       >  1:215948/1‑150 (MQ=255)                                                                                                          agagACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACCGCCTTCAAAGACGACACTCTCATCTTCTTCTTCCCAGGCAATGCCCGCAACGCTATCGCCTGCCCTCACTGCGGCttc                                    >  2:303446/1‑149 (MQ=255)                                                                                                          agagACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCGACGATATCGCGTGCTTTCAGTGCGGCttt                                    <  1:527366/150‑1 (MQ=255)                                                                                                            agaCGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCa                                  <  4:303464/150‑1 (MQ=255)                                                                                                                gTTAACATCGGTGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACCACAACCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATGTCGCGTGCTTTCAGTGCGGCTTTCAGttt                              <  2:464592/150‑1 (MQ=255)                                                                                                                 ttAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGtttt                             <  2:390079/150‑1 (MQ=255)                                                                                                                        cGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATAGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAg                      >  1:314034/1‑150 (MQ=255)                                                                                                                         gATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGc                     >  4:233026/1‑150 (MQ=255)                                                                                                                             ttGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAAc                 >  2:204447/1‑150 (MQ=255)                                                                                                                                gACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTc                                               <  2:333381/117‑1 (MQ=255)                                                                                                                                gACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTc                                               >  1:333381/1‑117 (MQ=255)                                                                                                                                           cAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTctg   <  4:48835/150‑1 (MQ=255)                                                                                                                                            aGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTctgc  <  2:256086/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                GTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAATCTTCTTCCGTTACCGGCGTGCCAACAGGCGTAACGCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAATGCCCGCAACGATATCGCGTGCTTTCAGTGCGGCTTTCAGTTTTTCATTAGCGAACAAGTCGAGTTTCTGC  >  NC_000913_3_bme_pgi/4455137‑4455425

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