Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,594,817	Δ1 bp	7.8%	coding (466/606 nt)	PP_2272 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,594,817	0	C	.	7.8%	157.7 / 6.8	51	coding (466/606 nt)	PP_2272	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base . (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

CAGAACGTACAGGGCTTCCCCGTCCAGTTGATCCTGGCCGACTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCTGGAGCAAATCCGTG  >  NC_002947/2594670‑2594949                                                                                                                                                     |                                                                                                                                     cAGAACGTACAGGGCTTCCCCGTCCAGTTGATCCTGGCCGACTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCg                                                                                                                                     >  6:246784/1‑149 (MQ=255)      cGTACAGGGCTTCCCCGTCCAGTTGATCCTGGCCGACTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACa                                                                                                                                <  5:246784/149‑1 (MQ=255)                  ccccGTCCAGTTGATCCTGGCCGACTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGaa                                                                                                                    >  8:15235/1‑149 (MQ=255)                    ccGTCCAGTTGATCCTGGCCGACTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAAtc                                                                                                                  >  1:98698/1‑149 (MQ=255)                                     gCCGACTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGc                                                                                                   <  8:261088/147‑1 (MQ=255)                                        gACTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCt                                                                                              <  6:323875/149‑1 (MQ=255)                                         aCTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTg                                                                                             >  8:127302/1‑149 (MQ=255)                                         aCTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTg                                                                                             <  7:155171/149‑1 (MQ=255)                                          cTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTg                                                                                             >  3:280523/1‑148 (MQ=255)                                              aTGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGgaaga                                                                                        <  4:280523/149‑1 (MQ=255)                                              aTGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGgaag                                                                                         >  7:3993/1‑148 (MQ=255)                                                 aCCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAAcg                                                                                     <  2:98698/149‑1 (MQ=255)                                                   cccGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACCAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAAcgcc                                                                                   >  7:98597/1‑149 (MQ=255)                                                                  cGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGcgcg                                                                    >  4:256653/1‑149 (MQ=255)                                                                     ttcgttcgACCTGAACCTGTCCCGCGCCACGTAGTGACGCATACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCATCCGCTGCGAAAGCGCCCTGGAAGAACGCCGCTCATTCAAGCGCGt                                                                   >  1:205605/1‑148 (MQ=255)                                                                         ttcgACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCg                                                             <  3:256653/149‑1 (MQ=255)                                                                           cgACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGc                                                            <  7:127302/148‑1 (MQ=255)                                                                             aCCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCCACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTg                                                         >  5:277331/1‑149 (MQ=255)                                                                               cTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAATAGGGCGTG‑GGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGAc                                                                                                >  3:70413/1‑107 (MQ=255)                                                                               cTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAATAGGGCGTG‑GGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGAc                                                                                                >  3:74753/1‑107 (MQ=255)                                                                               cTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAATAGGGCGTG‑GGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGAc                                                                                                <  4:74753/107‑1 (MQ=255)                                                                               cTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAATAGGGCGTG‑GGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGAc                                                                                                <  4:70413/107‑1 (MQ=255)                                                                               cTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGc                                                        >  3:17130/1‑148 (MQ=255)                                                                                tGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGc                                                        >  4:88625/1‑147 (MQ=255)                                                                                 gAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGaa                                                                                                                    <  3:282893/86‑1 (MQ=255)                                                                                 gAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGaa                                                                                                                    >  4:282893/1‑86 (MQ=255)                                                                                     cTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGaa                                                 <  3:304059/149‑1 (MQ=255)                                                                                     cTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGaa                                                 <  7:15235/149‑1 (MQ=255)                                                                                     cTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGaa                                                 <  8:3993/149‑1 (MQ=255)                                                                                        tCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTa                                              >  3:106309/1‑149 (MQ=255)                                                                                           cgcgCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTAtctc                                                              <  8:258553/130‑1 (MQ=255)                                                                                           cgcgCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTAtctc                                                              >  7:258553/1‑130 (MQ=255)                                                                                            gcgcCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTAcccc                                          >  2:261658/1‑149 (MQ=255)                                                                                               ccACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAg                                       <  2:205605/149‑1 (MQ=255)                                                                                                cACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGt                                      <  2:62052/149‑1 (MQ=255)                                                                                                   gTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTAc                                             <  1:13325/139‑1 (MQ=255)                                                                                                   gTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTAc                                             >  2:13325/1‑139 (MQ=255)                                                                                                         aCGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGGGAAAGCGACCTGGAATAACGCCGCTAATTCAAGAGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCa                               <  8:98597/148‑1 (MQ=255)                                                                                                           ggCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCacga                            >  4:53727/1‑149 (MQ=255)                                                                                                               aCTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCacga                            >  2:26729/1‑145 (MQ=255)                                                                                                               aCTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCAGCAAGAGAAGTACCCCGAGTACATCCacga                            <  1:26729/145‑1 (MQ=255)                                                                                                                   ccTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAacca                    >  1:215136/1‑149 (MQ=255)                                                                                                                       cACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCt                >  4:218671/1‑149 (MQ=255)                                                                                                                        aCGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCTg               <  4:17130/149‑1 (MQ=255)                                                                                                                           tttCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCTGGAg            <  3:218671/149‑1 (MQ=255)                                                                                                                                gcgGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACa                                   >  6:276277/1‑121 (MQ=255)                                                                                                                                gcgGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACa                                   <  5:276277/121‑1 (MQ=255)                                                                                                                                  ggCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCTGGAGCAAATcc     <  6:277331/149‑1 (MQ=255)                                                                                                                                    ccTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCTGGAGCAAATCCgt   <  2:215136/149‑1 (MQ=255)                                                                                                                                    ccTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCTGGAGCAAATCCgt   >  5:58162/1‑149 (MQ=255)                                                                                                                                     cTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCTGGAGCAAATCCgtg  >  2:62220/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                     CAGAACGTACAGGGCTTCCCCGTCCAGTTGATCCTGGCCGACTTCATGACCCCGCAAGAGCTGGCCGATTCGTTCGACCTGAACCTGTCCCGCGCCACGTA‑TGACGGCTACTCCCTGCACGTTTCCGCGGCCTTCGAAGAGGGCGTGCGGGACAAGGCATTCCGAATCCTCCGCTGCGAAAGCGACCTGGAAGAACGCCGCTCATTCAAGCGCGTATCTCGCCTGCAAGAGAAGTACCCCGAGTACATCCACGACGAAACCACCCTGGAGCAAATCCGTG  >  NC_002947/2594670‑2594949

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