Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,240,982	C→T	9.8%	G239G (GGG→GGA)	PP_3718 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,240,982	0	C	T	9.8%	139.6 / 4.8	51	G239G (GGG→GGA)	PP_3718	aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base T (3/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

TGGGGTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTG  >  NC_002947/4240835‑4241129                                                                                                                                                    |                                                                                                                                                    tGGGGTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCAccc                                                                                                                                                    >  8:59424/1‑149 (MQ=255)   gggTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACtccg                                                                                                                                                  >  3:185347/1‑149 (MQ=255)        aCGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACcccgcccg                                                                                                                                              <  3:240690/148‑1 (MQ=255)          gCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGGGTGATCATTCCGCCCGAgg                                                                                                                                           <  8:57948/149‑2 (MQ=255)          gCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGc                                                                                                                                           >  1:141191/1‑149 (MQ=255)            cAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGTCGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGcc                                                                                                                                          >  6:40427/1‑148 (MQ=255)            cAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGcc                                                                                                                                          >  1:126299/1‑148 (MQ=255)             aaaCAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCa                                                                                                                                        >  6:102803/1‑149 (MQ=255)                 aGTTCCTCGCTGGCGAAGAAAGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACACCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCGGCCCGAGCCCATGAt                                                                                                                                    <  5:278543/149‑1 (MQ=255)                   ttCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGAtcg                                                                                                                                  <  4:139319/149‑1 (MQ=255)                     ccTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGAtcggc                                                                                                                                >  8:229099/1‑149 (MQ=255)                         gCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGAtcggctcg                                                                                                                             <  7:259165/148‑1 (MQ=255)                             gCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCgaa                                                                                                                        <  3:64738/149‑1 (MQ=255)                               gaagaaGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCgaag                                                                                                                       >  8:212922/1‑148 (MQ=255)                                   aagGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgc                                                                                                                  <  2:203498/149‑1 (MQ=255)                                     ggTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACTCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgccg                                                                                                                <  4:185347/149‑1 (MQ=255)                                        gCCCAGGCGCCCGAACGAGGTCACCACTTCGGCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCa                                                                                                             >  1:226965/1‑149 (MQ=255)                                         cccAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCAc                                                                                                            <  5:102803/149‑1 (MQ=255)                                               cgcCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTc                                                                                                      <  7:229099/149‑1 (MQ=255)                                                 cccGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGGAGGTCTGGCACAGCTGCCACATGCGCGGGAAGAAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCgg                                                                                                    <  6:284323/149‑1 (MQ=255)                                                              accacTTGGTCGGCGACGAACAGGATGGCGTACGTCTGGCACAGCTGCCACATGCGGTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTGGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGagag                                                                                       <  8:285485/149‑1 (MQ=255)                                                              accacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgccgcc                                                                                                              >  7:2063/1‑126 (MQ=255)                                                              accacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATgccgcc                                                                                                              <  8:2063/126‑1 (MQ=255)                                                                cacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGa                                                                                     >  1:184656/1‑149 (MQ=255)                                                                 tgttCGTCGCCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGa                                                                                     >  4:100960/3‑148 (MQ=255)                                                                             acgaACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGaa                                                                        <  2:141191/149‑1 (MQ=255)                                                                               gaACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGAGATTAATTCCGCCACGTTGTCGGCACCCAGAGAGAGGATATTATCCTCGAAca                                                                      >  6:5248/1‑148 (MQ=255)                                                                                 aCAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCg                                                                    <  2:126299/149‑1 (MQ=255)                                                                                     gATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCgg                                                                   >  1:11045/1‑146 (MQ=255)                                                                                     gATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCgg                                                                   <  2:11045/146‑1 (MQ=255)                                                                                     gATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGcca                                                                >  3:138306/1‑149 (MQ=255)                                                                                           gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAgg                                                           >  3:290996/1‑148 (MQ=255)                                                                                           gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAgg                                                           >  5:83263/1‑148 (MQ=255)                                                                                              ggTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAg                                                        <  2:226965/148‑1 (MQ=255)                                                                                               gTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTc                                                      >  7:304597/1‑149 (MQ=255)                                                                                               gTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTc                                                      <  3:286333/149‑1 (MQ=255)                                                                                                    gCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGa                                                    <  4:138306/146‑1 (MQ=255)                                                                                                      acaGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTg                                                                                <  3:333510/116‑1 (MQ=255)                                                                                                      acaGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTg                                                                                >  4:333510/1‑116 (MQ=255)                                                                                                         gCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTcc                                            <  7:59424/149‑1 (MQ=255)                                                                                                             ccACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGc                                          <  3:100960/147‑1 (MQ=255)                                                                                                              cacaTGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTc                                       <  5:40427/149‑1 (MQ=255)                                                                                                                       tgcaaGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGATGTCGGCACCCCGCGAGAGGATCTTGTCCTCGAACTCGGCCCCCAGGAAGTCGAGGAACTCCGCCTCGTCCATgtc                              >  7:266696/4‑149 (MQ=255)                                                                                                                       tGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATgtc                              >  4:364915/1‑149 (MQ=255)                                                                                                                         gAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATgtcgt                            >  5:32947/1‑149 (MQ=255)                                                                                                                          aaGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATgtcgtc                           <  3:364915/149‑1 (MQ=255)                                                                                                                                  cTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGgcgcgc                   >  2:111125/1‑149 (MQ=255)                                                                                                                                   tCGGGCGGGGTGATCATTCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGcgcgcg                  >  2:35335/1‑149 (MQ=255)                                                                                                                                     gggCGGGGTGATCATTCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTGCATGTCGTGTGgcgggcgg                 >  8:145107/1‑148 (MQ=255)                                                                                                                                     gggCGGGGTGATCATTCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGAgg                 >  6:336191/1‑148 (MQ=255)                                                                                                                                     gggCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGt                <  7:212922/149‑1 (MQ=255)                                                                                                                                       gCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTAc              >  2:336201/1‑149 (MQ=255)                                                                                                                                          gAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACggg           <  4:290996/149‑1 (MQ=255)                                                                                                                                               aTCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttggg      <  8:304597/149‑1 (MQ=255)                                                                                                                                                   ccccGCCCGAGCGCATGATCGGCTCGGCGAAGAATGCGGGCAGGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGATGTCGAGGAACTCCGCCTCGTGCATGTCGTCTGGCGCGCGGTACgggttggggttg  <  8:266696/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGGGGTGACGCCAAACAGTTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCATGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTG  >  NC_002947/4240835‑4241129

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