Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,672,278	G→A	G78S (GGT→AGT)	hcaD →	putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,672,278	0	G	A	100.0%	119.3 / NA	34	G78S (GGT→AGT)	hcaD	putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/21);  total (13/21)

CTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGAT  >  NC_000913/2672203‑2672352                                                                            |                                                                           cTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAgtgt                                                                        >  2:478960/1‑80 (MQ=255)   ggAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAAc                                                                     >  2:478942/1‑81 (MQ=255)    gAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAAc                                                                     <  2:478911/80‑1 (MQ=255)       gATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCAt                                                                  >  1:478967/1‑80 (MQ=255)       gATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATc                                                                 <  2:478912/81‑1 (MQ=255)            cccGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAAc                                                            <  1:478913/81‑1 (MQ=255)                cAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTg                                                        >  1:478952/1‑81 (MQ=255)                 aGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTgg                                                       >  1:478961/1‑81 (MQ=255)                    tACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTc                                                    >  2:478965/1‑81 (MQ=255)                         cagGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGac                                               <  2:478914/81‑1 (MQ=255)                           ggTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGacac                                             <  2:478915/81‑1 (MQ=255)                               ttACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACg                                          <  1:478916/80‑1 (MQ=255)                                tACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGacacac                                           >  1:478959/1‑78 (MQ=255)                                tACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGacacac                                           >  2:478973/1‑78 (MQ=255)                                  cccGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGagag                                      <  1:478917/81‑1 (MQ=255)                                     gCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGt                                     >  2:478962/1‑79 (MQ=255)                                     gCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGt                                     >  2:478958/1‑79 (MQ=255)                                      cTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAg                                  <  1:478918/81‑1 (MQ=255)                                         aCTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAgtgt                               <  1:478919/81‑1 (MQ=255)                                          cTGGTGGCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGtt                              <  2:478920/81‑1 (MQ=255)                                             gtggCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTa                             <  2:478921/79‑1 (MQ=255)                                                gCAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCaa                        <  1:478922/81‑1 (MQ=255)                                                 cAGGAAAACAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAAc                       <  2:478923/81‑1 (MQ=255)                                                   ggAAAACAATGTTCATCTGCATTCCAGTGTAGCCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACgg                     <  2:478924/81‑1 (MQ=255)                                                     aaaaCAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGc                    >  2:478966/1‑80 (MQ=255)                                                       aaCAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGaa                 <  2:478926/81‑1 (MQ=255)                                                       aaCAATGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGaa                 <  2:478925/81‑1 (MQ=255)                                                          aaTGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGc              <  1:478927/81‑1 (MQ=255)                                                           aTGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCt             <  1:478928/81‑1 (MQ=255)                                                           aTGTTCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCt             <  1:478929/81‑1 (MQ=255)                                                              ttCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACCCCCTCTCTTTAGTGTTAACCAACGGCGAAAGCTgg           >  1:1016202/1‑80 (MQ=255)                                                              ttCATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGc          <  2:478930/81‑1 (MQ=255)                                                                cATCTGCATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCAc        >  1:478969/1‑81 (MQ=255)                                                                      cATTCCAGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGAt  <  2:478931/81‑1 (MQ=255)                                                                            |                                                                           CTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGAT  >  NC_000913/2672203‑2672352

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