Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	102,733	T→C	10.8%	L518P (CTG→CCG)	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	102,733	0	T	C	10.8%	99.1 / 4.8	37	L518P (CTG→CCG)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base T (20/13);  new base C (2/2);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCCTGCACCACATGCTGACCCGCATTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAGTTGGAGTTCTCGCTGTTCGATTTCGAG  >  NC_002947/102604‑102863                                                                                                                                  |                                                                                                                                   gCCTGCACCACATGCTGACCCGCATTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAGCTGGTGCTGGGAGCCGGAAGGCCTGGCACTGAtttg                                                                                                                              >  1:625316/1‑135 (MQ=255)  ccTGCACCACATGCTGACCCGCATTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCt                                                                                                                             <  1:228372/136‑1 (MQ=255)             tGCTGACCCGCATTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGa                                                                                                                  <  2:226381/136‑1 (MQ=255)                     cGCATTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGgcg                                                                                                          >  1:298417/1‑136 (MQ=255)                      gCATTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGgcgc                                                                                                         >  1:275408/1‑136 (MQ=255)                      gCATTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGgcgc                                                                                                         >  1:420583/1‑136 (MQ=255)                        aTTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGcgccg                                                                                                       <  1:463669/136‑1 (MQ=255)                        aTTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGcgccg                                                                                                       <  2:427079/136‑1 (MQ=255)                        aTTGAGCATGCCGGCGGGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGcgccg                                                                                                       <  2:645509/136‑1 (MQ=255)                         ttGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGcgccgc                                                                                                      <  1:564951/136‑1 (MQ=255)                          tGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCgccgcc                                                                                                     >  1:249911/1‑136 (MQ=255)                               aTGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGcc                                                                                                >  2:603992/1‑136 (MQ=255)                                 gCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGcccc                                                                                              <  2:388997/136‑1 (MQ=255)                                    ggCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAgg                                                                                           >  2:554877/1‑136 (MQ=255)                                    ggCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACCGATTTCt                                                                                                                             <  1:344810/102‑1 (MQ=255)                                    ggCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACCGATTTCt                                                                                                                             <  1:341419/102‑1 (MQ=255)                                    ggCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACCGATTTCt                                                                                                                             >  2:344810/1‑102 (MQ=255)                                    ggCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACCGATTTCt                                                                                                                             >  2:341419/1‑102 (MQ=255)                                         gtCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACctg                                                                                      >  1:105441/1‑136 (MQ=255)                                         gtCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACctg                                                                                      >  2:333080/1‑136 (MQ=255)                                                             ccTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGcgccgc                                                                                                      >  2:650221/1‑100 (MQ=255)                                                             ccTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGcgccgc                                                                                                      <  1:650221/100‑1 (MQ=255)                                                                         tCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCa                                                      >  1:534350/1‑136 (MQ=255)                                                                           gAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGt                                                    >  2:480047/1‑136 (MQ=255)                                                                           gAGCTGCCCAGCCAGTTCATGGAAAACTAGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGt                                                    >  2:29209/1‑136 (MQ=255)                                                                            tgCTGCCCATCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTc                                                   <  2:358463/135‑1 (MQ=255)                                                                               tGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCggg                                                >  2:54545/1‑136 (MQ=255)                                                                                    agccagTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCatga                                           >  2:301995/1‑136 (MQ=255)                                                                                    agccagTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCatga                                           >  1:69363/1‑136 (MQ=255)                                                                                      ccagTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCatgatg                                         <  2:418584/136‑1 (MQ=255)                                                                                                  aaaCTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAg                             <  2:420583/136‑1 (MQ=255)                                                                                                           cTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAGTTGGAGTtc                    <  2:380814/136‑1 (MQ=255)                                                                                                           cTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAGTTGGAGTtc                    <  1:603992/136‑1 (MQ=255)                                                                                                              ggAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAGTTGGAGTTCTCg                 >  1:505542/1‑136 (MQ=255)                                                                                                              ggAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAGTTGGAGTTCTCg                 >  1:248221/1‑136 (MQ=255)                                                                                                                        ggCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAGTTGGAGTTCTCGCTGTTCGAtt       >  1:549068/1‑136 (MQ=255)                                                                                                                             ggCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAGTTGGAGTTCTCGCTGTTCGATTTCGAg  >  1:249106/1‑136 (MQ=255)                                                                                                                                  |                                                                                                                                   GCCTGCACCACATGCTGACCCGCATTGAGCATGCCGGCGTGTCCGGCATCAACGGCGTGGCCTGGGATGCGGTCGAGCTGCCCAGCCAGTTCATGGAAAACTGGTGCTGGGAGCCGGAAGGCCTGGCACTGATTTCTGGTCACTACGAAACCGGCGCCGCCCTGCCCCAGGACCTGCTGGACAAGATGCTGGCCGCAAAAAACTTCCAGTCGGGCATGATGATGGTGCGCCAGTTGGAGTTCTCGCTGTTCGATTTCGAG  >  NC_002947/102604‑102863

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