Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,790,908	G→A	8.1%	L214L (CTG→TTG)	sstT ←	serine/threonine transporter SstT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,790,908	0	G	A	8.1%	140.8 / 2.7	49	L214L (CTG→TTG)	sstT	serine/threonine transporter SstT
Reads supporting (aligned to +/- strand):  ref base G (24/21);  new base A (2/2);  total (26/23)
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.44e-01

CGAGCTGCGGGTGAAGAATGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAGGGTTACACCATTGGACAGGTCTTCGACCACGGTACGG  >  NC_002947/2790762‑2791053                                                                                                                                                    |                                                                                                                                                  cGAGCTGCGGGTGAAGAATGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCcag                                                                                                                                                   >  6:24588/1‑147 (MQ=255)            tGAAGAATGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAg                                                                                                                                      <  5:140715/149‑1 (MQ=255)               agaaTGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCg                                                                                                                                   <  6:155089/149‑1 (MQ=255)                gaaTGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCcag                                                                                                                                                   >  2:53286/1‑132 (MQ=255)                gaaTGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCcag                                                                                                                                                   <  1:53286/132‑1 (MQ=255)                      ggTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAgcg                                                                                                                            >  8:147216/1‑149 (MQ=255)                               gCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGc                                                                                                                   >  8:126837/1‑149 (MQ=255)                                     aCGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACt                                                                                                             >  7:223028/1‑149 (MQ=255)                                           gCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGc                                                                                                         >  1:126577/1‑147 (MQ=255)                                             acaGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAgg                                                                                                     >  7:624/1‑149 (MQ=255)                                                 cagcGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTgg                                                                                                 <  2:126577/149‑1 (MQ=255)                                                   gcGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGaa                                                                                               <  7:126837/149‑1 (MQ=255)                                                     gTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACt                                                                                             <  1:181921/149‑1 (MQ=255)                                                                    ttGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGAt                                                                              >  8:116908/1‑149 (MQ=255)                                                                              tCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGa                                                                                                >  7:104069/1‑121 (MQ=255)                                                                              tCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGa                                                                                                <  8:104069/121‑1 (MQ=255)                                                                                ttCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGc                                                                  >  5:95243/1‑149 (MQ=255)                                                                                 tCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGc                                                                  >  1:15809/1‑148 (MQ=255)                                                                                   cAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCGGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCg                                                                 >  2:246850/1‑147 (MQ=255)                                                                                   cAGAACACGATCAGCCGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCGGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGaa                                                               <  1:246850/149‑1 (MQ=255)                                                                                      aaCACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTgg                                                                                                 >  1:47923/1‑112 (MQ=255)                                                                                      aaCACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTgg                                                                                                 <  2:47923/112‑1 (MQ=255)                                                                                        cacGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGa                                                          <  2:15809/149‑1 (MQ=255)                                                                                         acGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGAt                                                         <  7:116908/149‑1 (MQ=255)                                                                                           gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCTACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       <  8:223028/149‑1 (MQ=255)                                                                                           gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAg                                                                                >  3:39913/1‑124 (MQ=255)                                                                                           gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAg                                                                                <  4:39913/124‑1 (MQ=255)                                                                                           gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  4:127209/1‑149 (MQ=255)                                                                                           gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  1:38688/1‑149 (MQ=255)                                                                                           gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  2:48300/1‑149 (MQ=255)                                                                                           gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  2:44408/1‑149 (MQ=255)                                                                                            aTCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGAc                                                      <  5:24588/149‑1 (MQ=255)                                                                                            aTCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGAc                                                      <  2:38688/149‑1 (MQ=255)                                                                                                 cGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCAcgc                                                 <  8:624/149‑1 (MQ=255)                                                                                                        atcaCCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGAt                                          <  3:127209/149‑1 (MQ=255)                                                                                                        atcaCCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCACGCTGGAACTGACCACGCAAAAGATGCCCAGCGGGGCGAACCCGATGACCACGCGTACGaa                                          >  3:23337/1‑148 (MQ=255)                                                                                                        atcaCCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGCTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGAt                                          >  5:49819/1‑149 (MQ=255)                                                                                                           aCCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCa                                        <  7:147216/148‑1 (MQ=255)                                                                                                             cAGCGCCACGAACAGCATGCAGCCGAGTCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCGAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAgg                                      >  1:7405/1‑149 (MQ=255)                                                                                                             cAGCGCCACGAACAGCATGCAGCCGAGTCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAgg                                      >  1:3159/1‑149 (MQ=255)                                                                                                             cAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGGCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAggg                                     <  4:69489/149‑1 (MQ=255)                                                                                                                cgcCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAGGGTTa                                  >  7:234931/1‑149 (MQ=255)                                                                                                                    aCGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAgg                                      >  3:40024/1‑141 (MQ=255)                                                                                                                    aCGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAgg                                      <  4:40024/141‑1 (MQ=255)                                                                                                                             aTGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCg                                                            <  5:55441/110‑1 (MQ=255)                                                                                                                             aTGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCg                                                            >  6:55441/1‑110 (MQ=255)                                                                                                                                   ccGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAGGGTTACACCATTGGACAGGTCTTc               <  1:48300/149‑1 (MQ=255)                                                                                                                                             aCCGCCAACAGGTGCAGGTAGCCAAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       <  7:249221/100‑1 (MQ=255)                                                                                                                                             aCCGCCAACAGGTGCAGGTAGCCAAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       <  7:245606/100‑1 (MQ=255)                                                                                                                                             aCCGCCAACAGGTGCAGGTAGCCAAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  8:245606/1‑100 (MQ=255)                                                                                                                                             aCCGCCAACAGGTGCAGGTAGCCAAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  8:249221/1‑100 (MQ=255)                                                                                                                                                 ccAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAGGGTTACACCATTGGACAGGTCTTCGACCACGGTACgg  >  6:233573/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGAGCTGCGGGTGAAGAATGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAGGGTTACACCATTGGACAGGTCTTCGACCACGGTACGG  >  NC_002947/2790762‑2791053

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