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

ATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTG  >  NC_000913/2889038‑2889185                                                                         |                                                                            atcatcAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATccc                                                                     >  2:524183/1‑81 (MQ=255)     tcaATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcg                                                                 <  1:524116/81‑1 (MQ=255)      caATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcgc                                                                >  1:524178/1‑81 (MQ=255)       aaTGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcgc                                                                <  2:524117/80‑1 (MQ=255)       aaTGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCcgcgc                                                                >  1:524190/1‑80 (MQ=255)            ccTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACATGATCCCCGCGCGACGCt                                                          <  2:524118/81‑1 (MQ=255)             cTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGTACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTc                                                         >  2:524171/1‑81 (MQ=255)             cTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTc                                                         >  2:524188/1‑81 (MQ=255)             cTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTc                                                         >  1:524173/1‑81 (MQ=255)               tAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCa                                                       >  2:524180/1‑81 (MQ=255)                aaCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTcc                                                        <  2:524119/79‑1 (MQ=255)                 aCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCAc                                                      <  1:524120/80‑1 (MQ=255)                 aCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACt                                                     <  1:524121/81‑1 (MQ=255)                   ccAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTc                                                   <  2:524122/81‑1 (MQ=255)                     aGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCg                                                 <  1:524123/81‑1 (MQ=255)                        cAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCtt                                              <  1:524124/81‑1 (MQ=255)                         aaTACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCttt                                             <  2:524126/81‑1 (MQ=255)                         aaTACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCttt                                             <  1:524125/81‑1 (MQ=255)                         aaTACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCtt                                              >  1:524185/1‑80 (MQ=255)                           tACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTaa                                           <  1:524127/81‑1 (MQ=255)                             cGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTaaa                                          <  2:524128/80‑1 (MQ=255)                               aTCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGt                                       >  2:524191/1‑81 (MQ=255)                                  gCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtctc                                    <  2:524130/81‑1 (MQ=255)                                  gCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                                     >  1:524205/1‑80 (MQ=255)                                  gCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                                     <  1:524129/80‑1 (MQ=255)                                  gCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGtct                                     >  2:524182/1‑80 (MQ=255)                                       gTCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAccc                               >  2:524175/1‑81 (MQ=255)                                        tCCGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAcccc                              <  2:524131/81‑1 (MQ=255)                                         ccGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAAcccc                              >  1:524181/1‑80 (MQ=255)                                          cGGTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCa                            >  1:524176/1‑81 (MQ=255)                                           ggTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAc                           <  2:524132/81‑1 (MQ=255)                                           ggTGAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAc                           >  1:524189/1‑81 (MQ=255)                                              gAACTCATATGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCAcgcg                        >  1:524209/1‑81 (MQ=255)                                                     tatGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGc                 <  2:524133/81‑1 (MQ=255)                                                      atGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCg                >  2:524177/1‑81 (MQ=255)                                                      atGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCg                >  2:524186/1‑81 (MQ=255)                                                       tGGACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgt               >  1:524201/1‑81 (MQ=255)                                                        ggACGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCgtg              >  2:524198/1‑81 (MQ=255)                                                          aCGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTa            <  1:524134/81‑1 (MQ=255)                                                           cGGATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAt           >  2:524158/1‑81 (MQ=255)                                                            ggATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAtt          <  2:524135/81‑1 (MQ=255)                                                             gATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAttt         >  1:524193/1‑81 (MQ=255)                                                             gATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAttt         >  1:524192/1‑81 (MQ=255)                                                             gATCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTAttt         <  1:524136/81‑1 (MQ=255)                                                              aTCGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTg        <  2:524137/81‑1 (MQ=255)                                                                cGGTTCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGt      >  2:524204/1‑81 (MQ=255)                                                                    tCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGttttt  >  2:524207/1‑80 (MQ=255)                                                                    tCAAACTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTg  >  1:524187/1‑81 (MQ=255)                                                                         |                                                                            ATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTG  >  NC_000913/2889038‑2889185

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: ‑