Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,645,004	T→C	D60G (GAC→GGC)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,645,004	0	T	C	100.0%	137.7 / NA	39	D60G (GAC→GGC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/17);  total (22/17)

GTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGA  >  NC_000913/3644890‑3645132                                                                                                                   |                                                                                                                                 gTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCt                                                                                                                >  2:458571/1‑133 (MQ=255) gTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACCCCCGAACACATCGCCCGCTTCCGCCAGCGGCt                                                                                                                >  2:206872/1‑133 (MQ=255)   gCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCAAACACATAGCCAACATCCGCCAGCGGCTgg                                                                                                              >  1:256618/1‑133 (MQ=255)      gTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCag                                                                                                           >  2:140328/1‑133 (MQ=255)      gTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCag                                                                                                           >  1:155492/1‑133 (MQ=255)       tATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCaga                                                                                                          >  1:84036/1‑133 (MQ=255)        aTTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCagag                                                                                                         <  1:458571/133‑1 (MQ=255)            cGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCACCAACACAACGCCCGCGACCGCCAGTGGCTTGCAGAGAttt                                                                                                     >  2:298855/1‑133 (MQ=255)               cagcagCGGCGGGGTTTGTTCGTATGCTACACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTcc                                                                                                  <  2:215015/133‑1 (MQ=255)               cagcagCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTcc                                                                                                  <  1:557016/133‑1 (MQ=255)               cagcagCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTcc                                                                                                  <  2:155492/133‑1 (MQ=255)               cagcagCGGCAGGGTATGGTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTcc                                                                                                  <  1:239886/133‑1 (MQ=255)                agcagcGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTccc                                                                                                 <  1:141289/133‑1 (MQ=255)                   agcGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAgg                                                                                              >  1:314026/1‑133 (MQ=255)                    gcGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGgt                                                                                             <  1:316310/133‑1 (MQ=255)                                gTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGcccc                                                                                 <  1:423131/133‑1 (MQ=255)                                gTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGcccc                                                                                 >  2:423131/1‑133 (MQ=255)                                   cGTAGGCTTCTCGCACTTCCGGGCCATTCTTCCCCGAGTTCCGGTGGCTGACCGGGGAGAAGATACGCCCAAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTg                                                                              <  2:163540/133‑1 (MQ=255)                                   cGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTg                                                                              <  2:84036/133‑1 (MQ=255)                                     tAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTgcg                                                                            >  1:80345/1‑133 (MQ=255)                                     tAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTgcg                                                                            <  1:140328/133‑1 (MQ=255)                                           tCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGctacta                                                                      <  1:266097/133‑1 (MQ=255)                                            cACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCtactac                                                                     >  2:377073/1‑133 (MQ=255)                                                 aGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCt                                                                <  2:525964/133‑1 (MQ=255)                                                     ccGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCAc                                                            >  1:456221/1‑133 (MQ=255)                                                             tttttCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCGCTTGCGCCACTACGCGCTCCACGTTTTcgc                                                    >  1:240340/1‑133 (MQ=255)                                                             tttttCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTcgc                                                    >  2:153826/1‑133 (MQ=255)                                                               tttCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCgg                                                  >  1:104241/1‑133 (MQ=255)                                                                   aCCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGt                                              <  2:311968/133‑1 (MQ=255)                                                                       aGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATTGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACCCGCTCCACGTTTTCGCGGCAGTCGtt                                          >  1:474961/1‑133 (MQ=255)                                                                       aGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCGCCTTGCGCTACTTCGCCCTCCACGTATTCGAGGCAGTCGtt                                          >  2:253196/1‑133 (MQ=255)                                                                                  cTGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAg                               >  1:65902/1‑133 (MQ=255)                                                                                   tGACCGGGGAGAAGATACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGt                              >  2:102868/1‑133 (MQ=255)                                                                                             gaagaTACGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGa                    <  2:65902/133‑1 (MQ=255)                                                                                                   aCGCCCCAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACa              <  2:314026/133‑1 (MQ=255)                                                                                                        ccAACACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTc         >  1:314383/1‑133 (MQ=255)                                                                                                          aaCACATCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCg       >  1:62659/1‑133 (MQ=255)                                                                                                               aTCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGgaga  >  1:148410/1‑133 (MQ=255)                                                                                                               aTCGCCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGgaga  <  2:80345/133‑1 (MQ=255)                                                                                                                   |                                                                                                                                 GTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTAGGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGA  >  NC_000913/3644890‑3645132

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