Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,582,068	T→G	7.1%	L351R (CTG→CGG) ‡	PP_4910 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,582,068	0	T	G	7.1%	163.4 / 3.0	56	L351R (CTG→CGG) ‡	PP_4910	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base G (2/2);  total (26/30)
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

TTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCTTGG  >  NC_002947/5581921‑5582217                                                                                                                                                    |                                                                                                                                                       ttGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTgg                                                                                                                                                      >  1:1274813/1‑150 (MQ=255)        gCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCAGGGTGCTgcc                                                                                                                                               <  1:413121/150‑3 (MQ=255)                          gCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTgg                                                                                                                             <  1:278147/150‑1 (MQ=255)                              cTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                         <  2:963047/150‑1 (MQ=255)                              cTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGATGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                         <  1:549602/150‑1 (MQ=255)                                gCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGa                                                                                                                       <  1:867408/150‑1 (MQ=255)                                      ccTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCt                                                                                                                 >  1:1177424/1‑150 (MQ=255)                                       cTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCTTGGTGCTGTTGCACAGCCGCGAGTGTGGAGAGGAGTACCTg                                                                                                                >  2:589566/1‑150 (MQ=255)                                            aCGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCaa                                                                                                           >  2:585684/1‑150 (MQ=255)                                             cGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAAc                                                                                                          <  2:586229/150‑1 (MQ=255)                                               cacaCGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGc                                                                                                        <  1:421747/150‑1 (MQ=255)                                                  acGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCAt                                                                                                     >  2:347523/1‑150 (MQ=255)                                                      cTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGc                                                                                                 >  2:551861/1‑150 (MQ=255)                                                      cTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGc                                                                                                 <  1:585684/150‑1 (MQ=255)                                                          ccGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCtgtg                                                                                             >  2:958611/1‑150 (MQ=255)                                                             aGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCg                                                                                          >  2:56482/1‑150 (MQ=255)                                                                 gACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAATGGTGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACt                                                                                      <  1:1314714/150‑1 (MQ=255)                                                                    cGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGaa                                                                                        <  2:1042479/145‑1 (MQ=255)                                                                    cGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGaa                                                                                        >  1:1042479/1‑145 (MQ=255)                                                                    cGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCAAAACGCACACGCAGACCGACGCCAAGCCTGTTGTCCTGTGCCAGAGGCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGtt                                                                                   <  1:883875/150‑1 (MQ=255)                                                                           tACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGcc                                                                            >  1:922953/1‑150 (MQ=255)                                                                            aCCGTGATTGCATGGAGAACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGccc                                                                           >  2:1059677/1‑150 (MQ=255)                                                                                       cTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGc‑c                                                               >  2:383725/1‑150 (MQ=255)                                                                                         ggAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGc                                                                                                                            <  1:1227831/88‑1 (MQ=255)                                                                                         ggAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGc                                                                                                                            >  2:1227831/1‑88 (MQ=255)                                                                                         ggAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCt                                                             >  2:1166923/1‑150 (MQ=255)                                                                                           aGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGc                                                           <  1:1166923/150‑1 (MQ=255)                                                                                            gCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCa                                                          <  1:347523/150‑1 (MQ=255)                                                                                              cagctaccTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGa                                                        <  1:309787/144‑1 (MQ=255)                                                                                              aCGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGa                                                                                                                                   <  2:714975/76‑1 (MQ=255)                                                                                              aCGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGa                                                                                                                                   >  1:714975/1‑76 (MQ=255)                                                                                                   gCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATagag                                                   <  2:1274813/150‑1 (MQ=255)                                                                                                   gCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGACCGAg                                                   <  2:1308513/150‑1 (MQ=255)                                                                                                    ccTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAgg                                                  <  1:56482/150‑1 (MQ=255)                                                                                                     cTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAg                                                   <  1:556953/148‑1 (MQ=255)                                                                                                     cTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAg                                                   >  2:556953/1‑148 (MQ=255)                                                                                                        gCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGa                                                                                                               <  2:1403983/86‑1 (MQ=255)                                                                                                        gCTCTACGACGGCGAAACGCACACGCTGAACGAAGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGa                                                                                                               >  1:1403983/1‑86 (MQ=255)                                                                                                              cgacgGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAg                                        <  2:1177424/150‑1 (MQ=255)                                                                                                                 cgGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCgg                                     <  1:591494/150‑1 (MQ=255)                                                                                                                  ggCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                         >  1:475636/1‑66 (MQ=255)                                                                                                                  ggCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                         <  2:475636/66‑1 (MQ=255)                                                                                                                       aaCGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTTCGCGCCCTGGGCCACGC‑CCCCCATTGCGAGGTAGCTGACAGCGGCATTAc                               >  1:714006/1‑150 (MQ=255)                                                                                                                        aCGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTAcc                              <  1:383725/150‑1 (MQ=255)                                                                                                                            acacGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGc                          <  1:551861/150‑1 (MQ=255)                                                                                                                              acGCTGAACGACGGCAACACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCt                                                                          >  1:632709/1‑101 (MQ=255)                                                                                                                              acGCTGAACGACGGCAACACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCt                                                                          <  2:632709/101‑1 (MQ=255)                                                                                                                                   gAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCa                   >  2:805801/1‑150 (MQ=255)                                                                                                                                   gAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCa                   >  2:869542/1‑150 (MQ=255)                                                                                                                                         gaacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCa                                                                    >  2:1072362/4‑96 (MQ=255)                                                                                                                                         gaacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCa                                                                    <  1:1072362/93‑1 (MQ=255)                                                                                                                                         gaacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAgg             >  1:1292700/4‑150 (MQ=255)                                                                                                                                           gcagcaCCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGcc           <  1:958611/150‑1 (MQ=255)                                                                                                                                                aCCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCAtt                                 <  1:1419988/123‑1 (MQ=255)                                                                                                                                                aCCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCAtt                                 >  2:1419988/1‑123 (MQ=255)                                                                                                                                                  ccTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCtt    >  1:43983/1‑150 (MQ=255)                                                                                                                                                    tGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCGTTCGCCACGCTCCTG‑ACATTGTTGGAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCTTgg  <  1:885220/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       TTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGC‑CCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCTTGG  >  NC_002947/5581921‑5582217

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