Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,240,997	G→C	7.4%	I234M (ATC→ATG)	PP_3718 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,240,997	0	G	C	7.4%	165.6 / 3.5	54	I234M (ATC→ATG)	PP_3718	aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base C (2/2);  total (28/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 = 8.69e-01

TTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACGTGGTGGAT  >  NC_002947/4240853‑4241143                                                                                                                                                  |                                                                                                                                                   ttCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCGA‑TCAtgg                                                                                                                                                >  7:262992/1‑147 (MQ=255)     tCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGAtcggctc                                                                                                                                            >  1:280176/1‑149 (MQ=255)     tCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGAtcggctc                                                                                                                                            >  4:195104/1‑149 (MQ=255)        cTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCGA‑TCAtgggctcggg                                                                                                                                         <  8:262992/149‑2 (MQ=255)              aagaagGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCgaagaa                                                                                                                                   >  8:162208/1‑149 (MQ=255)              aagaagGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCgaagaa                                                                                                                                   >  7:20827/1‑149 (MQ=255)                   ggTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATgccg                                                                                                                              >  3:78262/1‑149 (MQ=255)                       cccAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCAc                                                                                                                          >  8:250547/1‑149 (MQ=255)                           ggCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTg                                                                                                                      >  7:315498/1‑149 (MQ=255)                           ggCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTg                                                                                                                      <  7:162208/149‑1 (MQ=255)                              gcCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCg                                                                                                                   >  2:83220/1‑149 (MQ=255)                              gcCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCg                                                                                                                   >  7:171075/1‑149 (MQ=255)                                 cGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGc                                                                                                                 >  2:311318/1‑148 (MQ=255)                                        ggTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCg                                                                                                         <  7:250547/149‑1 (MQ=255)                                            accacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCACTGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCgaga                                                                                                      <  1:83220/149‑1 (MQ=255)                                              cacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCgaga                                                                                                      <  3:256829/146‑1 (MQ=255)                                              cacTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCgaga                                                                                                      >  4:256829/1‑146 (MQ=255)                                                 ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCgaa                                                                                                                                      <  6:29523/111‑1 (MQ=255)                                                 ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCgaa                                                                                                                                      >  5:29523/1‑111 (MQ=255)                                                 ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATgccgcc                                                                                                                            >  5:310347/1‑121 (MQ=255)                                                 ttCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATgccgcc                                                                                                                            <  6:310347/121‑1 (MQ=255)                                                   cgtcgGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCtt                                                                                               <  2:280176/148‑1 (MQ=255)                                                     tcgGCGACGAACAGGATGTCTTAGGTCTGGCACAGCTGCCAGATGCGCTGGAAGTAGCCCTCGGGCGTAATGATCACCCCGCCCGAGCCCA‑TTATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTc                                                                                            <  2:306616/149‑1 (MQ=255)                                                       ggCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTcc                                                                                           <  8:171075/148‑1 (MQ=255)                                                        gCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTc                                                                                         <  3:136433/149‑1 (MQ=255)                                                          gacgaACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCg                                                                                        >  3:282280/1‑148 (MQ=255)                                                               aCAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCg                                                                                  >  5:138675/1‑149 (MQ=255)                                                                 aGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGc                                                                                <  1:311318/149‑1 (MQ=255)                                                                         gTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTc                                                                                         >  4:129370/1‑132 (MQ=255)                                                                             gTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATc                                                                                                 >  8:156501/1‑120 (MQ=255)                                                                             gTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATc                                                                                                 <  7:156501/120‑1 (MQ=255)                                                                              tCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAATAATGCCGCCACGTTGTCGGCACCCAGCgagatg                                                                                                    <  7:197567/116‑3 (MQ=255)                                                                              tCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAgg                                                                                                    >  8:197567/1‑116 (MQ=255)                                                                                  gCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGa                                                                  <  6:138675/146‑1 (MQ=255)                                                                                  gCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAgg                                                                <  8:315498/148‑1 (MQ=255)                                                                                  gCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAg                                                                 <  4:282280/147‑1 (MQ=255)                                                                                                    cTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATgt                                             <  6:74485/149‑1 (MQ=255)                                                                                                          gTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTg                                       <  5:132598/149‑1 (MQ=255)                                                                                                                   gggCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCgg                               >  4:251162/1‑148 (MQ=255)                                                                                                                       ggAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAAc                                                                                     <  5:13431/90‑1 (MQ=255)                                                                                                                       ggAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAAc                                                                                     >  6:13431/1‑90 (MQ=255)                                                                                                                        gAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGCCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGcc                                                       <  5:153858/119‑1 (MQ=255)                                                                                                                        gAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGCCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGcc                                                       >  6:153858/1‑119 (MQ=255)                                                                                                                         aaTGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggt                        >  6:221605/1‑149 (MQ=255)                                                                                                                           tGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttg                      >  6:270314/1‑149 (MQ=255)                                                                                                                             aTCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttggg                    >  4:130834/1‑149 (MQ=255)                                                                                                                                aCCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTAc                            <  3:158588/138‑1 (MQ=255)                                                                                                                                 ccccGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACgggttggggttg                >  5:38734/1‑149 (MQ=255)                                                                                                                                   ttcgccgAGCCGA‑TCATGGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTcc                                                          <  1:173829/101‑1 (MQ=255)                                                                                                                                   ttcgccgAGCCGA‑TCATGGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTcc                                                          >  2:173829/5‑105 (MQ=255)                                                                                                                                          aGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACgtgg       <  5:270314/149‑1 (MQ=255)                                                                                                                                           gCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACgtggt      <  5:221605/149‑1 (MQ=255)                                                                                                                                           gCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACgtggt      <  6:38734/149‑1 (MQ=255)                                                                                                                                                 tGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACGTGGTGGAt  >  6:354964/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TTCCTCGCTGGCGAAGAAGGTGCCCAGGCGCCCGAACGAGGTCACCACTTCGTCGGCGACGAACAGGATGTCGTAGGTCTGGCACAGCTGCCACATGCGCTGGAAGTAGCCCTCGGGCGGAATGATCACCCCGCCCGAGCCCA‑TGATCGGCTCGGCGAAGAATGCCGCCACGTTGTCGGCACCCAGCGAGAGGATCTTGTCCTCGAACTCGGCCACCAGGAAGTCGAGGAACTCCGCCTCGTCCATGTCGTCTGGCGCGCGGTACGGGTTGGGGTTGGAGACGTGGTGGAT  >  NC_002947/4240853‑4241143

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