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%	73.4 / NA	24	G428G (GGC→GGT)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/12);  total (12/12)

GACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGC  >  NC_000913/3643766‑3644022                                                                                                                                      |                                                                                                                            gACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTAc                                                                                                                        >  1:579291/1‑139 (MQ=255)      aaCTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGt                                                                                                                   >  1:190033/1‑139 (MQ=255)      aaCTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGt                                                                                                                   <  1:863291/139‑1 (MQ=255)                        aTCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTcat                                                                                                                <  1:42135/124‑1 (MQ=255)                        aTCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTcat                                                                                                                >  2:42135/1‑124 (MQ=255)                                      agttaACAGCGCGGTTTCGCCGTTAACAGGGCGGTCGACGTTACAAGTCAATGACGGGCCCGGTTTTTGAAAAGAAACATCACTTTTAGCAATCGGACGTACGCGGTCATCAATCCCCGCCCCGCCGCGCTTGTTTCCa                                                                                   <  2:726407/139‑1 (MQ=255)                                      agtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                   <  2:579291/139‑1 (MQ=255)                                         gaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACgg                                                                                <  1:602057/139‑1 (MQ=255)                                         gaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACgg                                                                                <  2:190033/139‑1 (MQ=255)                                              gcgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                           >  2:140805/1‑139 (MQ=255)                                                gcCGGTTTACCATTTCCCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGATCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACa                                                                         >  2:718258/1‑139 (MQ=255)                                                                   ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAg                                                                  <  1:787583/127‑1 (MQ=255)                                                                   ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAg                                                                  >  2:787583/1‑127 (MQ=255)                                                                   ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTAcc                                                      >  1:58890/1‑139 (MQ=255)                                                                    gCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcg                                                     >  2:550042/1‑139 (MQ=255)                                                                     cGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACcgc                                                    <  1:550042/139‑1 (MQ=255)                                                                                           aCGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                              <  2:488142/139‑1 (MQ=255)                                                                                               gACCGGTTTTTGCACAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCgcgc                                                                                            >  1:652292/1‑73 (MQ=255)                                                                                               gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCgcgc                                                                                            <  2:652292/73‑1 (MQ=255)                                                                                               gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGt                          >  1:212631/1‑139 (MQ=255)                                                                                                  cGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCg                                   >  1:697586/1‑127 (MQ=255)                                                                                                  cGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCg                                   <  2:697586/127‑1 (MQ=255)                                                                                                                     tCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGAt    >  1:113894/1‑139 (MQ=255)                                                                                                                       aGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGc  <  2:212631/139‑1 (MQ=255)                                                                                                                                   cTGACCTACGCAGTCATCCATCCACGCCCCGCCgcgc                                                                                            <  2:178862/37‑1 (MQ=255)                                                                                                                                      |                                                                                                                            GACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGC  >  NC_000913/3643766‑3644022

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