Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,202,981	C→A	7.0%	L430M (CTG→ATG)	PP_4583 →	putative Peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,202,981	0	C	A	7.0%	182.1 / 4.1	57	L430M (CTG→ATG)	PP_4583	putative Peptidase
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base A (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AAGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTT  >  NC_002947/5202835‑5203127                                                                                                                                                   |                                                                                                                                                   aaGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTg                                                                                                                                                  <  1:101366/149‑1 (MQ=255) aaGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTg                                                                                                                                                  >  3:194506/1‑149 (MQ=255)     gcTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGc                                                                                                                                              >  5:245252/1‑149 (MQ=255)       tCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGcc                                                                                                                                             >  5:15112/1‑148 (MQ=255)           ccGCCCGGCCAAGGGGCGACAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCACGCAACCCCCGTATCTCGGTGCATTCCATCCCGATGTCTATGTCAGTGAGCGCCTGTGGGCAACCGCTGCGGACGGCACCAGGGTGcggtccg                                                                                                                                        <  8:268475/149‑7 (MQ=255)             gCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGc                                                                                                                                      >  4:72602/1‑149 (MQ=255)             gCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGc                                                                                                                                      >  5:167145/1‑149 (MQ=255)                              cAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGTAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCa                                                                                                                     <  1:77389/149‑1 (MQ=255)                               aaCTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAg                                                                                                                    <  7:329453/149‑1 (MQ=255)                                           ccACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGc                                                                                                           >  4:244978/1‑146 (MQ=255)                                                  cgcCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGc                                                                                                 <  6:245252/149‑1 (MQ=255)                                                    cccAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCa                                                                                               <  7:72084/149‑1 (MQ=255)                                                       aGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActg                                                                                            <  4:62146/149‑1 (MQ=255)                                                          tggtgCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtac                                                                                         <  6:15112/149‑1 (MQ=255)                                                                  aaGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCActgtacctgtac                                                                                   <  4:194506/147‑1 (MQ=255)                                                                     cAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCt                                                                                                                                    >  2:121685/1‑95 (MQ=255)                                                                      aaaCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGTCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTAt                                                                             >  7:39074/1‑149 (MQ=255)                                                                      aaaCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTAt                                                                             >  5:76039/1‑149 (MQ=255)                                                                       aaCCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATg                                                                            >  2:24234/1‑149 (MQ=255)                                                                           ccGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTTGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGc                                                                        <  5:320187/149‑1 (MQ=255)                                                                            cGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGcc                                                                       <  2:177795/149‑1 (MQ=255)                                                                                 cTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATgg                                                                  <  6:167145/149‑1 (MQ=255)                                                                                    ggTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgccgcc                                                                                                                         <  1:49033/91‑1 (MQ=255)                                                                                    ggTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTAcgccgcc                                                                                                                         >  2:49033/1‑91 (MQ=255)                                                                                     gTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTAt                                                                             <  4:105735/134‑1 (MQ=255)                                                                                     gTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTAt                                                                             >  3:105735/1‑134 (MQ=255)                                                                                          ttCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCt                                                         >  4:67132/1‑149 (MQ=255)                                                                                          ttCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCt                                                         >  8:227145/1‑149 (MQ=255)                                                                                           tCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCa                                                                                               <  6:192499/110‑1 (MQ=255)                                                                                              aTGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAc                                                     >  3:70127/1‑149 (MQ=255)                                                                                                gCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAccc                                                   >  5:220852/1‑149 (MQ=255)                                                                                                 ccGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAccc                                                   <  4:70127/148‑1 (MQ=255)                                                                                                  cGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCt                                                 >  7:345664/1‑149 (MQ=255)                                                                                                  cGAGGACTATGGCTGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGAGCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGGGCGACTGTCCCTGTACGGCTATGGGGCCTATGGCGTAAGCCTCGACCCCt                                                 <  3:101029/149‑1 (MQ=255)                                                                                                    aTGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGTCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTgg                                               <  8:39074/149‑1 (MQ=255)                                                                                                       aCTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGtt                                             <  3:72602/148‑1 (MQ=255)                                                                                                        cTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTtct                                           <  1:181640/149‑1 (MQ=255)                                                                                                         tATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTtctc                                          <  7:302314/149‑1 (MQ=255)                                                                                                               aGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGc                                                                 <  6:360079/120‑1 (MQ=255)                                                                                                                   aGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATgcgcgc                                >  3:233485/1‑149 (MQ=255)                                                                                                                     cgcCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCt                              <  6:76039/149‑1 (MQ=255)                                                                                                                       ccTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGa                            <  5:42134/149‑1 (MQ=255)                                                                                                                       ccTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGa                            <  8:345664/149‑1 (MQ=255)                                                                                                                        cTGTGGGCAACCGCAGCGGACGGCCCCCTGGTGCCGATCAGCCTGGTACGCCGACCCCACGACCTGGGCAACACCGTGCCACTGTACCAGTACGGCTATGGTGCCTAAGGCGAAAGCCTCGACCCCAGGTTCACACAAGCGCGCCagag                           >  6:62432/1‑149 (MQ=255)                                                                                                                        cTGTGGGCAACCGCAGCGGACGGCACCCTAGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAg                           >  6:318337/1‑149 (MQ=255)                                                                                                                                  ccGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAg                 >  2:133786/1‑149 (MQ=255)                                                                                                                                    gCAGCGGACGGCACCCTGGTGCCGATCAGCCAGGTACGCCGCCGCCAGGACCTGGGAAAGACCCTGCCACTCTACCTGTACGGCTATGCTGCCTATGCCGAAAGCCTAGACCCCAGGTTATCACATGCGCGCCTGAGCCTCCTGGAgcg               >  8:341764/1‑149 (MQ=255)                                                                                                                                      aGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAc                                                     <  7:17815/109‑1 (MQ=255)                                                                                                                                      aGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGAc                                                     >  8:17815/1‑109 (MQ=255)                                                                                                                                        tgatcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCgcggc           >  1:90601/5‑149 (MQ=255)                                                                                                                                          gACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGt         >  5:120189/1‑149 (MQ=255)                                                                                                                                          atcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCctgctg                    <  3:310864/136‑1 (MQ=255)                                                                                                                                          atcGGCACCAGGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCctgctg                    >  4:310864/3‑138 (MQ=255)                                                                                                                                           aCGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTg        <  6:220852/149‑1 (MQ=255)                                                                                                                                               cACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCt    >  4:297717/1‑149 (MQ=255)                                                                                                                                               cACCCTGGTGCCGATCAGCCTGGTACGCCCCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCt    >  5:112629/1‑149 (MQ=255)                                                                                                                                                aCCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCtt   >  8:306700/1‑149 (MQ=255)                                                                                                                                                 cccTGGTGCCGATCAGCCTGGTACGCCCCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCttt  <  6:112629/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AAGCGCTCAACCGCCCGGCCCAGGTGCGCCAACTGACGTTGCCCACGGGCGCCCAGGTGGTGCTCAAGCAAACCCCGGTGCTCGGTGCATTCGATGCCGATGACTATGTCAGTGAGCGCCTGTGGGCAACCGCAGCGGACGGCACCCTGGTGCCGATCAGCCTGGTACGCCGCCGCCAGGACCTGGGCAAGACCGTGCCACTGTACCTGTACGGCTATGGTGCCTATGGCGAAAGCCTCGACCCCTGGTTCTCACATGCGCGCCTGAGCCTGCTGGAGCGCGGCGTGGCCTTT  >  NC_002947/5202835‑5203127

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