Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,646,189	G→C	A155G (GCA→GGA)	pepP ←	M24B family Xaa‑Pro dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,646,189	0	G	C	100.0%	70.0 / NA	21	A155G (GCA→GGA)	pepP	M24B family Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (11/10);  total (11/10)

CACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATG  >  CP000730/1646058‑1646318                                                                                                                                     |                                                                                                                                  cACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTacaa                                                                                                                          <  2:402234/141‑1 (MQ=255)             cATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatat                                                                                                              >  1:386097/1‑141 (MQ=255)                    tAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCata                                                                                                       <  2:39055/141‑1 (MQ=255)                             gatATGAGCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatc                                                                                              <  2:239364/139‑1 (MQ=255)                              ctATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatca                                                                                             >  2:96982/1‑141 (MQ=255)                                            tACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                                                               <  2:256378/141‑1 (MQ=255)                                                aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCtttt                                                                           <  2:226819/141‑1 (MQ=255)                                                    tATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTaacaaagtcgatgc                                                                       >  2:433991/1‑128 (MQ=255)                                                     aTCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAg                                                                                   <  2:406509/128‑1 (MQ=255)                                                     aTCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAg                                                                                   >  1:406509/1‑128 (MQ=255)                                                      tCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAAtt                                                                     <  1:205270/141‑1 (MQ=255)                                                       cGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                                                                    >  1:187036/1‑141 (MQ=255)                                                              tgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAAt                                                             >  2:111140/1‑141 (MQ=255)                                                                  ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAATAGTTTACATATATTCATATGTTTGATCAACATTATTAGTTGCTTTTTGTATTAAAGCATTTTCg                                                         >  2:128304/1‑141 (MQ=255)                                                                                   tAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtc                                        >  1:196267/1‑141 (MQ=255)                                                                                     gCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtctc                                      >  2:213368/1‑141 (MQ=255)                                                                                           tGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCGTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAttt                                <  1:368253/141‑1 (MQ=255)                                                                                                      atTGCCTTTAATTCTTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGt                      <  1:77369/141‑1 (MQ=255)                                                                                                         gCCTTTAATTC‑TTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTa                  >  1:139575/1‑141 (MQ=255)                                                                                                                  ttatttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAACTGCt          <  2:205762/138‑1 (MQ=255)                                                                                                                          cAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATg  >  2:15443/1‑141 (MQ=255)                                                                                                                                     |                                                                                                                                  CACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTC‑TTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATG  >  CP000730/1646058‑1646318

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑