Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,643,899	G→A	G428G (GGC→GGT)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,899	0	G	A	100.0%	119.5 / NA	39	G428G (GGC→GGT)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/20);  total (19/20)

GACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACA  >  NC_000913/3643766‑3644034                                                                                                                                      |                                                                                                                                        gACCGAACTCGTGGCAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCCAAGAACCATCAGCTTTACGCATCTGAACTAc                                                                                                                                    >  1:271679/1‑139 (MQ=255)    cGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                                 <  2:597676/139‑1 (MQ=255)      aaCTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGt                                                                                                                               <  1:725374/139‑1 (MQ=255)      aaCTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGt                                                                                                                               <  1:49520/139‑1 (MQ=255)            tGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                                 <  2:211045/131‑1 (MQ=255)            tGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                                 >  1:211045/1‑131 (MQ=255)                            cTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                                 >  2:306054/1‑115 (MQ=255)                            cTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCa                                                                                                                                 <  1:306054/115‑1 (MQ=255)                               cgtcgtGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTGCGCATCTGACCTACGCAGTCATCCATCCACGCTCCGCCGCGCtt                                                                                                      >  2:3798/1‑139 (MQ=255)                                  cgtGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtt                                                                                                   >  2:432260/1‑139 (MQ=255)                                      agtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                               <  1:583107/139‑1 (MQ=255)                                      agtgaACAGCGCCGGTTTACCATTTACCGGGCGGCTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                               >  1:703964/1‑139 (MQ=255)                                              gcgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTCTTTGCAAAGAACCATCAGCCTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                       >  1:575313/1‑139 (MQ=255)                                               cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACgg                                                                                            >  1:755333/1‑133 (MQ=255)                                               cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACgg                                                                                            <  2:755333/133‑1 (MQ=255)                                                 ccGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAg                                                                                    >  2:361685/1‑139 (MQ=255)                                                  cGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGa                                                                                   >  1:101617/1‑139 (MQ=255)                                                   ggTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                                  <  1:361685/139‑1 (MQ=255)                                                               aCCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                             <  1:554643/116‑1 (MQ=255)                                                               aCCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                             >  2:554643/1‑116 (MQ=255)                                                                   ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                             >  2:139029/1‑112 (MQ=255)                                                                   ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                             <  1:139029/112‑1 (MQ=255)                                                                   ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTAcc                                                                  >  1:331582/1‑139 (MQ=255)                                                                      ggTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcgcg                                                               <  2:101617/139‑1 (MQ=255)                                                                       gTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTcatc                                                                                                                           <  1:291121/78‑1 (MQ=255)                                                                       gTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTcatc                                                                                                                           >  2:291121/1‑78 (MQ=255)                                                                        ttGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCa                                                             <  2:657203/139‑1 (MQ=255)                                                                                        aaTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                               <  2:478119/89‑1 (MQ=255)                                                                                        aaTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCACGCCGCGCTTGTTTTCa                                                                                               >  1:478119/1‑89 (MQ=255)                                                                                         aTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCt                                                                                                                                      >  2:88795/1‑49 (MQ=255)                                                                                         aTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCt                                                                                                                                      <  1:88795/49‑1 (MQ=255)                                                                                              cgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                               >  1:624804/1‑83 (MQ=255)                                                                                              cgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                               <  2:624804/83‑1 (MQ=255)                                                                                              cgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                    <  2:124989/126‑1 (MQ=255)                                                                                              cgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                    >  1:124989/1‑126 (MQ=255)                                                                                                      ttttGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCgaaga                               <  1:432260/139‑1 (MQ=255)                                                                                                         tGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAAc                            >  1:495153/1‑139 (MQ=255)                                                                                                                          tttGCGCATCTGACCTACGCAGTCATCCATCCACGCTCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAg           <  1:3798/139‑1 (MQ=255)                                                                                                                          tttACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAg           <  1:622113/139‑1 (MQ=255)                                                                                                                                   cTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACa  >  2:405983/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                        GACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACA  >  NC_000913/3643766‑3644034

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