Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,889,110	T→C	K330K (AAA→AAG)	cysI ←	sulfite reductase, hemoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,889,110	0	T	C	100.0%	168.0 / NA	46	K330K (AAA→AAG)	cysI	sulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/24);  total (22/24)

CCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGG  >  NC_000913/2889032‑2889186                                                                               |                                                                             ccAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAActt                                                                            <  1:521458/81‑1 (MQ=255)  cAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAActt                                                                            >  2:521519/1‑80 (MQ=255)   aGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGa                                                                          <  1:521459/81‑1 (MQ=255)     ttATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATc                                                                        >  2:521512/1‑81 (MQ=255)     ttATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATc                                                                        <  2:521460/81‑1 (MQ=255)     ttATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGGCCGGGGAACTCATATGGACGGATCGGTTCAAACTTGATc                                                                        <  1:521461/81‑1 (MQ=255)       atcatcAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATcc                                                                       <  1:521462/80‑1 (MQ=255)        tcatcaATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATcccc                                                                     <  2:521463/81‑1 (MQ=255)           tcaATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcg                                                                  >  1:521501/1‑81 (MQ=255)             aaTGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcgcg                                                                <  2:521464/81‑1 (MQ=255)             aaTGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcgc                                                                 >  1:521509/1‑80 (MQ=255)             aaTGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcgc                                                                 >  1:521513/1‑80 (MQ=255)                 cccTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGc                                                            <  1:521465/81‑1 (MQ=255)                  ccTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCt                                                           >  2:521508/1‑81 (MQ=255)                     tAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCa                                                        >  1:521523/1‑81 (MQ=255)                        cccAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACtt                                                     >  2:521537/1‑81 (MQ=255)                          cAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTcc                                                   <  1:521466/81‑1 (MQ=255)                          cAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTcc                                                   <  2:521467/81‑1 (MQ=255)                           aGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCg                                                  >  1:521521/1‑81 (MQ=255)                            gCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGc                                                 <  1:521468/81‑1 (MQ=255)                              cAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCtt                                               <  2:521469/81‑1 (MQ=255)                               aaTACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCtt                                               >  2:521522/1‑80 (MQ=255)                                atgcGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTa                                             <  2:521471/78‑1 (MQ=255)                                aTACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCttt                                              >  1:521526/1‑80 (MQ=255)                                aTACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTa                                             <  1:521470/81‑1 (MQ=255)                                    gATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACg                                         >  1:521511/1‑81 (MQ=255)                                          cTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCaa                                   <  2:521472/81‑1 (MQ=255)                                              tCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAcccc                               <  1:521473/81‑1 (MQ=255)                                              tCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAcccc                               <  2:521474/81‑1 (MQ=255)                                                cGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAccccc                              <  1:521475/80‑1 (MQ=255)                                                 ggTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAc                            <  1:521476/81‑1 (MQ=255)                                                   tGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAcgc                          >  2:521514/1‑81 (MQ=255)                                                     aaCTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAcgcg                         <  2:521477/80‑1 (MQ=255)                                                      aCTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCt                       <  2:521478/81‑1 (MQ=255)                                                          atatGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTcc                     <  2:521479/79‑1 (MQ=255)                                                             tGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgt                <  2:521480/81‑1 (MQ=255)                                                              ggACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgtg               >  2:521533/1‑81 (MQ=255)                                                               gACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgtgt              <  1:521481/81‑1 (MQ=255)                                                                aCGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTa             >  1:521518/1‑81 (MQ=255)                                                                 cGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAt            <  2:521482/81‑1 (MQ=255)                                                                  ggATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAtt           >  1:521517/1‑81 (MQ=255)                                                                  ggATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAtt           >  1:521531/1‑81 (MQ=255)                                                                      cGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGt       >  1:521530/1‑81 (MQ=255)                                                                       ggTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGtt      >  1:521534/1‑81 (MQ=255)                                                                           cAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTgg  >  1:521528/1‑81 (MQ=255)                                                                           cAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTgg  >  1:521538/1‑81 (MQ=255)                                                                           cAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTgg  >  1:521525/1‑81 (MQ=255)                                                                           cAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTgg  >  2:521527/1‑81 (MQ=255)                                                                               |                                                                             CCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGG  >  NC_000913/2889032‑2889186

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