Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,746,541	G→C	9.8%	intergenic (+32/+67)	PP_5042 → / ← PP_5043	hypothetical protein/PhoPQ‑activated pathogenicity‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,746,541	0	G	C	9.8%	117.9 / 4.2	41	intergenic (+32/+67)	PP_5042/PP_5043	hypothetical protein/PhoPQ‑activated pathogenicity‑like protein
Reads supporting (aligned to +/- strand):  ref base G (20/17);  new base C (2/2);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

TAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTGA  >  NC_002947/5746395‑5746687                                                                                                                                                   |                                                                                                                                                   tAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgg                                                                                                                                                 >  2:696187/1‑150 (MQ=255)  aaaaGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcggg                                                                                                                                                <  2:1220258/150‑1 (MQ=255)    aaGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcgggcg                                                                                                                                              >  2:47591/1‑150 (MQ=255)         ggCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCAAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                         <  1:630424/150‑1 (MQ=255)              tgCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGAtt                                                                                                                                    <  1:713452/150‑1 (MQ=255)                cACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCt                                                                                                                                  >  2:190905/1‑150 (MQ=255)                cACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCt                                                                                                                                  >  2:1284187/1‑150 (MQ=255)                                  cccGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTgg                                                                                                                <  1:419044/150‑1 (MQ=255)                                       aTACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATg                                                                                                           <  1:190905/150‑1 (MQ=255)                                          ccTGCCGCAACCACCCCAAGGACCGGGCACGTCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAg                                                                                                        <  1:643061/150‑1 (MQ=255)                                           cTGCCGCAACCACCCCAAGGTTGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGc                                                                                                       >  2:489328/1‑150 (MQ=255)                                           cTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGc                                                                                                       >  2:543102/1‑150 (MQ=255)                                           cTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGc                                                                                                       >  2:155540/1‑150 (MQ=255)                                           cTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                        >  2:688989/1‑117 (MQ=255)                                           cTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                        <  1:688989/117‑1 (MQ=255)                                            tGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCa                                                                                                      >  1:1258371/1‑150 (MQ=255)                                             gCCGCAACCACCCCAAGGTTGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCAc                                                                                                     <  1:489328/150‑1 (MQ=255)                                                     caccCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCa                                                                                             <  1:1132495/150‑1 (MQ=255)                                                                   gCCACGGCCGGGGTATTGCGCCTACACGACCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGAGCGGGCGGCGCACGATTCCTGCGCAATTCGCAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTca                                                                               <  1:543102/150‑1 (MQ=255)                                                                             gggTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCgggg                                                                     <  1:47591/150‑1 (MQ=255)                                                                                  ttGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTgg                                                                                                                >  2:1181080/1‑102 (MQ=255)                                                                                  ttGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTgg                                                                                                                <  1:1181080/102‑1 (MQ=255)                                                                                     cgcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCg                                                             <  1:258736/150‑1 (MQ=255)                                                                                      gcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCgg                                                            >  2:353700/1‑150 (MQ=255)                                                                                            aaGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCa                                                      >  2:856588/1‑150 (MQ=255)                                                                                                       ggtTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATc                                           >  1:148474/1‑150 (MQ=255)                                                                                                        gtTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATcc                                          <  1:155540/150‑1 (MQ=255)                                                                                                        gtTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                        <  1:282502/56‑1 (MQ=18)                                                                                                        gtTGGGAGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                        >  2:282502/1‑56 (MQ=17)                                                                                                              cGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtggt                                    <  1:696187/150‑1 (MQ=255)                                                                                                                ttGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCtggtgg                                  >  1:1283705/1‑150 (MQ=255)                                                                                                                    gTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTcc                                                              <  1:724074/118‑1 (MQ=255)                                                                                                                    gTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTcc                                                              >  2:724074/1‑118 (MQ=255)                                                                                                                     tCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCt                             >  1:950993/1‑150 (MQ=255)                                                                                                                       ttGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCa                           <  1:856588/150‑1 (MQ=255)                                                                                                                       ttGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTAGATTCTGTAGCAATGGTGGCCTGGGCATAAAGAACGGAACCGACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCa                           >  1:797970/1‑150 (MQ=255)                                                                                                                               gggCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGt                   >  1:1059626/1‑150 (MQ=255)                                                                                                                                     aaaTCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATg                                                                                                           >  1:247469/1‑56 (MQ=255)                                                                                                                                     aaaTCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATg                                                                                                           <  2:247469/56‑1 (MQ=255)                                                                                                                                      aaTCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGc                                                       >  1:995632/1‑107 (MQ=255)                                                                                                                                      aaTCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGc                                                       <  2:995632/107‑1 (MQ=255)                                                                                                                                         cGAGCGCCGCGCGGGCGGCGCGCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCaaa         >  2:666545/1‑150 (MQ=255)                                                                                                                                              gccgcGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTATGCGGCAAAACGCt    >  1:43429/1‑150 (MQ=255)                                                                                                                                               ccgcGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTg   >  1:548004/1‑150 (MQ=255)                                                                                                                                                cgcgcgGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCCGa  <  1:1284187/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGCAACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGCTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCTGA  >  NC_002947/5746395‑5746687

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