Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,605,418	A→C	7.2%	Q263P (CAG→CCG)	PP_4077 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,605,418	0	A	C	7.2%	148.9 / 2.7	54	Q263P (CAG→CCG)	PP_4077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base C (2/2);  total (26/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.88e-01

AGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGT  >  NC_002947/4605273‑4605565                                                                                                                                                  |                                                                                                                                                    aGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCAAAAGCCTGCCCACAGGTg                                                                                                                                                 >  1:316455/1‑150 (MQ=255) aGCCCCGAGGAGGCCTGGCCGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGCAATACTTCGTGTGCCCGACGCCCTTCCACTTCTTCAGCATTACCGGCATGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTg                                                                                                                                                 <  2:127982/150‑1 (MQ=255)  tcccTGAGGAGGCCTTGCTGCCCTAACCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTCCCGGCCTGGAGCATCTGTGGCCCGCACAAGCCCGCCCACAGGTgg                                                                                                                                                <  2:341115/149‑1 (MQ=255)  gCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTgg                                                                                                                                                <  2:1047378/150‑1 (MQ=255)  gCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTgg                                                                                                                                                <  1:1278388/150‑1 (MQ=255)             gCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTGCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTc                                                                                                                                     <  2:1067764/150‑1 (MQ=255)               cTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCa                                                                                                                                   <  1:15181/150‑1 (MQ=255)                  gctgcCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCAttt                                                                                                                                <  2:331345/150‑1 (MQ=255)                          aCCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTc                                                                                                                        <  2:879598/150‑1 (MQ=255)                              gCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACt                                                                                                                    <  1:746414/150‑1 (MQ=255)                                  aaCGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCa                                                                                                                <  2:590261/150‑1 (MQ=255)                                                  aCCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTg                                                                                                >  1:703905/1‑150 (MQ=255)                                                   ccGTATTCTGCAGGAATACGTCGTGTGCCCGAAGCGCTTCCACTTCTTCAGGATTACCGGCCTGGAAAAGCTGTGGCCGACCCAGGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGttag                                                                                               <  2:306394/150‑3 (MQ=255)                                                            gCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACagaaga                                                                                      >  2:869211/1‑150 (MQ=255)                                                                       tCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGcc                                                                                                                 >  1:1227987/1‑112 (MQ=255)                                                                       tCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGcc                                                                                                                 <  2:1227987/112‑1 (MQ=255)                                                                            ttCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGc                                                                      >  2:1171289/1‑150 (MQ=255)                                                                               ccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTc                                                                                  <  1:739795/135‑1 (MQ=255)                                                                               ccGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTc                                                                                  >  2:739795/1‑135 (MQ=255)                                                                                     cgcTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCt                                                             <  2:703905/150‑1 (MQ=255)                                                                                       cTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGt                                                           >  2:900506/1‑150 (MQ=255)                                                                                            aCTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTTGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAcc                                                      >  1:1057687/1‑150 (MQ=255)                                                                                            aCTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAcc                                                      >  1:886530/1‑150 (MQ=255)                                                                                             cttcttCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCt                                                     <  2:6659/150‑1 (MQ=255)                                                                                             cttcttCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCt                                                     >  1:190872/1‑150 (MQ=255)                                                                                               tcttcAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGt                                                   <  2:886530/150‑1 (MQ=255)                                                                                                 ttcAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTg                                                                                                >  1:1136491/1‑103 (MQ=255)                                                                                                 ttcAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTg                                                                                                <  2:1136491/103‑1 (MQ=255)                                                                                                  tcAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCa                                                >  2:311566/1‑150 (MQ=255)                                                                                                    aGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAg                                              >  2:249282/1‑150 (MQ=255)                                                                                                        ttACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGcaca                                          >  2:1123127/1‑150 (MQ=255)                                                                                                          aCCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGt                                        <  1:869211/150‑1 (MQ=255)                                                                                                                cTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCTGCTGTCAACCTGTTCAAGCACAGTGCCGAg                                  >  1:841319/1‑150 (MQ=255)                                                                                                                cTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAg                                  <  1:249282/150‑1 (MQ=255)                                                                                                                        gCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGa                          >  1:1201333/1‑150 (MQ=255)                                                                                                                         cTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGcaca                                          >  2:861914/1‑133 (MQ=255)                                                                                                                         cTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGcaca                                          <  1:861914/133‑1 (MQ=255)                                                                                                                           gtgGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCt                       >  1:524155/1‑150 (MQ=255)                                                                                                                            tgGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTg                      <  2:1201333/150‑1 (MQ=255)                                                                                                                            tgGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTg                      <  1:1171289/150‑1 (MQ=255)                                                                                                                            tgGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTg                      <  2:116831/150‑1 (MQ=255)                                                                                                                                cGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAAAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCa                              >  2:704943/1‑138 (MQ=255)                                                                                                                                cGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAAAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCa                              <  1:704943/138‑1 (MQ=255)                                                                                                                                      aaGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg            >  2:271782/1‑150 (MQ=255)                                                                                                                                      aaGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg            >  2:294587/1‑150 (MQ=255)                                                                                                                                      aaGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAgg            >  1:540149/1‑150 (MQ=255)                                                                                                                                       tcgatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                       <  2:56853/102‑1 (MQ=255)                                                                                                                                       tcgatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAAc                                                       >  1:56853/5‑106 (MQ=255)                                                                                                                                        gCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCAAGATACCCTTCGCGTTTGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCt          >  1:24509/1‑150 (MQ=255)                                                                                                                                         gatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAg             >  1:523384/3‑146 (MQ=255)                                                                                                                                         gatGCCCACCTGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAg             <  2:523384/144‑1 (MQ=255)                                                                                                                                          cTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTg         <  1:1127119/149‑1 (MQ=255)                                                                                                                                          cTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTg         >  2:1127119/1‑149 (MQ=255)                                                                                                                                              ccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAg    <  2:316455/150‑1 (MQ=255)                                                                                                                                               cacaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTc                                                          >  1:925948/1‑95 (MQ=255)                                                                                                                                               cacaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTc                                                          <  2:925948/95‑1 (MQ=255)                                                                                                                                                acaGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGt  <  1:271782/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AGCCCTGAGGAGGCCTTGCTGCCCTACCCGCAGAACGTGTTCGATGGCTACCGTATTCTGCAGGAATACTTCGTGTTCCCGAAGCGCTTCCACTTCTTCAGCATTACCGGCCTGGAAAAGCTGTGGCCGGCACAAGCCTGCCCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGT  >  NC_002947/4605273‑4605565

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