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	90.0%	94.8 / 2.3	31	A155G (GCA→GGA)	pepP	M24B family Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base C (13/14);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAATAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTAT  >  CP000730/1646046‑1646313                                                                                                                                                |                                                                                                                             cAATAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCttt                                                                                                                          <  1:269405/148‑1 (MQ=255) cAATAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCttt                                                                                                                          <  1:213119/148‑1 (MQ=255)         ttATCACTTGCAACACCATGTGGTAATGCGCCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                                                                  >  2:75650/1‑148 (MQ=255)         ttATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                                                                  >  2:351297/1‑148 (MQ=255)                   cAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatata                                                                                                        <  1:183050/148‑1 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  1:65445/1‑148 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  1:134546/1‑148 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  2:277660/1‑148 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  1:118642/1‑148 (MQ=255)                    aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                       >  2:102149/1‑148 (MQ=255)                        ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCat                                                                                                   <  2:118683/148‑1 (MQ=255)                              ggTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGttt                                                                                             >  1:350482/1‑148 (MQ=255)                                            aTGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTa                                                                               <  2:317028/148‑1 (MQ=255)                                                  aGATGCTACACTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCATTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgct                                                                         >  2:388755/1‑148 (MQ=255)                                                  aGATGCTACAAGCGTATCGAATGAAGGACCATCTGCTCCTAATTCTAGCATTTTGCTTTCTACTATTGACTGTAATTCTTTTCAACCCCTCCCTCCCTTTACAACAGTTAAAATATATTCATCTTTTTCACTAACAattgacgccggt                                                                         >  2:354130/1‑138 (MQ=255)                                                            aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                                                               <  2:350126/148‑1 (MQ=255)                                                            aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                                                               <  2:229991/148‑1 (MQ=255)                                                            aTCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTa                                                               <  2:148388/148‑1 (MQ=255)                                                                             tCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTGAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  2:205908/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  1:84488/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  1:320497/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  2:315027/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCa                                              <  2:345106/148‑1 (MQ=255)                                                                             tCCATCTACTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCATTTTCAGTCATGTCGACCTTTACAACAGTTAAAATATAGTCATATGTTTAATCACCAATATCAGCTGCTTTTTGAAATAAAGCAATTTCGTCAGCa                                              <  2:114177/148‑1 (MQ=255)                                                                                 tCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCtt                                          >  1:332028/1‑148 (MQ=255)                                                                                               tAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtt                            <  2:63356/148‑1 (MQ=255)                                                                                                    ttttACTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATACCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAACATCTTTGACGTCTCTAATTTTATc                       <  1:91349/148‑1 (MQ=255)                                                                                                    ttttACTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATACCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAACATCTTTGACGTCTCTAATTTTATc                       <  1:378593/148‑1 (MQ=255)                                                                                                    ttttACTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATACCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAACATCTTTGACGTCTCTAATTTTATc                       <  1:139934/148‑1 (MQ=255)                                                                                                              tAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTa             >  2:86100/1‑148 (MQ=255)                                                                                                                         ttAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttat  >  1:111642/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                             CAATAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTAT  >  CP000730/1646046‑1646313

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