Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,587,389	G→A	P242L (CCG→CTG)	ugpQ ←	glycerophosphodiester phosphodiesterase, cytosolic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,587,389	0	G	A	100.0%	102.1 / NA	32	P242L (CCG→CTG)	ugpQ	glycerophosphodiester phosphodiesterase, cytosolic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/15);  total (17/15)

CAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTC  >  NC_000913/3587260‑3587518                                                                                                                                  |                                                                                                                                  cAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCa                                                                                                                          >  1:849641/1‑139 (MQ=255)        aaTCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAAt                                                                                                                   <  1:961182/139‑1 (MQ=255)        aaTCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAAt                                                                                                                   <  1:746953/139‑1 (MQ=255)             ggATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCAt                                                                                                              <  2:849641/139‑1 (MQ=255)              gATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATc                                                                                                             >  1:376005/1‑139 (MQ=255)                   ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAAt                                                                                                                            >  2:282841/1‑119 (MQ=255)                   ttCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAAt                                                                                                                            <  1:282841/119‑1 (MQ=255)                        tGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATg                                                                                                   <  2:376005/139‑1 (MQ=255)                                    atatGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAccc                                                                                       >  1:761479/1‑139 (MQ=255)                                        gctgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgc                                                                                   >  1:172622/1‑139 (MQ=255)                                          tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGc                                                                                                                <  1:505352/108‑1 (MQ=255)                                          tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGc                                                                                                                >  2:505352/1‑108 (MQ=255)                                          tgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAgcgc                                                                                 >  2:588191/1‑139 (MQ=255)                                            ccGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCgcag                                                                               <  2:147087/139‑1 (MQ=255)                                                            gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCcaca                                                                                          <  1:506607/112‑1 (MQ=255)                                                            gCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCcaca                                                                                          >  2:506607/1‑112 (MQ=255)                                                             cGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCt                                                              <  2:172622/139‑1 (MQ=255)                                                                         tCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGa                                                  <  2:761479/139‑1 (MQ=255)                                                                             cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCCACTCTGCTGCGCGCTGGGGTTTATTGACGGt                                              >  1:918892/1‑139 (MQ=255)                                                                             cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                                                <  2:199893/121‑1 (MQ=255)                                                                             cccAGCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTgcgcg                                                                >  1:199893/1‑121 (MQ=255)                                                                                 gCACGCCGCAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATaa                                          >  2:995609/1‑139 (MQ=255)                                                                                       cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                                                 <  1:750279/78‑1 (MQ=255)                                                                                       cgcAGCCTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCa                                                                                                 >  2:750279/1‑78 (MQ=255)                                                                                            ccTGACATTCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCg                               <  2:735695/139‑1 (MQ=255)                                                                                                    tCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg                       >  2:546399/1‑139 (MQ=255)                                                                                                    tCCGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCgg                       <  1:588191/139‑1 (MQ=255)                                                                                                      cGTTGAAAATTTTTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCg                     >  1:360788/1‑139 (MQ=255)                                                                                                      cGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCg                     >  2:140139/1‑139 (MQ=255)                                                                                                      cGTTGAAAATTATTGGGCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCg                     >  1:1013672/1‑139 (MQ=255)                                                                                                                    gggCCGTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCAt       <  2:443196/139‑1 (MQ=255)                                                                                                                         gTAAAGTTCAGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTc  >  2:133062/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                  CAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTC  >  NC_000913/3587260‑3587518

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