Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	101,802	A→C	10.4%	K208Q (AAG→CAG) ‡	prlC →	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	101,802	0	A	C	10.4%	162.4 / 7.2	56	K208Q (AAG→CAG) ‡	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base C (2/4);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCG  >  NC_002947/101660‑101949                                                                                                                                               |                                                                                                                                                    cAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGcc                                                                                                                                               >  5:102011/1‑149 (MQ=255) cAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGcc                                                                                                                                               <  3:54891/149‑1 (MQ=255)   aGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTg                                                                                                                                             >  3:272952/1‑149 (MQ=255)     cTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                           <  4:272952/149‑2 (MQ=255)     cTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                           <  4:141778/149‑2 (MQ=255)     cTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGa                                                                                                                                           >  4:17828/1‑149 (MQ=255)        aGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGAc                                                                                                                                          <  8:27384/147‑1 (MQ=255)           gAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGcag                                                                                                                                     <  8:159076/149‑3 (MQ=255)            aGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTgg                                                                                                                                    >  8:240364/1‑149 (MQ=255)                    aGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAc                                                                                                                            >  6:130923/1‑149 (MQ=255)                    aGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAc                                                                                                                            <  3:113928/149‑1 (MQ=255)                      ccttttGTCCATTCAGGTGGTGGCCGCCACCACGGCCTGGACCAAGCATGTTACCGACGAATCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAccc                                                                                                                          <  8:275819/146‑1 (MQ=255)                           tctcCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGaa                                                                                                                     >  5:190917/1‑149 (MQ=255)                               cAATCAGCTGCTGGCCGCCACCCAGGCCTGGACCTAGCATTTTACCGACGAAGCTGCGCTGGCCGGCCTGAACGATTCGGCCAAGGCGCAGATGGCCGCTGCAGGCCAGGCCAAGGGACTGGACGGCTGGTTGCTCACACTTCAATTca                                                                                                                 >  8:285693/1‑148 (MQ=255)                                  tCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                              <  5:130923/149‑1 (MQ=255)                                       tgctgGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                           >  8:278962/1‑114 (MQ=255)                                       tgctgGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCCTTGGCCTggg                                                                                                                                           <  7:278962/115‑2 (MQ=255)                                       tgctgGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCATGGACGGCTGGTTGATCACCCTTGAAATCCCCAGcta                                                                                                          >  2:175112/1‑148 (MQ=255)                                         ctgGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGctac                                                                                                         <  3:17828/147‑1 (MQ=255)                                               gCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGt                                                                                                 >  3:92473/1‑149 (MQ=255)                                                 cACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtga                                                                                               >  2:6749/1‑149 (MQ=255)                                                   cccAGGCGTGGACAAAGCATGTTACCGACGACGCTGCGCTGGCCGGCCTGACCGATGCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgatg                                                                                             <  6:190917/149‑1 (MQ=255)                                                    ccAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgatga                                                                                            <  7:40459/149‑1 (MQ=255)                                                    ccAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtgatga                                                                                            <  6:102011/149‑1 (MQ=255)                                                        gCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTa                                                                                        <  2:59951/149‑1 (MQ=255)                                                               ccAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGc                                                                                 >  3:159813/1‑149 (MQ=255)                                                               ccAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGc                                                                                 >  1:137090/1‑149 (MQ=255)                                                                     aTGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACcgc                                                                           >  8:209907/1‑149 (MQ=255)                                                                          aCCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCt                                                                      >  7:225124/1‑149 (MQ=255)                                                                           ccGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCt                                                                      <  7:240364/148‑1 (MQ=255)                                                                            cgacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCACGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGAGCACCCTTGAATTCCCCCGCTACTACGCGGTGAGGACCTACGCCAGCCACCGCGCCCTgc                                                                    >  4:157857/1‑149 (MQ=255)                                                                            cgacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgc                                                                    >  4:17661/1‑149 (MQ=255)                                                                              acgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgcgc                                                                  >  2:158514/1‑149 (MQ=255)                                                                                     tGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGc                                                           >  7:142538/1‑149 (MQ=255)                                                                                      gcgcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGc                                                           >  7:179397/1‑148 (MQ=255)                                                                                       cgcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCta                                                         >  4:112264/1‑149 (MQ=255)                                                                                        gcTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCtat                                                        <  8:179397/149‑1 (MQ=255)                                                                                          tGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGc                                                                                     <  3:280753/118‑1 (MQ=255)                                                                                          tGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGc                                                                                     >  4:280753/1‑118 (MQ=255)                                                                                                  cTGACGGTTTCGGCAAAGGCGCATATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGTTTGATCACCCTTTAATTCCCCAGATAATACGCGTTTATGACCTACGCCAGCGACCGGACCCTGCGCGAAGAGCTATACGCCGCCTa                                              <  1:175112/149‑1 (MQ=255)                                                                                                     aCCGATTCGGCCATGGCGCAGTTGGCCGGTGCAGCCCAGGCCATGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                            <  2:137090/148‑1 (MQ=255)                                                                                                        gATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCAc                                        <  4:159813/149‑1 (MQ=255)                                                                                                          ttCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCAccc                                      <  1:6749/149‑1 (MQ=255)                                                                                                               ccAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGACGCCTACTGCACCCGCGcc                                 >  2:206944/1‑149 (MQ=255)                                                                                                                 aaGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTc                               <  7:275669/149‑1 (MQ=255)                                                                                                                        aGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAg                        >  1:212577/1‑149 (MQ=255)                                                                                                                         gATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAgg                       <  8:225124/149‑1 (MQ=255)                                                                                                                         gATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAgg                       <  2:212577/149‑1 (MQ=255)                                                                                                                            ggCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGc                     <  2:231433/148‑1 (MQ=255)                                                                                                                             gCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGccc                   >  2:240163/1‑149 (MQ=255)                                                                                                                             gCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGcc                    <  8:142538/148‑1 (MQ=255)                                                                                                                             gCCGCTGCAGCCCAGGCAAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGcc                    <  3:17661/148‑1 (MQ=255)                                                                                                                                     cgtccAGGCCCTTGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCg           <  5:291417/146‑1 (MQ=255)                                                                                                                                     aGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCg           <  1:158514/149‑1 (MQ=255)                                                                                                                                         cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCgga         >  6:897/1‑146 (MQ=255)                                                                                                                                              cAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCg  <  7:209907/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    CAAGCTCAGCGAGCTGGGCAGCCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCG  >  NC_002947/101660‑101949

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