Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,672,181	T→A	P45P (CCT→CCA)	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,181	0	T	A	100.0%	124.7 / NA	37	P45P (CCT→CCA)	hcaD	putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (19/18);  total (19/18)

GAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGT  >  NC_000913/2672049‑2672317                                                                                                                                     |                                                                                                                                         gaaagaaaAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCActct                                                                                                                                   <  2:623182/140‑1 (MQ=255) gaaagaaaAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCActct                                                                                                                                   >  1:768091/1‑140 (MQ=255)            cGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATg                                                                                                                        <  2:393287/140‑1 (MQ=255)               tCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTg                                                                                                                     <  1:788731/140‑1 (MQ=255)                       gTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGa                                                                                                             >  2:99960/1‑140 (MQ=255)                           gTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTcc                                                                                                         >  1:427131/1‑140 (MQ=255)                             ggCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTcccc                                                                                                       >  1:554081/1‑140 (MQ=255)                             ggCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTcccc                                                                                                       >  2:29777/1‑140 (MQ=255)                              gCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCg                                                                                                      >  2:122167/1‑140 (MQ=255)                                gggCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCa                                                                                                    <  2:400832/140‑1 (MQ=255)                                  gCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGt                                                                                                  <  2:427131/140‑1 (MQ=255)                                       cggcggCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTAcag                                                                                             <  1:99960/140‑1 (MQ=255)                                                tGGCTGCGGCCTCGCNACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTa                                                                                    >  1:599098/1‑140 (MQ=255)                                                         ccTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAAc                                                                           <  2:211849/140‑1 (MQ=255)                                                          cTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACt                                                                          <  2:141971/140‑1 (MQ=255)                                                                  gCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAg                                                                  <  1:29777/140‑1 (MQ=255)                                                                    cAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGa                                                                <  2:262175/140‑1 (MQ=255)                                                                       cAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGaaaa                                                             >  1:133892/1‑140 (MQ=255)                                                                               cACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTc                                                     >  1:257138/1‑140 (MQ=255)                                                                                aCCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCa                                                    >  1:342304/1‑140 (MQ=255)                                                                                    gTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTACTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTg                                                <  2:479267/140‑1 (MQ=255)                                                                                      gAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCa                                              <  2:554081/140‑1 (MQ=255)                                                                                        gCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCAtt                                            <  2:342304/140‑1 (MQ=255)                                                                                                tGTTTTCCGATGAGCGGCATCTTCCTTATGAACGACCACCGCTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTa                                    <  2:68549/140‑1 (MQ=255)                                                                                                            aGCGACATCTTCCTTATGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACg                        >  1:51692/1‑140 (MQ=255)                                                                                                                       ccTTGTGAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGa             >  1:417884/1‑140 (MQ=255)                                                                                                                             gAACGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACg       <  2:257138/140‑1 (MQ=255)                                                                                                                              aaCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTcg                <  2:764521/130‑1 (MQ=255)                                                                                                                              aaCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTcg                >  1:764521/1‑130 (MQ=255)                                                                                                                               aCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTcg                >  1:772611/1‑129 (MQ=255)                                                                                                                               aCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTcg                <  2:772611/129‑1 (MQ=255)                                                                                                                               aCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGacaca         <  1:691554/136‑1 (MQ=255)                                                                                                                               aCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGacaca         >  2:691554/1‑136 (MQ=255)                                                                                                                               aCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACgag     >  2:308486/1‑140 (MQ=255)                                                                                                                               aCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACgag     >  1:573293/1‑140 (MQ=255)                                                                                                                               aCGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACgag     >  2:16205/1‑140 (MQ=255)                                                                                                                                cGACCACCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACgaga    >  2:494309/1‑140 (MQ=255)                                                                                                                                  accaccACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGt  <  2:229202/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                         GAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAAAACAATGTTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGT  >  NC_000913/2672049‑2672317

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