Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,670,928	T→C	C7C (TGT→TGC)	hcaC →	putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,670,928	0	T	C	100.0%	192.0 / NA	56	C7C (TGT→TGC)	hcaC	putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/28);  total (28/28)

GGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTC  >  NC_000913/2670792‑2671056                                                                                                                                         |                                                                                                                                 ggAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCGGAGAGTAATGTTCTGAAGAATAGAATTTATGCGTGcccc                                                                                                                               >  2:298468/1‑140 (MQ=255)  gAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCg                                                                                                                              >  2:439801/1‑140 (MQ=255)  gAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTCATGAATCGAATTTATGCGTGCCCCg                                                                                                                              >  1:334057/1‑140 (MQ=255)   aaCGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGt                                                                                                                             >  2:160096/1‑140 (MQ=255)       cgTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCg                                                                                                                         >  1:568947/1‑140 (MQ=255)        gTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCgg                                                                                                                        <  2:625141/140‑1 (MQ=255)                aaaGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCg                                                                                                                <  2:608342/140‑1 (MQ=255)                    ttCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCTGTTGCGGATGTGCCGGAgg                                                                                                            >  2:187445/1‑140 (MQ=255)                      tgCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTTCCCCGTTGCGGATGTGCCGGAgggg                                                                                                          >  1:349509/2‑139 (MQ=255)                         cGTCTGGAAGATGACAACTGGCGCTTGCTTGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTAAGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAg                                                                                                       <  1:387859/140‑1 (MQ=255)                         cGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAg                                                                                                       <  1:132804/140‑1 (MQ=255)                         cGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGCAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAg                                                                                                       <  1:835621/140‑1 (MQ=255)                               gAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCtctc                                                                                                 <  2:710735/140‑1 (MQ=255)                                 aGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCg                                                                                               <  1:160096/140‑1 (MQ=255)                                 aGATGACAACTGGCGCTTGCTGCAACGGGATCTCGACCTAGATCAAGCCGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCg                                                                                               <  2:246957/140‑1 (MQ=255)                                             gcgcTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCtctc                                                                                                 >  1:737840/1‑126 (MQ=255)                                             gcgcTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCtctc                                                                                                 <  2:737840/126‑1 (MQ=255)                                                cTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCTGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGcc                                                                                <  1:187445/140‑1 (MQ=255)                                                cTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGcc                                                                                <  1:491920/140‑1 (MQ=255)                                                     gagacagggATATCGTCCTGGACCATGCGGTAATCACATCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCa                                                                           <  1:273263/134‑1 (MQ=255)                                                          cGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGcc                                                                      <  1:439801/140‑1 (MQ=255)                                                             gATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGCGGCTCTCTGGAGCGTTACCTCGCCCGTCATCGCCCTg                                                                   >  1:263112/1‑140 (MQ=255)                                                                       tGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTGCTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTg                                                         <  1:298468/140‑1 (MQ=255)                                                                       tGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTTAGGCTCTCCGGAGCGAGACCTCGCCGGTCATCGCCCTGTTcctgtctc                                                         >  2:273263/1‑133 (MQ=255)                                                                       tGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTg                                                         >  1:787593/1‑140 (MQ=255)                                                                              aGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATTCCTCGCCCGTCTTCACCCTGTTCAACGTTGGCGGCGa                                                  >  1:413525/1‑140 (MQ=255)                                                                              aGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGa                                                  >  2:606151/1‑140 (MQ=255)                                                                              aGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGa                                                  >  1:786139/1‑140 (MQ=255)                                                                              aGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTGCAACGTTGGCGGCGa                                                  >  1:330915/1‑140 (MQ=255)                                                                               gCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAg                                                 >  2:129015/1‑140 (MQ=255)                                                                                             ccATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTaa                                   >  1:213172/1‑140 (MQ=255)                                                                                               aTAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGAGAGGCTCTCCGGATCGATACCTCGCCCGACATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACg                                 >  2:679274/1‑140 (MQ=255)                                                                                                    cTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTAt                                           <  2:851321/125‑1 (MQ=255)                                                                                                    cTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTAt                                           >  1:851321/1‑125 (MQ=255)                                                                                                       aGTGTACTGTTCTGATGCATCGTATTTATGCTTGCCCCGATGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGc                         <  2:263112/140‑1 (MQ=255)                                                                                                       aGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGc                         >  1:103314/1‑140 (MQ=255)                                                                                                       aGTGTACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGc                         >  2:610424/1‑140 (MQ=255)                                                                                                       aGTGTACTGTTCTGACGAATCGAATCTATCCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGc                         <  1:241543/140‑1 (MQ=255)                                                                                                       aGTGAACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGc                         >  2:725379/1‑140 (MQ=255)                                                                                                         tgtACTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAg                       <  2:568947/140‑1 (MQ=255)                                                                                                            aCTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCa                    >  1:533139/1‑140 (MQ=255)                                                                                                             cTGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCAt                   >  1:846052/1‑140 (MQ=255)                                                                                                              tGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCa                    >  2:493434/1‑138 (MQ=255)                                                                                                              tGTTCTGATGAATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCa                    <  1:493434/138‑1 (MQ=255)                                                                                                                       gaATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGt         <  1:606151/140‑1 (MQ=255)                                                                                                                       gaATCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGt         <  1:610424/140‑1 (MQ=255)                                                                                                                        aaTCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTc        <  2:786139/140‑1 (MQ=255)                                                                                                                        aaTCGAATCTATGCGTTCCCCGTTGCGGATTTGCCAGAGAGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTc        <  2:413525/140‑1 (MQ=255)                                                                                                                          tCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCaa      >  2:589428/1‑140 (MQ=255)                                                                                                                          tCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCaa      <  1:589428/140‑1 (MQ=255)                                                                                                                          tCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCaa      >  1:64791/1‑140 (MQ=255)                                                                                                                          tCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCaa      <  1:725379/140‑1 (MQ=255)                                                                                                                          tCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCaa      <  1:129015/140‑1 (MQ=255)                                                                                                                          tCGAATTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCaa      >  1:740889/1‑140 (MQ=255)                                                                                                                          tCGAATTTATGCGTGCCACGTTTCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCaa      <  2:349509/140‑1 (MQ=255)                                                                                                                             aaTTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCGATTAACGCTCGTTGCAGCCATGGTAATGCGTCAATGt   >  1:283242/1‑140 (MQ=255)                                                                                                                              aTTTATGCGTGCCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTc  <  2:64791/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                 GGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGTCCCGTTGCGGATGTGCCGGAGGGTGAGGCTCTCCGGATCGATACCTCGCCCGTCATCGCCCTGTTCAACGTTGGCGGCGAGTTTTATGCCATTAACGATCGTTGCAGCCATGGTAATGCGTCAATGTC  >  NC_000913/2670792‑2671056

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