Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,890,626	G→A	L425F (CTC→TTC)	cysJ ←	sulfite reductase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,890,626	0	G	A	93.8%	80.8 / ‑6.5	33	L425F (CTC→TTC)	cysJ	sulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (22/8);  total (25/8)
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

GTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACC  >  NC_000913/2890384‑2890867                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    gTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCCATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAAAACGCGGACTTCGCCCTCTTCTTACACGCGGTCAGGGAGGAAGATgg                                                                                                                                                                                                                                             >  1:335652/1‑251 (MQ=255)    gcgcTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAgg                                                                                                                                                                                                                                          <  1:301283/251‑1 (MQ=255)                           ggAAGTCTTCCGTAAAGTGCGGCTTACCAAAGAACAGCCCGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGGTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGAGAAGAAGCTGGAGGCACAACCGGCACGGGCGCGGCCt                                                                                                                                                                                                                   >  1:217427/1‑251 (MQ=255)                                         aaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACg                                                                                                                                                                                                    <  1:67753/251‑1 (MQ=255)                                            tGCGGAATAACAAAGAACAGCAAGTTTTTACCTGGCGATTCGTCGGCGGCGAGTTGCAGCATACAGACGCGGTACGGCGCAATAACGGTGCATGGGCCAATCATAATCACCGGGGTTTCTGGATTGGATGGCAGGCGGGAGTTATCGATATGTTCGGAAAGTACTAGGACGTAGACCTCCTCTTCAACACGGTCAGAGAGGAAGCAAGAGGCAACACCGGCCAGGGGGCGGCCTTCCAA‑GTAGTAACgagc                                                                                                                                                                                                 >  2:14220/1‑248 (MQ=255)                                                 aTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGcaccacac                                                                                                                                                                                            >  1:132330/1‑251 (MQ=255)                                                      cAAAGAACAGACAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGACAGGCGCAATACCGGTGCCTGGGCCAATCATAATCCCCGGGGTTTCTGGTTAGGCTGACCGGCGGACGTTATCGGTATGTTCGCTAAATGCGCAGACTTCGCCCTCTTCTTCCACGCGGTCGGAGACGAACCGGGAGGCACCACAGGCCCGTGCGCTGACTTCCAc‑gacgcatcgccccaccccgacc                                                                                                                                                                                       >  2:255129/1‑230 (MQ=255)                                                                       ttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGCTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGCAACGCACCACACCAACGGTGACGTGTACTTCG‑tt                                                                                                                                                                     >  1:88041/1‑251 (MQ=255)                                                                             ggCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGCGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTCCTTCG‑TTCTCGAc                                                                                                                                                               >  1:348603/1‑251 (MQ=255)                                                                               cgcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGGCAGCGAAGAAGCTGGAGGCACCACCGGCACGGACGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑Ttctagcctt                                                                                                                                                             >  1:68073/1‑245 (MQ=255)                                                                                     gtcgGCGGCGCGTTGCTGCATAGAGGCGCGGAACGGCGCAATACCGGCGCCTGGGCCAATCATAATCACCGGGGTTTCTGGCTTGGCTGGCCGGCGGAAGTTATCGTTATGTTCGATACATACGCGGACTTCGCCCTCCTCTTCCACGCGGCCAGAGCAGAAGCTGGAGGCACCACCGGCACGGGCGTGTCATTCCAC‑GGCGTAACACACCACACCGACGGGGGAGCGTACTTCG‑TTTTCGCCTTCCGcac                                                                                                                                                       >  2:163287/1‑249 (MQ=255)                                                                                            gcgcGTTGCTGAATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAAAACGCGGCCTTCGCCCTCTTCTTCCACGCGATCAGCAAAGAAGCTGGAGGCCCCAACGGCACGGGCGCTGGCTACCACAGT‑GTAACGCAACACACCCACGGCGGCGTGTACTTCG‑TTTTCGCCTTCCGCCttcaaggg                                                                                                                                                >  2:74914/1‑244 (MQ=255)                                                                                                     tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATgg                                                                                                                                        >  1:273814/1‑250 (MQ=255)                                                                                                      gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGa                                                                                                                                       <  1:327043/250‑1 (MQ=255)                                                                                                             ggCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCgaagaa                                                                                                                                                                                                                                                  >  1:182247/1‑138 (MQ=255)                                                                                                             ggCGCGGAACGGCGCAAGACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATGCGCGGACGTCGCCCTCTTCTTCCACGCGGTCAGCgaagaa                                                                                                                                                                                                                                                  <  2:182247/138‑1 (MQ=255)                                                                                                                   gAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTAGATAACTACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGTCGCGGCCTTCCAC‑GGCGTAACGCACCCCACCAACGGTGACGTGTACTTCG‑TTCTCTCCTTCCGCCTGCGAGGAGGCGAGGGAATAcgaccgcggcg                                                                                                                         >  2:380980/1‑241 (MQ=255)                                                                                                                       tgcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACg                                                                                                                                                                                       >  1:280763/2‑186 (MQ=255)                                                                                                                                                              gTTTCTGGATTGACTGGCAGGCGGAAGTTATAGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGCAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTACCGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGAAGTCCTCCTGCAAGGAGGCGATGTAATCCATACGAGGCGTCCGCGGTGGCAGCAGATTAATTAGTGCTTCGGCATCacgcgt                                                                              >  2:38821/1‑246 (MQ=255)                                                                                                                                                                     gATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCggg                                                                       <  1:187300/251‑1 (MQ=255)                                                                                                                                                                      aTTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCCCGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACgg                                                                                                                                                                                      >  2:85450/1‑140 (MQ=255)                                                                                                                                                                      aTTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACgg                                                                                                                                                                                      <  1:85450/140‑1 (MQ=255)                                                                                                                                                                      aTTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAAGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCgggg                                                                      >  1:300676/1‑251 (MQ=255)                                                                                                                                                                             ggcAAGCGGAAGTTATCGTTATGTTCGCTAACTACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACAAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGCGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCCGATTCATTAGCGCTTCGGCCTCAAGCCGTGACGGGgagcaccc                                                               >  2:158763/1‑246 (MQ=255)                                                                                                                                                                                aGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCAc                                                            >  1:382009/1‑251 (MQ=255)                                                                                                                                                                                     gAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGt                                                       >  1:119840/1‑251 (MQ=255)                                                                                                                                                                                                 aTGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGcg                                           >  1:165722/1‑251 (MQ=255)                                                                                                                                                                                                             tacccggtcTTCGCCCTCGTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTCCTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCgcgg                                  <  2:113154/240‑1 (MQ=255)                                                                                                                                                                                                             tACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCgcgg                                  >  1:113154/1‑248 (MQ=255)                                                                                                                                                                                                                gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGat                                                                                                       <  1:158312/176‑1 (MQ=255)                                                                                                                                                                                                                gcgGACTTCGCCCTCTTCTTCCACGCGGTAAGCGAAGAAGCTGGAGGAACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAAAGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGCCTTCCGCCTGCGCGGCGGCGATGGCATACAGACGCGGCGTCAGCGGGCGccgcagct                                                                                                       >  2:158312/1‑174 (MQ=255)                                                                                                                                                                                                                    aCTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTg                        >  1:250041/1‑251 (MQ=255)                                                                                                                                                                                                                         gCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACt                   >  1:254597/1‑251 (MQ=255)                                                                                                                                                                                                                          cccTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACtt                  >  1:270428/1‑251 (MQ=255)                                                                                                                                                                                                                                          ggTCAGGGAAGAAGCGGGAGGCCCCCCCGGCACGGGAGCGGCCTTCCAC‑GTCGGAACGCCCCCCTCCAACGGTGACGTGTACTTCGGTTC‑CGACTTCCGCCTGCGGGGAGGCGATGGAAGACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTGCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCccacc  <  2:270428/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    GTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAC‑GTCGTAACGCACCACACCAACGGTGACGTGTACTTCG‑TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACC  >  NC_000913/2890384‑2890867

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