Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,927,241	(C)6→7	coding (1599/2061 nt)	ptrB ←	oligopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,927,235	1	.	C	97.5%	130.0 / ‑3.1	40	coding (1605/2061 nt)	ptrB	oligopeptidase B
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (18/21);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

ACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCTTCGTACC  >  NC_000913/1927097‑1927365                                                                                                                                            |                                                                                                                                    aCCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGTCGTGGAATAATTTCGGGCGTTGATTTATTGGAACGCCCATCAACATG‑c                                                                                                                                    >  2:286717/1‑140 (MQ=255)   cccACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑ccc                                                                                                                                  >  2:578521/1‑140 (MQ=255)   cccACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑ccc                                                                                                                                  <  1:679654/140‑1 (MQ=255)          ttCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccc‑cccgca                                                                                                                          >  2:816415/1‑136 (MQ=255)          ttCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAACTGCGACTCCCATCCACATG‑cccc‑cccgca                                                                                                                          >  1:273397/1‑136 (MQ=255)            cAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACt                                                                                                                         >  1:63454/1‑140 (MQ=255)             aaaCTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAATAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTc                                                                                                                        >  2:643078/1‑140 (MQ=255)                  cACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑ccc                                                                                                                                  <  2:3194/125‑1 (MQ=255)                  cACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑ccc                                                                                                                                  >  1:3194/1‑125 (MQ=255)                  cACCAGTGGTAAGAGCAATTGATTCGTCAGGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGcc                                                                                                                   <  1:343126/140‑1 (MQ=255)                  cACCAGTGGTAAGAGCAATTGAATCATCATGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGcc                                                                                                                   <  1:341768/140‑1 (MQ=255)                    ccAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCa                                                                                                                 <  1:15820/140‑1 (MQ=255)                    ccAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCa                                                                                                                 >  1:757248/1‑140 (MQ=255)                          gTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCAt                                                                                                           <  2:461388/140‑1 (MQ=255)                            aaGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATaa                                                                                                         <  2:757248/140‑1 (MQ=255)                             agagGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAAc                                                                                                        >  2:759519/1‑140 (MQ=255)                             agagGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAAc                                                                                                        >  1:7806/1‑140 (MQ=255)                                  aaTTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGGGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCCCATAACaaatc                                                                                                   >  2:829464/1‑138 (MQ=255)                                    ttGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGa                                                                                                 >  2:507492/1‑140 (MQ=255)                                         tcatcaAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGgag                                                                                            >  2:750231/1‑140 (MQ=255)                                           atcaAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGagag                                                                                          <  2:639620/140‑1 (MQ=255)                                                 cATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCG‑CCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                                                   <  2:24685/140‑1 (MQ=255)                                                  aTCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                                                   <  2:7806/140‑1 (MQ=255)                                                    cGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGccc                                                                                 <  1:759519/140‑1 (MQ=255)                                                      ttgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAg                                                                               >  1:809914/1‑140 (MQ=255)                                                       tgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGt                                                                              <  1:28184/140‑1 (MQ=255)                                                        gttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGtt                                                                             >  1:172934/1‑140 (MQ=255)                                                                        gcGGTACCTGGGCGATACCGCCGTGGAATAACTCCGGGCCTTGATTAAATGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCa                                                             <  2:255708/139‑1 (MQ=255)                                                                        aCGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCa                                                             <  2:172934/140‑1 (MQ=255)                                                                         cGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAg                                                            >  2:470388/1‑140 (MQ=255)                                                                            tACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCa                                                         <  1:578521/140‑1 (MQ=255)                                                                            tACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCa                                                         <  2:63454/140‑1 (MQ=255)                                                                             aCCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCAt                                                        <  1:816415/140‑1 (MQ=255)                                                                                    cGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCAGAGCCCAGGTTTAACAATGCATCGCAGGCATCAAGAtc                                                 >  1:255708/1‑139 (MQ=255)                                                                                        aaCGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCa                                                         <  1:2313/128‑1 (MQ=255)                                                                                        aaCGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCa                                                         >  2:2313/1‑128 (MQ=255)                                                                                        aaCGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCa                                             <  2:134720/140‑1 (MQ=255)                                                                                             cGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCAATAGCCCATTTTTAACAATGCATCGCAGGCATCAAGATAATCATTaaa                                        >  2:348386/1‑140 (MQ=255)                                                                                                           gggCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTttcttct                          >  1:633563/1‑140 (MQ=255)                                                                                                            ggCGTTGATTAATTGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTttcttctt                         >  2:147273/1‑140 (MQ=255)                                                                                                                       aTTGCAACGCCCATCAAAATGCCCCC‑CCGCACTCCAGCCAATCGCATATCAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTc              <  2:273397/140‑1 (MQ=255)                                                                                                                        ttGCAACGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTcc             <  1:470388/140‑1 (MQ=255)                                                                                                                            aaCGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGGATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCt         <  2:339090/140‑1 (MQ=255)                                                                                                                             aCGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCTTTAAACGCATTGTTATTCTTCTGAAATTTTCCGTCtt        >  1:382160/1‑140 (MQ=255)                                                                                                                             aCGCCCATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCtt        >  1:70991/1‑140 (MQ=255)                                                                                                                                 ccATCAACATGCCCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCTTCGTa    <  1:819414/140‑1 (MQ=255)                                                                                                                                   aTCAACATG‑CCCCACCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCTTCGTAcc  >  2:740790/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                    ACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CCCC‑CCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCTTCGTACC  >  NC_000913/1927097‑1927365

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