Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,455,044	A→G	L189P (CTG→CCG)	fbp ←	fructose‑1,6‑bisphosphatase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,455,044	0	A	G	100.0%	182.3 / NA	51	L189P (CTG→CCG)	fbp	fructose‑1,6‑bisphosphatase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/25);  total (26/25)

GTTGGTGGATTTATCTTCTTCCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAT  >  NC_000913/4454914‑4455176                                                                                                                                   |                                                                                                                                     ggtggtggATTTATCTTCTTCCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGaaaa                                                                                                                             >  2:610193/3‑140 (MQ=255) gTTGGTGGATTTATCTTCTTCCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGaaaa                                                                                                                             >  2:387230/1‑140 (MQ=255) gTTGGTGGATTTATCTTCTTCCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGaaaa                                                                                                                             >  2:497488/1‑140 (MQ=255) gTTGGTGGATTTATCTTCTTCCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGaaaa                                                                                                                             >  1:38508/1‑140 (MQ=255)              tcttctTCCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAgg                                                                                                                <  1:425767/140‑1 (MQ=255)                   ttcCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGaa                                                                                                                  <  1:315081/133‑1 (MQ=255)                   ttcCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGaa                                                                                                                  >  2:315081/1‑133 (MQ=255)                        ggCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCgtaagta                                                                                                      <  1:610193/140‑1 (MQ=255)                                ttAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTg                                                                                              <  1:488289/140‑1 (MQ=255)                                    tGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAaca                                                                                          >  1:473000/1‑140 (MQ=255)                                      tACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACAcc                                                                                        <  1:306082/140‑1 (MQ=255)                                       aCTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCg                                                                                       <  1:294336/140‑1 (MQ=255)                                       aCTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCg                                                                                       >  2:294336/1‑140 (MQ=255)                                         ttcttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCa                                                                                     >  1:80638/1‑140 (MQ=255)                                         ttcttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCa                                                                                     >  2:265604/1‑140 (MQ=255)                                         ttcttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCa                                                                                     >  2:463948/1‑140 (MQ=255)                                         ttcttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCa                                                                                     >  2:152163/1‑140 (MQ=255)                                         ttcttcACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTACTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCa                                                                                     >  2:550564/1‑140 (MQ=255)                                              cACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCgg                                                                                <  1:356042/140‑1 (MQ=255)                                              cACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCgg                                                                                <  2:38508/140‑1 (MQ=255)                                                   cGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGgtggtg                                                                           >  2:335795/1‑140 (MQ=255)                                                     ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCa                                                                                     <  2:315868/128‑1 (MQ=255)                                                     ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCa                                                                                     >  1:315868/1‑128 (MQ=255)                                                     ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggtgg                                                                             <  1:515647/136‑1 (MQ=255)                                                     ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCggtgg                                                                             >  2:515647/1‑136 (MQ=255)                                                     ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCgg                                                                                >  2:215268/1‑133 (MQ=255)                                                     ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCgg                                                                                <  1:215268/133‑1 (MQ=255)                                                     ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTa                                                                         >  1:350654/1‑140 (MQ=255)                                                     ttCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTa                                                                         >  2:160953/1‑140 (MQ=255)                                                        ggAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAAc                                                                      <  1:385156/140‑1 (MQ=255)                                                        ggAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAAc                                                                      <  2:350654/140‑1 (MQ=255)                                                             cTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCa                                                                 <  1:152163/140‑1 (MQ=255)                                                             cTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCa                                                                 <  1:265604/140‑1 (MQ=255)                                                              ttAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCAt                                                                <  1:335795/140‑1 (MQ=255)                                                                aaTGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATgg                                                              >  1:443525/1‑140 (MQ=255)                                                                    tAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTaga                                                          >  2:417895/1‑140 (MQ=255)                                                                     aGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTagag                                                         <  1:417895/140‑1 (MQ=255)                                                                      gTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAgagg                                                        <  1:463948/140‑1 (MQ=255)                                                                        ttCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAg                                                                   <  2:479687/127‑1 (MQ=255)                                                                        ttCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAg                                                                   >  1:479687/1‑127 (MQ=255)                                                                            ttCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGt                                                  >  1:216001/1‑140 (MQ=255)                                                                              cGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGtat                                                <  2:443525/140‑1 (MQ=255)                                                                                        gTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTacc                                      <  2:80638/140‑3 (MQ=255)                                                                                              tttGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCg                                <  1:387230/140‑1 (MQ=255)                                                                                                 gCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCaa                             >  2:320188/1‑140 (MQ=255)                                                                                                       ctcCGGGAAGCGCATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCtgtt                       <  1:160953/140‑1 (MQ=255)                                                                                                                 cgcATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAgg             >  1:587217/1‑140 (MQ=255)                                                                                                                 cgcATCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAgg             >  1:523748/1‑140 (MQ=255)                                                                                                                    aTCCGTTCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTg          >  1:215855/1‑140 (MQ=255)                                                                                                                          tCCTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGaa    <  1:497488/140‑1 (MQ=255)                                                                                                                            cTGGCACGGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAt  <  1:320188/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     GTTGGTGGATTTATCTTCTTCCTGGCAGAATTTAATGTACTTCTTCACCCCGTTCGGAAACTTAATGTAGTTTCCTTCGTTGATGGAGTAGGTTTTGCCTTTCTCCGGGAAGCGCATCCGTTCCTGGCACAGGCAGAAAACGCCGAGCGAAGGATCGTAAGTAAAGGCGTGAACACCGCATCCGGTGGTGTAAACCAGCATGGTAGAGGAGCCGTATACCACGTAACCTGCCGCAACCTGTTTGTTACCAGGCTGGAGGAAAT  >  NC_000913/4454914‑4455176

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