Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,793,707	C→G	9.1%	intergenic (+16/+111)	PP_2445 → / ← PP_2446	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,793,707	0	C	G	9.1%	126.5 / 4.8	44	intergenic (+16/+111)	PP_2445/PP_2446	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (11/29);  new base G (2/2);  total (13/31)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGG  >  NC_002947/2793561‑2793842                                                                                                                                                    |                                                                                                                                         tCGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTc                                                                                                                                        >  7:334176/1‑149 (MQ=255)    ccGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTcgcg                                                                                                                                     >  6:120010/1‑149 (MQ=255)     cGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCgcgc                                                                                                                                    <  5:21862/149‑1 (MQ=255)        aGGCTTTGGCGCCAATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCGAGGCGcccc                                                                                                                                 <  6:4049/149‑1 (MQ=255)        aGGCTTTGGCGCCAATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGcccc                                                                                                                                 <  4:198959/149‑1 (MQ=255)           cTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCctct                                                                                                                              <  2:20024/149‑1 (MQ=255)           cTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCctct                                                                                                                              <  2:25911/149‑1 (MQ=255)               ggCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACGTCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGc                                                                                                                           <  2:201516/148‑1 (MQ=255)                gcgcCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCgg                                                                                                                         <  3:306686/149‑1 (MQ=255)                 cgcCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGa                                                                                                                        <  8:324701/149‑1 (MQ=255)                  gcCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGa                                                                                                                        >  5:348559/1‑148 (MQ=255)                  gcCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGa                                                                                                                        >  2:202291/1‑148 (MQ=255)                    ccATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑Caa                                                                                                                    <  5:377008/149‑1 (MQ=255)                    ccATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑Caa                                                                                                                    <  8:21491/149‑1 (MQ=255)                      aTGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACTTCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGc                                                                                                                  <  1:2106/149‑1 (MQ=255)                         tACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCg                                                                                                               <  5:264595/149‑1 (MQ=255)                         ggtacatgggCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGccc                                                                                                                <  1:250415/141‑1 (MQ=255)                            ccAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAGAGGGGCGCGAGGCGCCCCtttttaca                                                                                                                         >  7:138469/1‑133 (MQ=255)                             cAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGc                                                                                                           >  3:302290/1‑149 (MQ=255)                              aaGGGCTATGTGTATGCGGCCATGGCGTTCTCGCCGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGc                                                                                                           <  1:293465/149‑1 (MQ=255)                              aaGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGca                                                                                                          >  2:197634/1‑149 (MQ=255)                               aGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCg                                                                                                                          >  2:334717/1‑133 (MQ=255)                               aGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCg                                                                                                                          <  1:334717/133‑1 (MQ=255)                                gggCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGcaca                                                                                                        <  5:199786/149‑1 (MQ=255)                                 ggCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAg                                                                                                       <  2:48674/149‑1 (MQ=255)                                 ggCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAg                                                                                                       <  1:202291/149‑1 (MQ=255)                                    tATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTa                                                                                                    >  1:166716/1‑149 (MQ=255)                                              cGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACCAAGCCCGCTGCACAGGTATTCGCGTg                                                                                            <  1:197634/148‑1 (MQ=255)                                              cGGCCATGGCGTTCTAGACGGCTATCGTGGTCCTCAACATCCTGGCCCGTCGAGGGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTgg                                                                                           <  2:330247/148‑1 (MQ=255)                                                    tGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGc                                                                                    <  5:346019/149‑1 (MQ=255)                                                         ttCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAAcgc                                                                               <  5:120010/149‑1 (MQ=255)                                                              gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAAt                                                                           <  3:354441/148‑1 (MQ=255)                                                               aCGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCTGA‑CAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATAc                                                                         <  8:62006/149‑1 (MQ=37)                                                               aCGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATAc                                                                         <  4:302290/149‑1 (MQ=255)                                                               aCGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATCTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATAc                                                                         <  1:336937/149‑1 (MQ=255)                                                                   cTATCGAGGTCCTCAACTTCCTGGCGCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                     <  4:68671/149‑1 (MQ=255)                                                                   cTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                     <  6:348559/149‑1 (MQ=255)                                                                                 aaCATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtg                                                                      <  6:202336/134‑1 (MQ=255)                                                                                                       gTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATc                                 >  8:240325/1‑149 (MQ=255)                                                                                                            aaGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCg                                        <  7:280790/137‑1 (MQ=255)                                                                                                            aaGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCg                                        >  8:280790/1‑137 (MQ=255)                                                                                                                                       tAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTTCACAGGTATTCGCGTGGCACTGGCAACGGCAATAACTGTGCAGCGGGCATGTCAGTCAGTGGGCTTCGGGAAATCCTGACGTTAGCTCACCCCCGCATGACGCCgg  >  8:223294/1‑149 (MQ=255)                                                                                                                                             ggggCGCGAGGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATa                                                                          >  4:92890/1‑71 (MQ=37)                                                                                                                                             ggggCGCGAGGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATa                                                                          <  3:92890/71‑1 (MQ=37)                                                                                                                                                    |                                                                                                                                         TCGCCGAAGGCTTTGGCGCCCATGTACCCAAGGGCTATGTGTATGCGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGG‑CTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGA‑CAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGG  >  NC_002947/2793561‑2793842

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