Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,790,924	G→A	8.5%	L208L (CTC→CTT)	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,924	0	G	A	8.5%	134.1 / 2.5	47	L208L (CTC→CTT)	sstT	serine/threonine transporter SstT
Reads supporting (aligned to +/- strand):  ref base G (24/19);  new base A (2/2);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

AGAATGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAGGGTTACACCATTGGACAGGTCTTCGACCACGGTACGG  >  NC_002947/2790776‑2791053                                                                                                                                                      |                                                                                                                                  agaaTGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCg                                                                                                                                   <  6:155089/149‑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                                                                                                 <  2:47923/112‑1 (MQ=255)                                                                        aaCACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTgg                                                                                                 >  1:47923/1‑112 (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                                                       >  2:48300/1‑149 (MQ=255)                                                                             gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  2:44408/1‑149 (MQ=255)                                                                             gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  4:127209/1‑149 (MQ=255)                                                                             gATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  1:38688/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                                                       >  8:249221/1‑100 (MQ=255)                                                                                                                               aCCGCCAACAGGTGCAGGTAGCCAAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       <  7:249221/100‑1 (MQ=255)                                                                                                                               aCCGCCAACAGGTGCAGGTAGCCAAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       >  8:245606/1‑100 (MQ=255)                                                                                                                               aCCGCCAACAGGTGCAGGTAGCCAAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGa                                                       <  7:245606/100‑1 (MQ=255)                                                                                                                                   ccAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAGGGTTACACCATTGGACAGGTCTTCGACCACGGTACgg  >  6:233573/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                  AGAATGCGGTAATGCCGCTTTCACGCAGGCACAGCAGCGTCAGCGGGTAGGGGTTGCGGCGGATCTTCCAGAACACGATCAGCGGGTTCATCACCAGCGCCACGAACAGCATGCAGCCGA‑TCAGTACCGCCAGCAGGTGCAGGTAGCCGAGCAGCGCGTCCAGGCCCGACTGGGCCAGGGTGGAACTGACCAGGCCAAAGATGCCCAGCGGGGCGAAGCGGATGACCACGCGTACGATCAGGGTTACACCATTGGACAGGTCTTCGACCACGGTACGG  >  NC_002947/2790776‑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: ‑