Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,531,596	(G)6→7	coding (134/762 nt)	cysZ →	sulfate transporter, sulfite inhibited

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,531,590	1	.	G	94.1%	179.6 / 3.6	51	coding (128/762 nt)	cysZ	sulfate transporter, sulfite inhibited
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (19/29);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑T‑GGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTG  >  NC_000913/2531453‑2531723                                                                                                                                            |                                                                                                                                        cATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑T‑g                                                                                                                                        <  1:267093/139‑1 (MQ=255)      tGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑T‑ggggg‑g                                                                                                                                  <  1:1011571/139‑1 (MQ=255)      tGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑T‑ggggg‑g                                                                                                                                  <  1:1098072/139‑1 (MQ=255)         ttATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑T‑ggggg‑ggcg                                                                                                                               <  1:900504/139‑4 (MQ=255)                 tcatcaTTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtgg                                                                                                                       <  1:523594/139‑1 (MQ=255)                  catcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtggt                                                                                                                      >  2:729111/1‑139 (MQ=255)                  catcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtggt                                                                                                                      >  1:34228/1‑139 (MQ=255)                   atcatTCACATCTGCCGCACGCAGCTGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGAAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAACATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtggtg                                                                                                                     <  1:550265/139‑1 (MQ=255)                   atcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtggtg                                                                                                                     <  2:1000357/139‑1 (MQ=255)                   atcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtggtg                                                                                                                     <  2:685329/139‑1 (MQ=255)                   atcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtggtg                                                                                                                     <  2:446113/139‑1 (MQ=255)                   atcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtggtg                                                                                                                     <  1:729111/139‑1 (MQ=255)                   atcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCtggtg                                                                                                                     <  1:552058/139‑1 (MQ=255)                        tCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGctct                                                                                                                >  2:331195/1‑139 (MQ=255)                         cacaTCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCtt                                                                                                               >  2:384228/1‑139 (MQ=255)                                 ccccACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGc                                                                                                       >  2:3295/1‑139 (MQ=255)                                  cccACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCt                                                                                                      <  2:1010329/139‑1 (MQ=255)                                      cGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGAt                                                                                                  >  1:421196/1‑139 (MQ=255)                                         aGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTc                                                                                               >  1:913300/1‑139 (MQ=255)                                                tttACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATcc                                                                                        <  2:32619/139‑1 (MQ=255)                                                 ttACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATccc                                                                                       <  2:860028/139‑1 (MQ=255)                                                   aCTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGa                                                                                     <  1:600343/139‑1 (MQ=255)                                                     tATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGAct                                                                                   <  1:266920/139‑1 (MQ=255)                                                       ttttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGActct                                                                                 >  2:648348/1‑139 (MQ=255)                                                        tttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACtctc                                                                                <  2:421196/139‑1 (MQ=255)                                                          tGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCAt                                                                              <  1:575485/139‑1 (MQ=255)                                                          tGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCAt                                                                              >  1:1092812/1‑139 (MQ=255)                                                                 ggCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTAc                                                                       >  1:86999/1‑139 (MQ=255)                                                                        aGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCgg                                                                <  2:106964/139‑1 (MQ=255)                                                                         gCTCGTCTCGCAACCTGGGATACGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑T‑GGGGGCG‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGa                                                               <  1:648348/139‑1 (MQ=255)                                                                                        tGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt                                                <  1:331195/139‑1 (MQ=255)                                                                                         gggATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTg                                               >  2:388794/1‑139 (MQ=255)                                                                                          ggATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGa                                              <  1:336946/139‑1 (MQ=255)                                                                                              tCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTACCGGACTGGCTACAATGGCTGAGTTa                                          >  1:54212/1‑139 (MQ=255)                                                                                                ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGACG‑GGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATc                                        >  2:566728/1‑139 (MQ=255)                                                                                                ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGa                                                               >  1:1033050/1‑116 (MQ=255)                                                                                                ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGa                                                               <  2:1033050/116‑1 (MQ=255)                                                                                                ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATc                                        >  1:114097/1‑139 (MQ=255)                                                                                                 gCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCt                                       >  1:914193/1‑139 (MQ=255)                                                                                                             tttACCGCTGCTGGTCAATATTTTTTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCctct                           <  1:712267/139‑1 (MQ=255)                                                                                                              ttACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTg                          <  1:888872/139‑1 (MQ=255)                                                                                                               tACCGCTGCTGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTgg                         <  1:388794/139‑1 (MQ=255)                                                                                                                   gctgctGGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGt                     <  2:387826/139‑1 (MQ=255)                                                                                                                      gctgGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGa                                                                                     <  2:1071604/72‑1 (MQ=255)                                                                                                                      gctgGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGa                                                                                     >  1:1071604/1‑72 (MQ=255)                                                                                                                      gctgGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCAt                  >  2:159507/1‑139 (MQ=255)                                                                                                                       ctgGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGTCTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTAt                                         <  1:782885/115‑1 (MQ=255)                                                                                                                       ctgGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTAt                                         >  2:782885/1‑115 (MQ=255)                                                                                                                        tgGTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCAtct                <  1:3295/139‑1 (MQ=255)                                                                                                                         ggTCAATATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCAtctc               <  1:55376/139‑1 (MQ=255)                                                                                                                           tCAATCTTTTGTTGA‑T‑GGGGG‑GACGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGGTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtg             <  2:806374/139‑1 (MQ=255)                                                                                                                             aaTATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtgtg           <  1:566728/139‑1 (MQ=255)                                                                                                                             aaTATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtgtg           <  1:384228/139‑1 (MQ=255)                                                                                                                             aaTATTTTGTTGA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtgtg           <  1:159507/139‑1 (MQ=255)                                                                                                                                      ttgA‑TGGGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAgtg  >  2:588094/1‑139 (MQ=255)                                                                                                                                            |                                                                                                                                        CATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGA‑T‑GGGGG‑G‑CGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTG  >  NC_000913/2531453‑2531723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑