Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,645,118	A→G	V22A (GTC→GCC)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,645,118	0	A	G	100.0%	117.6 / NA	33	V22A (GTC→GCC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/18);  total (15/18)

CCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGC  >  NC_000913/3644991‑3645247                                                                                                                                |                                                                                                                                   ccccAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGga                                                                                                                        >  1:478255/1‑140 (MQ=255)            cGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAAATTTAGaa                                                                                                             >  1:23460/1‑140 (MQ=255)               ccACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAAc                                                                                                          <  2:44395/140‑1 (MQ=255)                    tCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAgg                                                                                                     >  1:443082/1‑140 (MQ=255)                    tCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAgg                                                                                                     <  2:478255/140‑1 (MQ=255)                        ccAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAAtttt                                                                                                                 >  2:162543/1‑124 (MQ=255)                        ccAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAAtttt                                                                                                                 <  1:162543/124‑1 (MQ=255)                        ccAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAgg                                                                                                     <  1:198882/136‑1 (MQ=255)                        ccAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAgg                                                                                                     >  2:198882/1‑136 (MQ=255)                         cAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAAtt                                                                                                <  1:437553/140‑1 (MQ=255)                               cTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑gg                                                                                         <  2:142231/140‑1 (MQ=255)                                   cAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGt                                                                                     <  2:90695/140‑1 (MQ=255)                                                          tGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTcgcg                                                              <  1:270382/140‑1 (MQ=255)                                                            ccccTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCa                                                            >  1:433759/1‑140 (MQ=255)                                                            ccccTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCa                                                            >  2:61610/1‑140 (MQ=255)                                                                    gcTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGt                                                    >  1:471910/1‑140 (MQ=255)                                                                                  cACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTaaa                                                >  2:273749/1‑130 (MQ=255)                                                                                  cACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTaaa                                                >  2:94513/1‑130 (MQ=255)                                                                                  cACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTaaa                                                <  1:94513/130‑1 (MQ=255)                                                                                  cACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTaaa                                                <  1:273749/130‑1 (MQ=255)                                                                                  cACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAgg                                      >  2:399017/1‑140 (MQ=255)                                                                                        ttCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAATCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAgcgc                                <  1:375022/140‑1 (MQ=255)                                                                                         tCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCt                               <  1:97089/140‑1 (MQ=255)                                                                                                   cGTTCAATGCCTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGgcgc                     <  1:61610/140‑1 (MQ=255)                                                                                                            ccTTAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATaa            <  1:375053/140‑1 (MQ=255)                                                                                                               tAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtgt         >  1:90646/1‑140 (MQ=255)                                                                                                               tAGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtgt         <  2:90646/140‑1 (MQ=255)                                                                                                                aGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAg                         >  2:157530/1‑124 (MQ=255)                                                                                                                aGTCACGGCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAg                         <  1:157530/124‑1 (MQ=255)                                                                                                                      ggCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtgtgtg      <  2:425624/136‑1 (MQ=255)                                                                                                                      ggCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtgtgtg      >  1:425624/1‑136 (MQ=255)                                                                                                                      ggCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAt           >  2:144188/1‑131 (MQ=255)                                                                                                                      ggCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAt           <  1:144188/131‑1 (MQ=255)                                                                                                                      ggCTGGAACGGCATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGc  >  2:273533/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                   CCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAA‑GGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGC  >  NC_000913/3644991‑3645247

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