Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,332,701	A→T	7.6%	L135Q (CTG→CAG)	PP_2930 ←	L‑serine dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,332,701	0	A	T	7.6%	150.3 / 3.2	53	L135Q (CTG→CAG)	PP_2930	L‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base T (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 = 9.26e-01

AGCACAGGCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCG  >  NC_002947/3332554‑3332841                                                                                                                                                     |                                                                                                                                              aGCACAGGCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCAcag                                                                                                                                              >  3:210315/1‑149 (MQ=255)     caGGCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGt                                                                                                                                          <  4:132228/149‑1 (MQ=255)     caGGCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGt                                                                                                                                          <  2:68631/149‑1 (MQ=255)      aGGCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACTGCtgtg                                                                                                                                         >  6:54385/1‑149 (MQ=255)      aGGCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCC‑CTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTgg                                                                                                                                        <  2:45713/149‑1 (MQ=255)       ggCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTgg                                                                                                                                        <  4:210315/149‑1 (MQ=255)       ggCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTgg                                                                                                                                        >  2:296189/1‑149 (MQ=255)        gCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTggg                                                                                                                                       <  1:19943/149‑1 (MQ=255)          acaTCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACTGCTGTgggcc                                                                                                                                     <  5:54385/149‑3 (MQ=255)                ccccAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGaa                                                                                                                              >  5:230665/1‑149 (MQ=255)                 cccAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGCCACCACCCCCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCCCAGCAGTGGGTCGT‑CGAAt                                                                                                                             >  4:17927/1‑149 (MQ=255)                   cAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATgg                                                                                                                           <  6:230665/149‑1 (MQ=255)                       cATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTg                                                                                                                        >  4:110173/1‑148 (MQ=255)                               ttGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAgg                                                                                                                >  2:32045/1‑148 (MQ=255)                                     cccTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGACCCACTGCTGTggg                                                                                                                                       >  7:6625/1‑120 (MQ=255)                                     cccTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGACCCACTGCTGTggg                                                                                                                                       <  8:6625/120‑1 (MQ=255)                                      ccTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAg                                                                                                                 >  8:107560/1‑140 (MQ=255)                                      ccTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAg                                                                                                                 <  7:107560/140‑1 (MQ=255)                                      ccTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTc                                                                                                        >  7:273053/1‑149 (MQ=255)                                          caccacGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGc                                                                                                    <  2:243058/149‑1 (MQ=255)                                          caccacGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGc                                                                                                    <  2:35453/149‑1 (MQ=255)                                                        gCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCtgt                                                                                      <  1:32045/149‑1 (MQ=255)                                                            gccgccCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCTACGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGcc                                                                                   >  1:89531/1‑149 (MQ=255)                                                            gccgccCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGc                                                                                    >  8:306173/1‑147 (MQ=255)                                                                 ccccGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCg                                                                             >  6:168539/1‑149 (MQ=255)                                                                    cgccgACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTg                                                                          <  3:110173/149‑1 (MQ=255)                                                                     gccgACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTgg                                                                         >  4:157027/1‑149 (MQ=255)                                                                          cTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAgg                                                                                >  7:83289/1‑137 (MQ=255)                                                                          cTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAgg                                                                                <  8:83289/137‑1 (MQ=255)                                                                              caacaCCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGcc                                                                <  5:168539/149‑1 (MQ=255)                                                                               aacaCCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGGCTACCAGTGTCTGTGCCAGGGCGTTGGCCTATTGc                                                                 >  6:14437/1‑147 (MQ=255)                                                                                 caccacCGCGTCCGGCTTCAGGCCGGCCT‑CGGACACTTAATCCACCAGGCACGCATGCCCAGACCAAAGCAGAGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTATGCAAGGGCGTTGTCCTCTTGcact                                                              >  4:66022/1‑145 (MQ=255)                                                                                   ccaccGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGGGTTGGCCTCTTGCCAgg                                                             >  6:67734/1‑147 (MQ=255)                                                                                       cgcgTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGc                                                       <  8:273053/149‑1 (MQ=255)                                                                                          gTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCtt                                                    >  3:171784/1‑149 (MQ=255)                                                                                            ccGTCTTCAGGCCGGCCT‑CGGCCACTTCTTGCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGc                                                  <  5:78267/149‑1 (MQ=255)                                                                                                     ggccggccT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCaccaccacc                                         >  4:132715/1‑149 (MQ=255)                                                                                                          gccT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGc                                    <  4:171784/149‑1 (MQ=255)                                                                                                                   cACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAgtgt                                                                                          <  7:162548/87‑1 (MQ=255)                                                                                                                   cACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAgtgt                                                                                          >  8:162548/1‑87 (MQ=255)                                                                                                                   cACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTACACCACCACCTTGGCGTCTTc                              >  6:244082/1‑147 (MQ=255)                                                                                                                     cTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGTGCGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGTACTGCCTTGCACCACCACCTTGGCGTCTTCCAg                           >  1:254397/1‑149 (MQ=255)                                                                                                                     cTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGTACTGCCTTGCACCACCACCTTGGCGTCTTCCa                            <  2:254397/147‑1 (MQ=255)                                                                                                                        cATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTg                        >  2:238053/1‑148 (MQ=255)                                                                                                                            caccaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTgcagca                   <  3:132715/149‑1 (MQ=255)                                                                                                                            caccaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTgcagc                    >  3:105568/1‑148 (MQ=255)                                                                                                                            caccaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGACTTGCACCACCAACTTGGAGTCTTCAAGGTgcagc                    >  4:45690/1‑148 (MQ=255)                                                                                                                             accaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCaa                  <  3:17927/149‑1 (MQ=255)                                                                                                                              ccaGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGGCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAAc                 <  3:45690/149‑1 (MQ=255)                                                                                                                                      gCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTg           <  7:306173/147‑1 (MQ=255)                                                                                                                                      gCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGCGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTg           <  1:296189/147‑1 (MQ=255)                                                                                                                                          gCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGt     >  6:182946/1‑149 (MQ=255)                                                                                                                                           cccGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGtt    <  3:157027/149‑1 (MQ=255)                                                                                                                                             cGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTACACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCg  >  1:237706/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              AGCACAGGCACATCGCCCCAGCCATTGCGCTTGAGCCCCTCCACCACGCCACTGAGCAAGCCGCCCCCGCCGACTGACAACACCACCGCGTCCGGCTTCAGGCCGGCCT‑CGGCCACTTCATCCACCAGGCTGGCATGCCCGGCCCACAGCAGTGGGTCGT‑CGAATGGGTGGATGAAGGCATCGTTCGGGCCTACCAGTGTCTGTGCCAGGGCGTTGGCCTCTTGCCAGGAACTGCCTTGCACCACCACCTTGGCGTCTTCCAGGTGCAGCAACTGCTTGGCCCGTTCG  >  NC_002947/3332554‑3332841

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