Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,244,009	G→T	8.0%	R14L (CGC→CTC)	PP_1980 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,244,009	0	G	T	8.0%	157.1 / 4.1	50	R14L (CGC→CTC)	PP_1980	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base T (2/2);  total (25/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.06e-01

TGATTCGCGGGCAAGCCCGCCCCTATAACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGATGGCTGTGGTG  >  NC_002947/2243862‑2244157                                                                                                                                                    |                                                                                                                                                     tGATTCGCGGGCAAGCCCGCCCCTATAACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTgcgc                                                                                                                                                     >  6:278948/1‑149 (MQ=255)  gATTCGCGGGCAAGCCCGCCCCTATAACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGcgcg                                                                                                                                                    >  8:152576/1‑149 (MQ=255)   aTTCGCGGGCAAGCCCGCCCCTATAACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTgcgc                                                                                                                                                     <  8:15260/147‑1 (MQ=255)   aTTCGCGGGCAAGCCCGCCCCTATAACGACTGCATCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGt                                                                                                                                                   <  5:306060/149‑1 (MQ=255)    ggggcgGGCTTGCCCGCCCCTATAACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTc                                                                                                                                                  >  4:284204/4‑149 (MQ=255)            cAAGCCCGCCCCTATAACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGa                                                                                                                                          <  1:194249/149‑1 (MQ=255)             aaGCCCGCCCCTATAACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAg                                                                                                                                         <  7:95761/149‑1 (MQ=255)                         ataACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGc                                                                                                                             >  6:111347/1‑149 (MQ=255)                         ataACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGc                                                                                                                             >  3:278153/1‑149 (MQ=255)                          taACGACTGCCTCAGGCTCACGGGAGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACgg                                                                                                                              <  4:5370/147‑1 (MQ=255)                              gACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTgg                                                                                                                        <  6:81020/149‑1 (MQ=255)                                 tGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGcgcg                                                                                                                                                    <  8:143592/118‑1 (MQ=255)                                 tGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGcgcg                                                                                                                                                    >  7:143592/1‑118 (MQ=255)                                   ccTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATg                                                                                                                   >  7:44820/1‑149 (MQ=255)                                    cTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGc                                                                                                                  <  4:278153/149‑1 (MQ=255)                                    cTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGc                                                                                                                  >  2:66148/1‑149 (MQ=255)                                              cGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGcc                                                                                                        <  7:152576/149‑1 (MQ=255)                                               gggggCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGAc                                                                                                                                <  6:292715/124‑1 (MQ=255)                                               gggggCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGAc                                                                                                                                >  5:292715/1‑124 (MQ=255)                                                 gggCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCcaacaa                                                                                                            >  5:13571/1‑142 (MQ=255)                                                 gggCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCcaacaa                                                                                                            <  6:13571/142‑1 (MQ=255)                                                       gACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGt                                                                                                >  5:312260/1‑148 (MQ=255)                                                       gACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTc                                                                                               <  5:2924/149‑1 (MQ=255)                                                             gCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTgg                                                                                         >  6:255692/1‑149 (MQ=255)                                                              ccccGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGc                                                                                        <  7:101551/149‑1 (MQ=255)                                                              ccccGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGc                                                                                        <  5:278948/149‑1 (MQ=255)                                                               cccGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCt                                                                                       >  7:239386/1‑149 (MQ=255)                                                                ccGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCtt                                                                                      <  3:284204/149‑1 (MQ=255)                                                                       tCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCt                                                                               >  6:28936/1‑149 (MQ=255)                                                                         gctgctAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTgc                                                                             <  6:312260/149‑1 (MQ=255)                                                                          ctgctAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTgcg                                                                            >  7:161686/1‑149 (MQ=255)                                                                                aaTCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTgg                                                                       <  5:111347/148‑1 (MQ=255)                                                                                       ttttGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTc                                                               >  6:220249/1‑149 (MQ=255)                                                                                                 aGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTg                                                     <  6:313829/149‑1 (MQ=255)                                                                                                     aGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGc                                                 <  8:44820/149‑1 (MQ=255)                                                                                                       cGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTg                                               <  5:28936/149‑1 (MQ=255)                                                                                                          aaaTGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGAc                                            <  8:239386/149‑1 (MQ=255)                                                                                                                 cTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGcc                                      >  3:281847/1‑148 (MQ=255)                                                                                                                    ggACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAg                                   >  3:314875/1‑148 (MQ=255)                                                                                                                          ggCAACACCCTTCGTCATCGACCTGCTCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGgc                            >  2:314598/1‑149 (MQ=255)                                                                                                                          ggCAACACCCTTCGTCATCGACCTGCTCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGgc                            >  2:319431/1‑149 (MQ=255)                                                                                                                           gCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGgc                            <  5:220249/148‑1 (MQ=255)                                                                                                                             aaCACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGcc                         >  8:196629/1‑149 (MQ=255)                                                                                                                               cacCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCtt                       >  3:169958/1‑149 (MQ=255)                                                                                                                                acCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCtt                       >  2:285483/1‑148 (MQ=255)                                                                                                                                  ccTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGAc                    >  5:239705/1‑149 (MQ=255)                                                                                                                                  ccTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGAc                    >  2:72496/1‑149 (MQ=255)                                                                                                                                      cGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCgg                 <  1:285483/148‑1 (MQ=255)                                                                                                                                          aTCGACCTGCTCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGAt            <  1:319431/149‑1 (MQ=255)                                                                                                                                          aTCGACCTGCTCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGAt            <  1:314598/149‑1 (MQ=255)                                                                                                                                          aTCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGAt            >  2:69176/1‑149 (MQ=255)                                                                                                                                          aTCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGAt            <  1:66148/149‑1 (MQ=255)                                                                                                                                          aTCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGAt            >  1:59559/1‑149 (MQ=255)                                                                                                                                             gACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGATgg          <  8:161686/148‑1 (MQ=255)                                                                                                                                                cTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGATGGCtgt      >  5:192432/1‑149 (MQ=255)                                                                                                                                                    gcgTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGATGGCTgtggtg  >  4:218844/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGATTCGCGGGCAAGCCCGCCCCTATAACGACTGCCTCAGGCTCACGGGGGCTGGACATTGCCCCGCAGATCGCTGCTAAATCGCCTTTTGCGCCAAGGAAGCGAAAATGAGCTGGGACCTGGCAACACCCTTCGTCATCGACCTGCGCGTCGGCAGCGAGGACATCGACGGCCTTGGCCATGCCAACAATGCCGTCTACGTCACCTGGCTTGAGCGCTGCGCCTGGCGCCACTCCCAGCGTCTGGGGCTGGACCTGGCCGAGTACCGGCGCCTTGACCGGGCGATGGCTGTGGTG  >  NC_002947/2243862‑2244157

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