Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,266,113	A→T	15.1%	intergenic (+30/+29)	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,113	0	A	T	15.1%	87.1 / 9.6	40	intergenic (+30/+29)	PP_2863/PP_2864	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base T (2/4);  total (17/23)
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.14e-01

AAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCT  >  NC_002947/3265973‑3266258                                                                                                                                             |                                                                                                                                                    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                                                                                        >  2:262924/1‑100 (MQ=255)                                                                                                       cTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTggg                                                                                        <  1:262924/100‑1 (MQ=255)                                                                                                                       ttGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                     >  2:269677/1‑149 (MQ=255)                                                                                                                       ttGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                     <  1:269677/149‑1 (MQ=255)                                                                                                                         gttccGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTg                                                              <  3:64802/103‑1 (MQ=21)                                                                                                                         gttccGGCCCTTTCGCGGCTTTACCCGCGAAGAGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTg                                                              >  4:64802/4‑106 (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)                                                                                                                                             |                                                                                                                                                    AAGACCCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCG‑‑ACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATGCT  >  NC_002947/3265973‑3266258

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