Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,845,970	Δ1 bp	7.5%	coding (146/330 nt)	PP_2497 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,845,970	0	A	.	7.5%	163.9 / 7.4	53	coding (146/330 nt)	PP_2497	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base . (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

GCAGAGGTTGTGGGCCAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAATGACA  >  NC_002947/2845822‑2846114                                                                                                                                                     |                                                                                                                                                  gCAGAGGTTGTGGGCCAGGACATGACGGTCCAGCGACTGGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGa                                                                                                                                                   <  7:276304/149‑1 (MQ=255)        ttGTGGGCCAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGtgctgc                                                                                                                                            <  4:324151/149‑1 (MQ=255)               ccAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCG‑GTGCAGCACTCGGcc                                                                                                                                    <  7:51446/149‑1 (MQ=255)               ccAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCG‑GTGCAGCACTCGGcc                                                                                                                                    <  7:46200/149‑1 (MQ=255)                cAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGtgc                                                                                                                                               >  1:214343/1‑138 (MQ=255)                cAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGtgc                                                                                                                                               <  2:214343/138‑1 (MQ=255)                 aGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAc                                                                                                                                          >  7:144674/1‑142 (MQ=255)                 aGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAc                                                                                                                                          <  8:144674/142‑1 (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                                                                                                                    <  1:324055/117‑1 (MQ=255)                                                               aaCCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCagag                                                                                                                    >  2:324055/1‑117 (MQ=255)                                                               aaCCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACa                                                                                    <  5:218705/149‑1 (MQ=255)                                                                   accGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAgcgc                                                                                >  5:106700/1‑149 (MQ=255)                                                                   accGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAgcgc                                                                                >  2:132753/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                                                                           >  2:259518/1‑131 (MQ=255)                                                                                          cgAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATg                                                                           <  1:259518/131‑1 (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                              >  8:46956/1‑127 (MQ=255)                                                                                                                                           gATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTggg                              <  7:46956/127‑1 (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)                                                                                                                                                     |                                                                                                                                                  GCAGAGGTTGTGGGCCAGGACATGACGGTCCAGCGACTCGCCGAGCTGCGCGCCAATCGCACAACCACCGACGCTTCCGGCCAGGCCGCCGAGAATGGCGCCGGATACCGAGCCGAGGGTGATGCCGACAGGGCCGGCGATTAGGCCGAGTGCTGCAC‑CGGCCTTGCCGCCTGCCAGAGCAGCGCTGACGCCACGTGCGAGGCCTCCGACAGCGCCGATGGCCGCGCCGGCTTTCATGGCGTGATGGAATGAAGCGACTTTTGGGGAGTGACAGCGGGGGCACTGCAATGACA  >  NC_002947/2845822‑2846114

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