Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	107,336	G→C	10.4%	L401L (CTG→CTC)	PP_0101 →	sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	107,336	0	G	C	10.4%	96.7 / 2.3	37	L401L (CTG→CTC)	PP_0101	sulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (19/14);  new base C (2/2);  total (21/16)
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.85e-01

TCAAACTGGTCGACTTCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGATACCGTACCGGCGGGCACCACGCTG  >  NC_002947/107188‑107483                                                                                                                                                      |                                                                                                                                                    tCAAACTGGTCGACTTCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTg                                                                                                                                                     >  6:217794/1‑149 (MQ=255)         gTCGACTTCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCa                                                                                                                                              <  5:217794/148‑1 (MQ=255)         gTCGACTTCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTCGCGTTc                                                                                                                                               >  6:171275/1‑147 (MQ=255)         gTCGACTTCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTCGCGTTc                                                                                                                                               >  6:166523/1‑147 (MQ=255)            gACTTCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAggc                                                                                                                                          >  1:165979/1‑149 (MQ=255)              cTTCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCAACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGACGCTGGCGTTCAAggcgc                                                                                                                                        >  5:277966/1‑148 (MQ=255)               ttCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAggcggc                                                                                                                                       >  8:116826/1‑149 (MQ=255)                        ttCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGAGTTCAAGGCGGCGCGGCTg                                                                                                                                >  4:19180/1‑147 (MQ=255)                           cgcgGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGa                                                                                                                            <  5:120449/148‑1 (MQ=255)                                 cTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTCGCGTTCAAGGCGGCGCGGCTGAAGATCAACCt                                                                                                                     <  5:171275/149‑1 (MQ=255)                                 cTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTCGCGTTCAAGGCGGCGCGGCTGAAGATCAACCt                                                                                                                     <  5:166523/149‑1 (MQ=255)                                      ccGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGa                                                                                                                <  4:113094/149‑1 (MQ=255)                                                               aCCTTCGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGGTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATgg                                                                                        <  7:232571/148‑1 (MQ=255)                                                                    cgcgGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACt                                                                                    <  7:98770/147‑1 (MQ=255)                                                                             tGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTg                                                                         >  2:1543/1‑149 (MQ=255)                                                                               cgcTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTg                                                                         <  2:10774/147‑1 (MQ=255)                                                                                    gCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCgcggc                                                                  <  1:1543/149‑1 (MQ=255)                                                                                        ttATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGa                                                                <  4:30956/147‑1 (MQ=255)                                                                                        ttATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGa                                                                >  3:30956/1‑147 (MQ=255)                                                                                           tCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTc                                                           >  7:218604/1‑149 (MQ=255)                                                                                             tttACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCt                                                         >  5:43862/1‑149 (MQ=255)                                                                                             tttACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCt                                                         >  5:48726/1‑149 (MQ=255)                                                                                              ttACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCt                                                         >  5:93510/1‑148 (MQ=255)                                                                                                 ccGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATgg                                                                                                <  7:20116/106‑1 (MQ=255)                                                                                                 ccGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATgg                                                                                                >  8:20116/1‑106 (MQ=255)                                                                                                  cGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAgg                                                    >  7:86816/1‑149 (MQ=255)                                                                                                         ctgACCGGTGTACTGCTGGGCTTTGCCCTGACCTGTGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCgg                                              >  1:45973/1‑149 (MQ=255)                                                                                                         ctgACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCg                                               <  3:19180/147‑1 (MQ=255)                                                                                                                      tgctgGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTg                                >  5:174911/1‑149 (MQ=255)                                                                                                                        ctgGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCt                              <  6:93510/149‑1 (MQ=255)                                                                                                                        ctgGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCt                              >  4:3919/1‑149 (MQ=255)                                                                                                                                tGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGAtac                        <  6:43862/147‑1 (MQ=255)                                                                                                                                tGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGAtac                        <  6:48726/147‑1 (MQ=255)                                                                                                                                   ccTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGAtaccgtac                   >  5:135179/1‑149 (MQ=255)                                                                                                                                         cT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGATACCGTACcggcgg             <  2:165979/149‑1 (MQ=255)                                                                                                                                              cTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGATACCGTACCGGCGGGcac         >  7:160921/1‑149 (MQ=255)                                                                                                                                                gAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGATACCGTACCGGCGGGcacca       >  1:320135/1‑149 (MQ=255)                                                                                                                                                     tGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGATACCGTACCGGCGGGCACCACGCTg  >  4:137126/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TCAAACTGGTCGACTTCAAGGCGTTCCGCGGGCTGGGCCGTTATGGGCGCATGCCGATGTTTACCTACGCGGCGACTGCGCTGGCGATTATCTTTACCGACCTGCTGACCGGTGTACTGCTGGGCTTTGCCCTGACCT‑TGCTGAAGCTGGCGTTCAAGGCGGCGCGGCTGAAGATCAACCTGGTGAGCCTGGCCAAGGATGGGCACATGGAACTTCGCCTCAGTGGCGCGGCGACCTTCCTGAAGGTGCCGGCGCTGACCCAGGTGCTGGATACCGTACCGGCGGGCACCACGCTG  >  NC_002947/107188‑107483

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑