Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,746,485	G→C	8.5%	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,485	0	G	C	8.5%	127.3 / 3.4	47	intergenic (+32/+67)	PP_5042/PP_5043	hypothetical protein/PhoPQ‑activated pathogenicity‑like protein
Reads supporting (aligned to +/- strand):  ref base G (26/17);  new base C (2/2);  total (28/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 = 4.10e-01

TAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCT  >  NC_002947/5746339‑5746629                                                                                                                                                     |                                                                                                                                                  tAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCCGGCA‑ACCACCCCAAGGTCGGGCCACTGCCGGGGTATTGCGCCTACAAGCCAAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGAcgcgc                                                                                                                                                  >  1:88938/1‑149 (MQ=255)   aaaGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCgcggg                                                                                                                                               <  5:76904/149‑1 (MQ=255)          gCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGC‑t                                                                                                                                       <  1:358348/149‑1 (MQ=255)                cACCATTTATGACAAGCGCCCGGATACCTGCC‑GCATACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCAAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGC‑TCGAta                                                                                                                                  >  1:56592/1‑148 (MQ=39)                cACCATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTc                                                                                                                                 >  2:337701/1‑149 (MQ=255)                cACCATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCCCGCGCGGCGC‑TCGa                                                                                                                                    <  2:56592/146‑1 (MQ=255)                cACCATTTATGACAAGCGCCCGGATACCAGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGAGTc                                                                                                                                 >  8:144143/1‑149 (MQ=255)                 aCCATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCt                                                                                                                                >  1:349245/1‑149 (MQ=255)                 aCCATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCt                                                                                                                                >  3:139618/1‑149 (MQ=255)                  ccATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCtg                                                                                                                               >  3:271674/1‑149 (MQ=255)                   cATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCtgt                                                                                                                              >  7:211781/1‑149 (MQ=255)                    aTTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGGTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCtgtg                                                                                                                             >  7:141771/1‑149 (MQ=255)                    aTTTAAGACAAGCGCCCGGATACATGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGAttatctg                                                                                                                             >  2:101319/1‑144 (MQ=255)                     tttATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCATGGTCGGGCCACGGCGGGGGTTTTGCGCTTACATGCCCAAGGTGGTTGGGGTTTGTGTCTTGGGCGGGGGCTTAATTCTAGCGCCGCGCGTGCGGCGC‑TCGATTCTGTgg                                                                                                                            <  2:88938/149‑1 (MQ=255)                          gACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATg                                                                                                                       >  3:96156/1‑149 (MQ=255)                           aaaaGCGCACGGATACCTGCC‑TCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCATACAAGCCCAAGGTGGTTGGGCGTTGATTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGc                                                                                                                      <  4:229298/147‑1 (MQ=255)                            cAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGa                                                                                                                         >  3:232070/1‑145 (MQ=255)                                  cccGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGAtt                                                                                                                                  <  3:59828/130‑1 (MQ=255)                                  cccGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGAtt                                                                                                                                  >  4:59828/1‑130 (MQ=255)                                     ggATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGAATAAATCGACCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCCAGCTCGGc                                                                                                            <  1:337701/149‑1 (MQ=255)                                            cgCC‑CCA‑CCCCCCCCATGGGGGGGCCCCGGGGGGGTTTTTGCGCCCACAACCCCCTGGGTGGTGGGGGTTGTTGCTTGGGCGGGGTCATTAATGCTCGGGCGCGGGGGGGGCGCGTCGATTCTGTGGCGATGCCGAGCTGGGCATGAAg                                                                                                      <  2:349245/148‑1 (MQ=255)                                                   a‑aCCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGaa                                                                                               <  7:241147/149‑1 (MQ=255)                                                   a‑aCCACCCCAAGGTCGGGCAACGGCCGGGGTATTGCGCCAACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGCCATAAATCGAGCGCCGCGGGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGaa                                                                                               <  2:143149/149‑1 (MQ=255)                                                       caccCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCAAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TAGATTCTGTGGCGAAGCCGAGCTGGCCATGAAGCACTGaa                                                                                               >  7:121167/1‑146 (MQ=255)                                                         ccccAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGa                                                                                                                                    <  2:270733/107‑1 (MQ=255)                                                         ccccAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGa                                                                                                                                    >  1:270733/1‑107 (MQ=255)                                                         ccccAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCAc                                                                                          <  4:271674/149‑1 (MQ=255)                                                          cccAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTAGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCAcc                                                                                         >  7:71702/1‑149 (MQ=255)                                                                gTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTcc                                                                                    >  4:93512/1‑148 (MQ=255)                                                                     gCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTAGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTc                                                                              <  8:71702/149‑1 (MQ=255)                                                                      ccACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCAttt                                                                               >  6:318716/1‑147 (MQ=255)                                                                        aCGGCCGGGGAATTGCGCCTACAACCCCAAGGTGGTGGGGCGTTGAGGCTTCCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTcagc                                                                           <  7:144143/149‑1 (MQ=255)                                                                                      gcgcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTc                                                             >  3:113012/1‑149 (MQ=255)                                                                                      gcgcCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTc                                                             >  3:360849/1‑149 (MQ=255)                                                                                          cTACAAGCCCAAGGGGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGc                                                         <  3:51268/149‑1 (MQ=255)                                                                                           tACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGGGGGGAAGGTCCGGc                                                         >  4:195474/1‑148 (MQ=255)                                                                                            aCAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGc                                                         >  1:165018/1‑147 (MQ=255)                                                                                                 cccAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTc                                                             >  3:338987/1‑138 (MQ=255)                                                                                                    aaGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAAACGAGCGCCGCGCGGGCGGCGC‑TCGAATCTCTGGCGATGCCGAGCTGGTCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCAGGCAAGCAccccc                                               >  1:112836/1‑149 (MQ=255)                                                                                                      ggtggtTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCTCGCGGGCGGCTC‑TCGATTCTGTGGCGATGCCGAGCTTGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCt                                             <  6:375174/149‑1 (MQ=255)                                                                                                                      gTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATg                             <  4:96156/149‑1 (MQ=255)                                                                                                                       tCTTGCGCCGGGGCATAAATCGAGGGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATg                             <  8:211781/148‑1 (MQ=255)                                                                                                                         ttGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTc                          <  2:112836/149‑1 (MQ=255)                                                                                                                                        aaTCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGGCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCgg           <  5:318716/149‑1 (MQ=255)                                                                                                                                          tCGAGCGCCGCCCGCGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCGGGTGGATGCTCAGGATAAGTCTGCGGc          >  4:88286/1‑148 (MQ=255)                                                                                                                                           cGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCCGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCa         >  5:283297/1‑148 (MQ=255)                                                                                                                                            gAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGACATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCAGCgttttg       >  8:335876/1‑144 (MQ=255)                                                                                                                                                 ccgcGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCt  >  6:359761/1‑149 (MQ=255)                                                                                                                                                 ccgcGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCt  >  5:11345/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TAAAAGCCGGCTGTGCACCATTTATGACAAGCGCCCGGATACCTGCC‑GCA‑ACCACCCCAAGGTCGGGCCACGGCCGGGGTATTGCGCCTACAAGCCCAAGGTGGTTGGGCGTTGAGTCTTGCGCCGGGGCATAAATCGAGCGCCGCGCGGGCGGCGC‑TCGATTCTGTGGCGATGCCGAGCTGGGCATGAAGCACTGAACCCACCGCTCCCATTTCAGCCAGCGGGGAAGCTCCGGCAAGCACCCCCCTGATCCGCTGGTGGATGCTCAGGATAAGTCTGCGGCAAAACGCT  >  NC_002947/5746339‑5746629

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