Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,582,067	C→A	7.4%	L351M (CTG→ATG) ‡	PP_4910 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,582,067	0	C	A	7.4%	163.9 / 2.2	55	L351M (CTG→ATG) ‡	PP_4910	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base A (2/2);  total (25/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.73e-01

CTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCTT  >  NC_002947/5581920‑5582215                                                                                                                                                    |                                                                                                                                                     cTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGGCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCAccc                                                                                                                                                      >  2:201664/1‑148 (MQ=255) cTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCAccc                                                                                                                                                      <  1:201664/148‑1 (MQ=255)  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)                                                                                        cTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGcc                                                             >  2:383725/1‑150 (MQ=255)                                                                                          ggAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGc                                                                                                                         >  2:1227831/1‑88 (MQ=255)                                                                                          ggAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGc                                                                                                                         <  1:1227831/88‑1 (MQ=255)                                                                                          ggAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCt                                                           >  2:1166923/1‑150 (MQ=255)                                                                                            aGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGc                                                         <  1:1166923/150‑1 (MQ=255)                                                                                             gCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCa                                                        <  1:347523/150‑1 (MQ=255)                                                                                               cagctaccTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGa                                                      <  1:309787/144‑1 (MQ=255)                                                                                               aCGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGa                                                                                                                                <  2:714975/76‑1 (MQ=255)                                                                                               aCGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGa                                                                                                                                >  1:714975/1‑76 (MQ=255)                                                                                                    gCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATagag                                                 <  2:1274813/150‑1 (MQ=255)                                                                                                    gCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGACCGAg                                                 <  2:1308513/150‑1 (MQ=255)                                                                                                     ccTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAgg                                                <  1:56482/150‑1 (MQ=255)                                                                                                      cTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAg                                                 >  2:556953/1‑148 (MQ=255)                                                                                                      cTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAg                                                 <  1:556953/148‑1 (MQ=255)                                                                                                         gCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGa                                                                                                            <  2:1403983/86‑1 (MQ=255)                                                                                                         gCTCTACGACGGCGAAACGCACACGCTGAACGAAGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGa                                                                                                            >  1:1403983/1‑86 (MQ=255)                                                                                                               cgacgGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAg                                      <  2:1177424/150‑1 (MQ=255)                                                                                                                  cgGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCgg                                   <  1:591494/150‑1 (MQ=255)                                                                                                                   ggCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                      <  2:475636/66‑1 (MQ=255)                                                                                                                   ggCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAg                                                                                                                      >  1:475636/1‑66 (MQ=255)                                                                                                                        aaCGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTTCGCGCCCTGGGCCACGCCCCCCATTGCGAGGTAGCTGACAGCGGCATTAc                             >  1:714006/1‑150 (MQ=255)                                                                                                                         aCGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTAcc                            <  1:383725/150‑1 (MQ=255)                                                                                                                             acacGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGc                        <  1:551861/150‑1 (MQ=255)                                                                                                                               acGCTGAACGACGGCAACACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCt                                                                       >  1:632709/1‑101 (MQ=255)                                                                                                                               acGCTGAACGACGGCAACACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCt                                                                       <  2:632709/101‑1 (MQ=255)                                                                                                                                    gAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCa                 >  2:869542/1‑150 (MQ=255)                                                                                                                                    gAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCa                 >  2:805801/1‑150 (MQ=255)                                                                                                                                          gaacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCa                                                                 >  2:1072362/4‑96 (MQ=255)                                                                                                                                          gaacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCa                                                                 <  1:1072362/93‑1 (MQ=255)                                                                                                                                          gaacagcaCCAGGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAgg           >  1:1292700/4‑150 (MQ=255)                                                                                                                                            gcagcaCCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGcc         <  1:958611/150‑1 (MQ=255)                                                                                                                                                 aCCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCAtt                               <  1:1419988/123‑1 (MQ=255)                                                                                                                                                 aCCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCAtt                               >  2:1419988/1‑123 (MQ=255)                                                                                                                                                   ccTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCtt  >  1:43983/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CTTGCCGTGCAGCTGGGCTCGCAGGAGCAACTGCGTCGCCTGCCACGCACACGCCTGGCCGAGCTGACCGAGCGCTACCGTGATTGCCTGGAGCACGCGGCCTCGCTCTACGACGGCGAAACGCACACGCTGAACGACGGCAGCACCCTGGTGCTGTTCCACAGCCGCGAGTGTGGCGAGGACTACCTGACCAACGCCATTTGCTGTGGCGAACTGTTGCGCGCCCTGGGCCACGCCCTGCAGATCGAGGTAGCTGACAGCGGCATTACCCTGCAACTGCAGCTAGGCCTGGTCTT  >  NC_002947/5581920‑5582215

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: ‑