Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,793,755	A→C	19.0%	intergenic (+64/+63)	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,755	0	A	C	19.0%	49.5 / 5.2	26	intergenic (+64/+63)	PP_2445/PP_2446	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/13);  new base C (3/2);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGA‑GGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGACGC  >  NC_002947/2793607‑2793902                                                                                                                                                     |                                                                                                                                                     ggCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCa                                                                                                                                                      <  1:295532/149‑1 (MQ=255)              ccGACGGCTATCCAGCGCCTCAACCTCCTGGCCCGGCGAGCGCGTTCGCAGCGTGAGGCGGCTGAAGGATGATGTAAAAAGGCGCGCCGCGCGCCCCTCTTGCGGACAAGCCCGCTGCCCAGGTATTCGCGTGGCACTGGCAACGCGaa                                                                                                                                         <  8:295456/148‑1 (MQ=255)               cGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGaa                                                                                                                                         >  8:206410/1‑148 (MQ=255)                gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAAt                                                                                                                                        <  3:108018/148‑1 (MQ=255)                gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATa                                                                                                                                       <  4:75099/149‑1 (MQ=255)                   ggCTATTGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGTCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCt                                                                                                                                    <  3:172399/149‑1 (MQ=255)                   ggCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATAcca                                                                                                                                    >  5:211336/1‑148 (MQ=38)                   ggCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATACCt                                                                                                                                    <  6:272494/149‑1 (MQ=38)                     cTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAACAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                                                                                  <  5:55577/149‑1 (MQ=255)                     cTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                                                                                  <  5:156223/149‑1 (MQ=255)                     cTATCGAGGTCCTCAACATCCTGGCCCATCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                                                                                  <  8:227219/149‑1 (MQ=255)                        tCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgtg                                                                                                                                 >  7:303059/1‑147 (MQ=255)                                cTCACCACCCTCGCCCGTCGTGCGCGTGTGTAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCtt                                                                                                                       <  4:236396/149‑1 (MQ=38)                                cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAAt                                                                                                                                        <  1:128943/132‑1 (MQ=37)                                cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATa                                                                                                                                       >  2:301992/1‑133 (MQ=255)                                cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATa                                                                                                                                       <  1:301992/133‑1 (MQ=255)                                cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCtt                                                                                                                       >  3:236396/1‑149 (MQ=255)                                   aaCATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAACGTTGATGTAAAAAGGGGCGCCTAGCGCCCCTCTTGAGGACAAGCCAGCTGCACAGGTATTCCCGTGGCACTg                                                                                                                                                   >  7:152397/1‑118 (MQ=255)                                             cccGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAACCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTggg                                                                                                          <  7:206410/149‑1 (MQ=255)                                                                  gTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCAGCCCAAGTATTAGCGTGGCACAGGCAAAGCGAATAACAGTGCCCCGGGCGTGTCAGGCACTGGGCGTAGGGAAATCCTGA‑GGTTa                                                                                    >  2:318813/1‑149 (MQ=255)                                                                                          aaaaGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTAATCGCGAGGCCCTGTCACCGCGAATACCtgtac                                                                                                                                >  7:207725/1‑80 (MQ=37)                                                                                                                    ctgACACGCCCGCTGCACAGGTATTCGCGGGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGA‑GGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCgg                                  <  8:303059/147‑1 (MQ=255)                                                                                                                         aaGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGA‑GGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGcc                             >  3:218259/1‑149 (MQ=25)                                                                                                                              cGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGA‑GGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTTGCTTGGCGGGTGCCTGCgg                        >  8:134286/1‑149 (MQ=255)                                                                                                                                      atgTAGTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAATTCCTGA‑GGTTAGCTCACAGCGGCATTCCGCGGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGt                                <  5:198108/131‑1 (MQ=255)                                                                                                                                              gcgTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGA‑GGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAggtg               >  4:102536/1‑142 (MQ=255)                                                                                                                                              gcgTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGA‑GGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAggtg               <  3:102536/142‑1 (MQ=255)                                                                                                                                                     aCTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGGTTAGCTAACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGAcgc  >  1:256588/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGA‑GGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGACGC  >  NC_002947/2793607‑2793902

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