Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,266,116	G→C	17.0%	intergenic (+33/+26)	PP_2863 → / ← PP_2864	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,266,116	0	G	C	17.0%	100.4 / 12.1	41	intergenic (+33/+26)	PP_2863/PP_2864	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/18);  new base C (2/5);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

AAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACCCA  >  NC_002947/3265973‑3266264                                                                                                                                                |                                                                                                                                                       aaGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCcgc                                                                                                                                                    >  3:251782/1‑148 (MQ=255)     cccGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGCTTTACCCGCGaaga                                                                                                                                               <  8:186815/149‑3 (MQ=255)        gCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCCTTTCGCGGCTTTACCCGCGaagaggc                                                                                                                                            <  1:14652/149‑6 (MQ=255)          aaGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCg                                                                                                                                          <  8:31528/149‑1 (MQ=38)                     ccAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGa                                                                                                                               >  6:171283/1‑149 (MQ=255)                      cAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGAt                                                                                                                              >  4:61379/1‑149 (MQ=255)                            gggtgCTGGGGGATTGGGCTGCTGCACCAGGTGTTCTGCGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACAt                                                                                                                         <  8:121263/146‑1 (MQ=255)                             tgtgCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGa                                                                                                                       <  4:251782/149‑1 (MQ=255)                             tgtgCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGa                                                                                                                       <  7:177656/149‑1 (MQ=255)                               tgCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAc                                                                                                                      <  5:171283/148‑1 (MQ=255)                                gCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAcc                                                                                                                     >  8:283047/1‑148 (MQ=255)                                   gggggATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAAcc                                                                                                                 >  7:299468/1‑149 (MQ=255)                                    ggggATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAc                                                                                                                      >  8:106302/1‑143 (MQ=255)                                    ggggATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGAc                                                                                                                      <  7:106302/143‑1 (MQ=255)                                              ctgctgcTCCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCaa                                                                                                       <  2:143550/148‑1 (MQ=255)                                                    cACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCtt                                                                                                 <  7:283047/148‑1 (MQ=255)                                                     aCCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTgg                                                                                               <  8:299468/149‑1 (MQ=255)                                                          gtgtTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAAc                                                                                          >  3:22479/1‑149 (MQ=255)                                                                tGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATc                                                                                                                             >  8:33364/1‑108 (MQ=255)                                                                  ggCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCAc                                                                                <  3:61379/149‑1 (MQ=255)                                                                         tGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAg                                                                          >  6:195911/1‑148 (MQ=255)                                                                           gCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAg                                                                          <  5:195911/146‑1 (MQ=255)                                                                           gCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCt                                                                        <  4:22479/148‑1 (MQ=255)                                                                             tGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACa                                                                                                                                     >  2:283086/1‑87 (MQ=255)                                                                             tGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACa                                                                                                                                     <  1:283086/87‑1 (MQ=255)                                                                               tttCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGc                                                                   >  2:259509/1‑149 (MQ=255)                                                                                    tGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATc                                                                                                                             >  4:247680/1‑88 (MQ=255)                                                                                      aCGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTCCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTCGGCAACCG‑‑ACGCACTACCAGCTGGTTCGACCGCt                                                            <  8:141053/149‑1 (MQ=255)                                                                                          gCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCaa                                                        <  2:44993/149‑1 (MQ=255)                                                                                                tCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAAcccc                                                                                                               >  2:273345/1‑90 (MQ=255)                                                                                                       cTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTggg                                                                                              <  1:262924/100‑1 (MQ=255)                                                                                                       cTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTggg                                                                                              >  2:262924/1‑100 (MQ=255)                                                                                                                       ttGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                           <  1:269677/149‑1 (MQ=255)                                                                                                                       ttGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                           >  2:269677/1‑149 (MQ=255)                                                                                                                         gttccGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTg                                                                    >  4:64802/4‑106 (MQ=21)                                                                                                                         gttccGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTg                                                                    <  3:64802/103‑1 (MQ=21)                                                                                                                          tACCCGCCTCATCGGCGGTAAACCCGCGTAATGGCCGGTAAATGTCGTATACATGTACAAACCCACCTTCAACCGCTTTGGCAAACGTCACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                           <  1:259509/148‑1 (MQ=255)                                                                                                                             cGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAgg                     <  2:252781/149‑1 (MQ=255)                                                                                                                               gccctttCGCGGCTTTACCCGCGAAGAGGCCGGTACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggc                   >  6:16345/6‑149 (MQ=38)                                                                                                                                          ggTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCt        <  6:112072/149‑1 (MQ=255)                                                                                                                                               aGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGAcc    >  2:163688/1‑148 (MQ=255)                                                                                                                                                gCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACCCa  >  7:167605/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                       AAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCTGACCCA  >  NC_002947/3265973‑3266264

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