Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,845,975	T→A	7.8%	A47A (GCA→GCT)	PP_2497 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,845,975	0	T	A	7.8%	149.4 / 3.3	51	A47A (GCA→GCT)	PP_2497	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base A (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 = 6.68e-01

TTGTGGGCCAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAATGACATGGGG  >  NC_002947/2845829‑2846119                                                                                                                                                   |                                                                                                                                                  ttGTGGGCCAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGtgctgc                                                                                                                                                 <  4:324151/149‑1 (MQ=255)        ccAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCG‑GTGCAGCACTCGGcc                                                                                                                                         <  7:51446/149‑1 (MQ=255)        ccAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCG‑GTGCAGCACTCGGcc                                                                                                                                         <  7:46200/149‑1 (MQ=255)          aGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAc                                                                                                                                               <  8:144674/142‑1 (MQ=255)          aGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAc                                                                                                                                               >  7:144674/1‑142 (MQ=255)            gACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGTCCGAGTGCTGCAC‑CGGCCTTgc                                                                                                                                     <  3:126666/149‑1 (MQ=255)            gACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTgc                                                                                                                                     <  4:192634/149‑1 (MQ=255)                      tCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGc                                                                                                                              >  3:30965/1‑146 (MQ=255)                        cAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCagag                                                                                                                         >  2:95436/1‑149 (MQ=255)                                              ccAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAg                                                                                                   <  3:345390/149‑1 (MQ=255)                                                  tCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGAGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCt                                                                                               <  7:132716/149‑1 (MQ=255)                                                     cacaACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGagca                                                                                                                       >  6:361857/1‑122 (MQ=255)                                                     cacaACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGagca                                                                                                                       <  5:361857/122‑1 (MQ=255)                                                        aaCCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCagag                                                                                                                         >  2:324055/1‑117 (MQ=255)                                                        aaCCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCagag                                                                                                                         <  1:324055/117‑1 (MQ=255)                                                        aaCCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACa                                                                                         <  5:218705/149‑1 (MQ=255)                                                            accGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAgcgc                                                                                     >  2:132753/1‑149 (MQ=255)                                                            accGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAgcgc                                                                                     >  5:106700/1‑149 (MQ=255)                                                              cGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGcc                                                                                    >  8:123894/1‑148 (MQ=255)                                                                       cGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGACGAGTGCTGCCC‑CGGCCTTGCCGCCTGCCAGCGCAGCGATGACGCCCCGTGCGAGGCCTCCGCCAGCGCCGATGGCcgcg                                                                          >  7:340924/1‑149 (MQ=255)                                                                        ggCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCg                                                                                                     >  3:315937/1‑121 (MQ=255)                                                                        ggCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCg                                                                                                     <  4:315937/121‑1 (MQ=255)                                                                         gCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGcc                                                                        >  1:70586/1‑149 (MQ=255)                                                                               ccgccgAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCttt                                                                  <  1:95436/149‑1 (MQ=255)                                                                                   cgAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATg                                                                                <  1:259518/131‑1 (MQ=255)                                                                                   cgAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATg                                                                                >  2:259518/1‑131 (MQ=255)                                                                                       aaTGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCg                                                           >  2:350217/1‑148 (MQ=255)                                                                                       aaTGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCg                                                           <  1:350217/148‑1 (MQ=255)                                                                                           gcgcCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATg                                                      <  2:176516/149‑1 (MQ=255)                                                                                                ggATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATg                                                      <  8:143349/144‑1 (MQ=255)                                                                                                ggATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATg                                                      >  7:143349/1‑144 (MQ=255)                                                                                                  aTACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGa                                                >  8:119713/1‑148 (MQ=255)                                                                                                   tACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGaa                                               <  7:119713/148‑1 (MQ=255)                                                                                                      cgagccgagAGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGa                                           >  3:200786/1‑149 (MQ=255)                                                                                                       gagccgagGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGAc                                          <  6:106700/149‑1 (MQ=255)                                                                                                            gagGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGCCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTg                                     >  7:210149/1‑149 (MQ=255)                                                                                                               ggTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTggg                                   <  7:123894/148‑1 (MQ=255)                                                                                                                gTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTggg                                   >  6:91304/1‑147 (MQ=255)                                                                                                                        gACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGc                         >  5:105101/1‑149 (MQ=255)                                                                                                                         aCAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCg                        <  6:105101/149‑1 (MQ=255)                                                                                                                               ccGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGa                             >  3:1656/1‑138 (MQ=255)                                                                                                                               ccGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATCAAGCGACTTTTGGGGAGTGa                             <  4:1656/138‑1 (MQ=255)                                                                                                                                cGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCAc                 >  8:215735/1‑149 (MQ=255)                                                                                                                                 ggCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACt                <  5:91304/149‑1 (MQ=255)                                                                                                                                   cggcaaGGCCG‑GTGCAGCACTCGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTg               >  5:269549/6‑148 (MQ=255)                                                                                                                                    gATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTggg                                   <  7:46956/127‑1 (MQ=255)                                                                                                                                    gATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTggg                                   >  8:46956/1‑127 (MQ=255)                                                                                                                                      caaGGCCG‑GTGCAGCACTCGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAAt           >  5:84057/3‑149 (MQ=255)                                                                                                                                        aGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAATg          <  1:132753/148‑1 (MQ=255)                                                                                                                                         ggCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAATGa         >  7:37949/1‑148 (MQ=255)                                                                                                                                          gCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAATGACa       >  3:269635/1‑149 (MQ=255)                                                                                                                                               gTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAATGACATgggg  <  2:70586/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TTGTGGGCCAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAATGACATGGGG  >  NC_002947/2845829‑2846119

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