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%	137.4 / NA	37	K330K (AAA→AAG)	cysI	sulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/19);  total (18/19)

CAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGG  >  NC_000913/2889033‑2889186                                                                              |                                                                             cAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTg                                                                           <  2:519889/81‑1 (MQ=255)  aGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGa                                                                          >  1:519931/1‑81 (MQ=255)  aGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGa                                                                          >  1:519926/1‑81 (MQ=255)   gTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGAt                                                                         >  2:519941/1‑81 (MQ=255)       tcatcaATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATccc                                                                      <  2:519891/80‑1 (MQ=255)       tcatcaATGCCCTTAACCCAGCCAATACGATCGACTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATcccc                                                                     <  1:519890/81‑1 (MQ=255)          tcaATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcg                                                                  <  1:519892/81‑1 (MQ=255)           caATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcgc                                                                 >  1:519929/1‑81 (MQ=255)            aaTGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcgc                                                                 <  1:519893/80‑1 (MQ=255)             aTGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGa                                                               >  1:519945/1‑81 (MQ=255)              tGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGAc                                                              <  1:519898/81‑1 (MQ=255)              tGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGAc                                                              <  2:519896/81‑1 (MQ=255)              tGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGAc                                                              <  1:519897/81‑1 (MQ=255)              tGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGAc                                                              <  1:519895/81‑1 (MQ=255)              tGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGAc                                                              <  2:519894/81‑1 (MQ=255)                        ccAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTc                                                    <  1:519899/81‑1 (MQ=255)                         cAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTcc                                                   >  2:519943/1‑81 (MQ=255)                         cAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTcc                                                   >  1:519935/1‑81 (MQ=255)                           gCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGc                                                 >  1:519927/1‑81 (MQ=255)                            ccAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCt                                                >  1:519940/1‑81 (MQ=255)                              aaTACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCTGCGCGACGCTCCACTTCCGCttt                                              <  2:519900/81‑1 (MQ=255)                              aaTACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGAT‑CCCGCGCGACGCTCCACTTCCGCttt                                              >  1:519936/1‑80 (MQ=255)                                   gATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACg                                         >  2:519932/1‑81 (MQ=255)                                    aTCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGt                                        >  1:519933/1‑81 (MQ=255)                                     tCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtc                                       <  2:519901/81‑1 (MQ=255)                                       gCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                                      >  2:519942/1‑80 (MQ=255)                                        ccTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtctc                                     >  2:519947/1‑80 (MQ=255)                                            gTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAccc                                >  2:519934/1‑81 (MQ=255)                                               cGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAccccc                              <  2:519902/80‑1 (MQ=255)                                                ggTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAc                            >  1:519944/1‑81 (MQ=255)                                                ggTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAc                            >  2:519938/1‑81 (MQ=255)                                                  tGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAcgc                          >  2:519937/1‑81 (MQ=255)                                                    aaCTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACgcgc                        <  2:519903/81‑1 (MQ=255)                                                         atatGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCa                    <  2:519904/80‑1 (MQ=255)                                                          tatGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGc                  <  1:519905/81‑1 (MQ=255)                                                                  gATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAttt          <  2:519906/81‑1 (MQ=255)                                                                           aaaCTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTgg  <  2:519907/80‑1 (MQ=255)                                                                              |                                                                             CAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGG  >  NC_000913/2889033‑2889186

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