Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,625,160	G→A	R89H (CGC→CAC)	serB →	phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,625,160	0	G	A	100.0%	143.0 / NA	43	R89H (CGC→CAC)	serB	phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/21);  total (22/21)

GGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAA  >  NC_000913/4625024‑4625289                                                                                                                                         |                                                                                                                                  ggTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCt                                                                                                                                >  1:595464/1‑140 (MQ=255)     gTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGcc                                                                                                                            >  1:112716/1‑140 (MQ=255)      tGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGccc                                                                                                                           >  2:798455/1‑140 (MQ=255)             gATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTgg                                                                                                                              >  2:399498/1‑130 (MQ=255)             gATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTgg                                                                                                                              <  1:399498/130‑1 (MQ=255)             gATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAgc                                                                                                                    >  1:116840/1‑140 (MQ=255)             gATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAgc                                                                                                                    >  1:398278/1‑140 (MQ=255)              aTAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAgcg                                                                                                                   <  2:595464/140‑1 (MQ=255)              aTAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTCGCCCACGAAgcg                                                                                                                   <  2:137855/140‑1 (MQ=255)                        tCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATg                                                                                                         >  2:159609/1‑140 (MQ=255)                            aCCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGc                                                                                                     >  2:203053/1‑140 (MQ=255)                            aCCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGc                                                                                                     <  1:798455/140‑1 (MQ=255)                                aaTACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCg                                                                                                 <  1:159609/140‑1 (MQ=255)                                aaTACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCg                                                                                                 >  1:198521/1‑140 (MQ=255)                                  tACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCt                                                                                               <  2:458088/140‑1 (MQ=255)                                       gagCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGa                                                                                          >  2:491902/1‑140 (MQ=255)                                       gagCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGGGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGa                                                                                          >  2:750282/1‑140 (MQ=255)                                        agCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGaa                                                                                         <  2:116840/140‑1 (MQ=255)                                            aaCTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATc                                                                                     >  1:504357/1‑140 (MQ=255)                                            aaCTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATc                                                                                     >  1:687168/1‑140 (MQ=255)                                                          aTGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCAc                                                                       <  1:164086/140‑1 (MQ=255)                                                                aTTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCggg                                                                 <  1:491902/140‑1 (MQ=255)                                                                      gccgccTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCt                                                           >  1:545349/1‑140 (MQ=255)                                                                       ccgccTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTg                                                          >  2:106529/1‑140 (MQ=255)                                                                       ccgccTGGTGCGTGGAAGATTATCAGGTGACTCGTCTGGCAGGTTCACTCACCGCCCGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCAACTGCGCACGCCGGGTTTGCTg                                                          <  2:603276/140‑1 (MQ=255)                                                                         gccTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGt                                                        >  1:745073/1‑140 (MQ=255)                                                                          ccTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTg                                                       >  1:579839/1‑140 (MQ=255)                                                                           cTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGa                                                      <  2:687168/140‑1 (MQ=255)                                                                           cTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGa                                                      <  2:745073/140‑1 (MQ=255)                                                                                 cGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGata                                                >  1:850422/1‑140 (MQ=255)                                                                                   tGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATg                                              <  2:545349/140‑1 (MQ=255)                                                                                   tGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATg                                              >  1:262192/1‑140 (MQ=255)                                                                                        gATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTc                                         <  2:112716/140‑1 (MQ=255)                                                                                        gATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTc                                         <  1:165826/140‑1 (MQ=255)                                                                                        gATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTc                                         <  1:32124/140‑1 (MQ=255)                                                                                             tCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATg                                                    >  2:44464/1‑124 (MQ=255)                                                                                             tCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATg                                                    <  1:44464/124‑1 (MQ=255)                                                                                                ggTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCa                                 >  2:21444/1‑140 (MQ=255)                                                                                                    aTTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCa                             >  2:795143/1‑140 (MQ=255)                                                                                                      tCGTCTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGa                           <  1:795143/140‑1 (MQ=255)                                                                                                          cTGGCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGa                       <  1:16980/140‑1 (MQ=255)                                                                                                            ggCAGGTTCACTCACCGCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAAt                     <  2:398278/140‑1 (MQ=255)                                                                                                                            gCACGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGc     >  1:485370/1‑140 (MQ=255)                                                                                                                               cGGGCTACACACCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCaa  <  1:573059/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  GGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAA  >  NC_000913/4625024‑4625289

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