Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,890,403	T→C	Q499R (CAG→CGG)	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,403	0	T	C	100.0%	147.7 / NA	41	Q499R (CAG→CGG)	cysJ	sulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/21);  total (20/21)

TCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACC  >  NC_000913/2890272‑2890537                                                                                                                                    |                                                                                                                                       tCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAgg                                                                                                                                <  2:90873/140‑1 (MQ=255)  ccAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGa                                                                                                                               >  1:1170597/1‑140 (MQ=255)   cAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGaa                                                                                                                              <  2:574804/140‑1 (MQ=255)   cAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGaa                                                                                                                              <  2:355670/140‑1 (MQ=255)   cAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGaa                                                                                                                              <  1:876707/140‑1 (MQ=255)    aGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAg                                                                                                                             >  1:900075/1‑140 (MQ=255)     gcgcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTTATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGt                                                                                                                            >  1:475339/1‑140 (MQ=255)     gcgcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTa                                                                                                                                    >  1:179639/1‑132 (MQ=255)     gcgcCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTa                                                                                                                                    <  2:179639/132‑1 (MQ=255)                 cgcgCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGt                                                                                                                >  1:832167/1‑140 (MQ=255)                  gcgcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTg                                                                                                               <  2:1170597/140‑1 (MQ=255)                   cgcCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGc                                                                                                              >  1:1146836/1‑140 (MQ=255)                                 aGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGaa                                                                                                <  2:832167/140‑1 (MQ=255)                                     tGTTTTGTACGTAAACTTTTTCTTTTTGACCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGc                                                                                            <  2:820562/140‑1 (MQ=255)                                     tGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGc                                                                                            <  1:197183/140‑1 (MQ=255)                                             aCGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGttttt                                                                                    >  1:17472/1‑140 (MQ=255)                                             aCGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGttttt                                                                                    >  1:1030979/1‑140 (MQ=255)                                               gTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTAc                                                                                  <  1:806405/140‑1 (MQ=255)                                                   aCTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTgg                                                                              >  1:1129443/1‑140 (MQ=255)                                                   aCTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTgg                                                                              >  1:678249/1‑140 (MQ=255)                                                              ttGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTcgg                                                                   <  2:768087/140‑1 (MQ=255)                                                                   cgcgcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGgcggcg                                                              <  2:900075/140‑1 (MQ=255)                                                                   cgcgcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGgcggcg                                                              >  1:816630/1‑140 (MQ=255)                                                                     cgcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGcgcg                                                            >  1:658460/1‑140 (MQ=255)                                                                     cgcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGcgcg                                                            >  1:2351/1‑140 (MQ=255)                                                                      gcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGt                                                           <  2:816630/140‑1 (MQ=255)                                                                      gcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGt                                                           <  2:1146836/140‑1 (MQ=255)                                                                      gcgACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGt                                                           <  2:2351/140‑1 (MQ=255)                                                                        gACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTtg                                                         <  2:475339/140‑1 (MQ=255)                                                                                        tACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCgg                                         >  1:2953/1‑140 (MQ=255)                                                                                          cGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGaa                                       <  2:564847/140‑1 (MQ=255)                                                                                              tCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGgc                                   >  1:433607/1‑140 (MQ=255)                                                                                                aGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGgcgc                                 >  1:1127621/1‑140 (MQ=255)                                                                                                       cATCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCg                          >  1:918197/1‑140 (MQ=255)                                                                                                        aTCTTTGACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCgg                         >  1:1070372/1‑140 (MQ=255)                                                                                                              gACGTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTg                   >  2:1147542/1‑140 (MQ=255)                                                                                                                 gTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGc                <  2:658460/140‑1 (MQ=255)                                                                                                                 gTAGCGCTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGc                <  2:1127621/140‑1 (MQ=255)                                                                                                                     cgcTGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAAt            >  1:656723/1‑140 (MQ=255)                                                                                                                        tGCCACTCCACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCAt         <  2:433607/140‑1 (MQ=255)                                                                                                                               ccACCCGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCAcc  <  2:632407/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       TCCAGCGCCACAGCTCCGCGCCCTGTTCGCGCAGTTTGTCTTGTACGTAAACTTTTTCTTTTTGATCGCGCGACCAGGCAAGATCGATACGTGTCAGCACGCCATCTTTGACGTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACC  >  NC_000913/2890272‑2890537

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