Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,891,465	T→C	N145S (AAC→AGC)	cysJ ←	sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,891,465	0	T	C	100.0%	160.4 / NA	43	N145S (AAC→AGC)	cysJ	sulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/26);  total (17/26)

GGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGTTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGATTTGTTTGAATTTATAGTCGCCCGCGT  >  NC_000913/2891328‑2891609                                                                                                                                          |                                                                                                                                                 ggTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTg                                                                                                                                      <  1:103736/150‑1 (MQ=255)    aTTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGgcg                                                                                                                                   <  1:382902/150‑1 (MQ=255)      tCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCt                                                                                                                                 <  2:144438/150‑1 (MQ=255)        aaCATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCtttt                                                                                                                              >  1:386085/1‑151 (MQ=255)          cATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCtttttt                                                                                                                            <  2:597059/151‑1 (MQ=255)          cATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCtttttt                                                                                                                            <  2:512177/151‑1 (MQ=255)            tCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTgg                                                                                                                          >  2:54522/1‑151 (MQ=255)                    gACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAgg                                                                                                                  >  1:139470/1‑151 (MQ=255)                                cAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAgc                                                                                                      >  2:36473/1‑151 (MQ=255)                                   gCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGcgcg                                                                                                   <  1:36474/151‑1 (MQ=255)                                   gCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGcgcg                                                                                                   <  2:48535/151‑1 (MQ=255)                                    cGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGa                                                                                                  <  2:139474/151‑1 (MQ=255)                                       tCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACg                                                                                                <  1:329902/150‑1 (MQ=255)                                        caccacCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACgg                                                                                               <  2:249067/150‑1 (MQ=255)                                           caccCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGcttc                                                                                           >  1:3071/1‑151 (MQ=255)                                               cAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTcc                                                                                       >  1:245806/1‑151 (MQ=255)                                                 gTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCcgg                                                                                     <  2:245864/151‑1 (MQ=255)                                                    ccGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCcggcg                                                                                   <  1:567293/150‑1 (MQ=255)                                                         aGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGttccc                                                                             <  1:221962/151‑1 (MQ=255)                                                          gCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGttcccc                                                                            <  1:192436/151‑1 (MQ=255)                                                               cTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGtt                                                                                <  1:180422/142‑1 (MQ=255)                                                               cTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGtt                                                                                >  2:180423/1‑142 (MQ=255)                                                                   cGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTg                                                                   <  2:646609/151‑1 (MQ=255)                                                                                    tGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACg                                                   >  2:404952/1‑150 (MQ=255)                                                                                     ggCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGAt                                                 <  2:128991/151‑1 (MQ=255)                                                                                       cAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCGTTGGCGCTTTTTTGGAGAACAGGATCGTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGTCGTCACTACGa                                                  <  1:179032/148‑1 (MQ=255)                                                                                            aaaTTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGt                                          <  2:535736/151‑1 (MQ=255)                                                                                            aaaTTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGt                                          <  2:456607/151‑1 (MQ=255)                                                                                               ttCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGtttt                                       >  2:541564/1‑151 (MQ=255)                                                                                                 cATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGttttt                                      <  1:626895/150‑1 (MQ=255)                                                                                                       agCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGc                               >  1:31982/1‑151 (MQ=255)                                                                                                        gCTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCg                              <  2:324199/151‑1 (MQ=255)                                                                                                         cTATCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCg                              >  2:30486/1‑150 (MQ=255)                                                                                                          tATCGCCGAGGCTAAACACGGCAAACGCGGGGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGCAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATTGGCCGTTTTTCGCTGGc                               >  1:156211/1‑148 (MQ=255)                                                                                                           aTCGCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGAtt                           >  1:341103/1‑151 (MQ=255)                                                                                                              gCCGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGAtttgt                        >  1:543281/1‑151 (MQ=255)                                                                                                               ccGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGAtttgt                        <  1:647940/150‑1 (MQ=255)                                                                                                                cGAGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGAtttgttt                      >  1:467077/1‑151 (MQ=255)                                                                                                                  aGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGAtttgtttg                     >  1:651619/1‑150 (MQ=255)                                                                                                                  aGGCTAAACACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGATTTGTTTGa                    <  1:386711/151‑1 (MQ=255)                                                                                                                        aaCACGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGATTTGTTTGAATTtat              <  2:590707/151‑1 (MQ=255)                                                                                                                           acGGCAAACGCGGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGATTTGTTTGAATTTATAGt           <  2:245807/151‑1 (MQ=255)                                                                                                                                    gcgGTGCTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGATTTGTTTGAATTTATAGTCGCCCGCGt  >  1:469962/1‑151 (MQ=255)                                                                                                                                          |                                                                                                                                                 GGTATTCAACATCGGCATCGACACGGTCGAGCAGGCGTTCACCACCCAGTTCCGCCAGCTTGCTGTCGAAATCTTTCCCGGACTGGCAGAAAAATTCATAAGAGCTATCGCCGAGGCTAAACACGGCAAACGCGGTGTTTTCCAGCTTTGGCGCTTTTTTGGAGAACAGGAACTTATGCAGCGCGACGGCTTCTTCCGGCGGTTCCCCTTCCCCTTGCGTTGACGTCACTACGATGAGCAGTTTTTCGCTGGCGATTTGTTTGAATTTATAGTCGCCCGCGT  >  NC_000913/2891328‑2891609

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