Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,359,451	C→G	7.3%	G169R (GGC→CGC)	PP_3834 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,359,451	0	C	G	7.3%	169.7 / 2.9	55	G169R (GGC→CGC)	PP_3834	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base G (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CTCATCGATAAGGCATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCACACCGTCGACGATCAGCAGGCTG  >  NC_002947/4359308‑4359599                                                                                                                                                |                                                                                                                                                     cTCATCGATAAGGCATTCGTGACGCCCCGACAACGCACCCCTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCtcgt                                                                                                                                                 <  7:342184/149‑1 (MQ=255)    aTCGATAAGGCATTCGTGATGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGCCCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCtcgtc                                                                                                                                                >  6:277774/1‑147 (MQ=255)    aTCGATAAGGCATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGGTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCGCtcgtc                                                                                                                                                >  1:57347/1‑147 (MQ=255)    aTCGATAAGGCATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCtcgtc                                                                                                                                                <  4:149696/147‑1 (MQ=255)    aTCGATAAGGCATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGATGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCGCtcgtc                                                                                                                                                >  1:62812/1‑147 (MQ=255)       gATAAGGCATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGTTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCGCTCGTCGAGc                                                                                                                                            <  2:62812/148‑1 (MQ=255)       gATAAGGCATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCGCTCGTCGAGc                                                                                                                                            <  2:57347/148‑1 (MQ=255)           aGGCATTCGTGACGCCCCGACAACGCACCACTTCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGc                                                                                                                                       <  5:169640/149‑1 (MQ=255)              cATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACt                                                                                                                                    <  8:304793/149‑1 (MQ=255)              cATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACt                                                                                                                                    >  3:265868/1‑149 (MQ=255)                         ccccGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACa                                                                                                                         <  1:238677/149‑1 (MQ=255)                           ccGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAg                                                                                                                        >  1:108327/1‑148 (MQ=255)                                aaCGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGt                                                                                                                  <  2:108327/149‑1 (MQ=255)                                     accacTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAAc                                                                                                             >  1:230617/1‑149 (MQ=255)                                       cacTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTACCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCg                                                                                                           >  4:2275/1‑149 (MQ=255)                                        acTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCgg                                                                                                          <  4:265868/149‑1 (MQ=255)                                         cTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGc                                                                                                         <  2:230617/149‑1 (MQ=255)                                            ccGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGaa                                                                                                      <  4:162133/149‑1 (MQ=255)                                                        cgcAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACTCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGgc                                                                                                   <  7:242417/140‑1 (MQ=255)                                                        cgcAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGgc                                                                                                   >  8:242417/1‑140 (MQ=255)                                                           aaTTGGACATGGCGCACGGCGCGGCTACCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTg                                                                                       <  3:2275/149‑1 (MQ=255)                                                             ttGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTg                                                                                                                                   <  2:314848/103‑1 (MQ=255)                                                                      gcgcACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAgaga                                                                            >  8:84788/1‑149 (MQ=255)                                                                            ggCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAAct                                                                      >  6:222402/1‑149 (MQ=255)                                                                                 ggCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCt                                                                                  <  1:143042/132‑1 (MQ=255)                                                                                 ggCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCt                                                                                  >  2:143042/1‑132 (MQ=255)                                                                                        gACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAAc                                                                       >  8:25617/1‑136 (MQ=255)                                                                                        gACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAAc                                                                       <  7:25617/136‑1 (MQ=255)                                                                                             gTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGcgcg                                                     <  7:84788/149‑1 (MQ=255)                                                                                              ttGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAActtc                                                                    >  5:145695/1‑133 (MQ=255)                                                                                              ttGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAActtc                                                                    <  6:145695/133‑1 (MQ=255)                                                                                                 cTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGttt                                                 >  4:24541/1‑149 (MQ=255)                                                                                                    aCCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTc                                                <  8:276344/147‑1 (MQ=255)                                                                                                    aCCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTc                                                >  7:276344/1‑147 (MQ=255)                                                                                                    aCCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCg                                              >  2:73574/1‑149 (MQ=255)                                                                                                      ccAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTg                                                                 >  8:117832/1‑128 (MQ=255)                                                                                                      ccAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTg                                                                 <  7:117832/128‑1 (MQ=255)                                                                                                      ccAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCa                                            <  1:73574/149‑1 (MQ=255)                                                                                                       cAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCgg                                                    >  5:319593/1‑140 (MQ=255)                                                                                                       cAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCgg                                                    <  6:319593/140‑1 (MQ=255)                                                                                                             gATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAgg                                                                                                    <  4:290967/86‑1 (MQ=255)                                                                                                             gATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAgg                                                                                                    >  3:290967/1‑86 (MQ=255)                                                                                                               tCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCTGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCgg                                   <  3:24541/149‑1 (MQ=255)                                                                                                                     cATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCAtt                              >  5:352003/1‑148 (MQ=255)                                                                                                                     cATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCAtt                              >  8:210885/1‑148 (MQ=255)                                                                                                                        cGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAg                           >  5:28564/1‑148 (MQ=255)                                                                                                                         gACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGc                          <  6:352003/148‑1 (MQ=255)                                                                                                                          aCCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCgg                                                    <  3:98332/121‑1 (MQ=255)                                                                                                                          aCCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGATCAGAGACGAACTTCTTGGTGTCCAGCGCgg                                                    >  4:98332/1‑121 (MQ=255)                                                                                                                            cATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGACCTTCTTGGTGTCCAGCGCGGGTTCCGCACGTCCCCGGGCATTGAGCacac                      >  4:58391/1‑149 (MQ=255)                                                                                                                            cATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCacac                      <  5:222402/149‑1 (MQ=255)                                                                                                                                     ttCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCACACCGTcgacga             >  1:30457/1‑149 (MQ=255)                                                                                                                                       ccGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCACACCGTc                  <  7:132987/142‑1 (MQ=255)                                                                                                                                       ccGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCACACCGTc                  >  8:132987/1‑142 (MQ=255)                                                                                                                                             tGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCACACCGTCGACGATCAGCAgg     >  8:321805/1‑149 (MQ=255)                                                                                                                                               cccTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCACACCGTCGACGATcagcag      >  1:149023/1‑146 (MQ=255)                                                                                                                                                ccTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCACACCGTCGACGATCAGCAGGCTg  >  7:230379/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     CTCATCGATAAGGCATTCGTGACGCCCCGACAACGCACCACTGCCGCGGGTCTTGCGCAATTGGACATGGCGCACGGCGCGGCTGCCGACCAGTTGCTCACCCAGTTCGATCACGCCATCGACCATGGTGTGTTCCGGGCTGCCCTCGTCGAGCCGGGCACTGGTCAGCAACAGTACGGTGCAACCGGCGAAGGCAGCATGCCCCTGCAGCTCAGAGACGAACTTCTTGGTGTCCAGCGCGGTTTCCGCACGTACCCGGGCATTGAGCACACCGTCGACGATCAGCAGGCTG  >  NC_002947/4359308‑4359599

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