Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,672,123	G→A	G26E (GGG→GAG)	hcaD →	phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,672,123	0	G	A	100.0%	176.1 / NA	51	G26E (GGG→GAG)	hcaD	phenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/22);  total (29/22)

GGTTTAGCGATTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAA  >  NC_000913/2671987‑2672256                                                                                                                                         |                                                                                                                                      ggTTTAGCGATTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGt                                                                                                                                     >  2:970156/1‑139 (MQ=255)            tCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAg                                                                                                                          >  2:572847/1‑139 (MQ=255)                     tCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTg                                                                                                                 >  2:725518/1‑139 (MQ=255)                                   ctgggctggATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCg                                                                                                   >  1:413282/1‑139 (MQ=255)                                                  aaGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGAAATCTTCCTTATGa                                                                                    <  2:274349/139‑1 (MQ=255)                                                     gAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACg                                                                                 >  2:894696/1‑139 (MQ=255)                                                       aacACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGAt                                                                                                        >  2:839446/3‑114 (MQ=255)                                                       aacACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGAt                                                                                                        <  1:839446/112‑1 (MQ=255)                                                       aGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGAc                                                                                                 <  1:1103584/121‑1 (MQ=255)                                                       aGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGAc                                                                                                 >  2:1103584/1‑121 (MQ=255)                                                       aGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGAc                                                                               >  1:654842/1‑139 (MQ=255)                                                       aGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGAc                                                                               >  1:220297/1‑139 (MQ=255)                                                        gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGc                                                                                                                         <  1:716925/96‑1 (MQ=255)                                                        gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGc                                                                                                                         >  2:716925/1‑96 (MQ=255)                                                        gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCttttaacggcct                                                                              >  2:753488/1‑129 (MQ=255)                                                        gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGAcc                                                                              >  2:986252/1‑139 (MQ=255)                                                        gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGAcc                                                                              >  2:808461/1‑139 (MQ=255)                                                        gCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGAcc                                                                              >  1:436490/1‑139 (MQ=255)                                                         cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCg                                                                                                          <  1:834295/110‑1 (MQ=255)                                                         cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCg                                                                                                          >  2:834295/1‑110 (MQ=255)                                                         cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGAc                                                                                                 <  1:126083/119‑1 (MQ=255)                                                         cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGAc                                                                                                 >  2:126083/1‑119 (MQ=255)                                                         cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCt                                                                             >  1:972239/1‑139 (MQ=255)                                                         cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCt                                                                             >  1:767626/1‑139 (MQ=255)                                                         cACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCt                                                                             >  2:637077/1‑139 (MQ=255)                                                               gaaagaaaAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCActc                                                                       >  1:389749/1‑139 (MQ=255)                                                                  agaaaAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCg                                                                    <  1:986252/139‑1 (MQ=255)                                                                    aaaaaaCGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGaa                                                                  <  2:389749/139‑1 (MQ=255)                                                                           gATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATg                                                           <  2:83514/139‑1 (MQ=255)                                                                           gATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATg                                                           <  1:725518/139‑1 (MQ=255)                                                                            aTCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGt                                                          <  2:1037909/139‑1 (MQ=255)                                                                                   ttGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGa                                                   <  1:894696/139‑1 (MQ=255)                                                                                        ggTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGa                                                                                    >  1:837692/1‑101 (MQ=255)                                                                                        ggTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGa                                                                                    <  2:837692/101‑1 (MQ=255)                                                                                        ggTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGAtt                                              >  2:843502/1‑139 (MQ=255)                                                                                        ggTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTGTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGAtt                                              >  1:1061264/1‑139 (MQ=255)                                                                                         gTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATg                                                           >  1:787215/1‑125 (MQ=255)                                                                                         gTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATg                                                           <  2:787215/125‑1 (MQ=255)                                                                                           ggCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTccc                                           >  1:526232/1‑139 (MQ=255)                                                                                           ggCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTccc                                           <  1:314902/139‑1 (MQ=255)                                                                                            gCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTcccc                                          <  2:413282/139‑1 (MQ=255)                                                                                             cGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCg                                         <  1:401156/139‑1 (MQ=255)                                                                                                gCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCGACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAg                                      >  1:713985/1‑139 (MQ=255)                                                                                                       gcggcAATGGCTGCGGCCTCGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTAcagc                               <  2:767626/139‑1 (MQ=255)                                                                                                                          cGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCg                                                                                                   <  2:295848/52‑1 (MQ=255)                                                                                                                          cGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCg                                                                                                   >  1:295848/1‑52 (MQ=255)                                                                                                                          cGCTACGCCAGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACtg            >  2:524569/1‑139 (MQ=255)                                                                                                                                   aGCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGa   <  2:139466/139‑1 (MQ=255)                                                                                                                                    gCAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGaa  <  2:654842/139‑1 (MQ=255)                                                                                                                                     cAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCt                                                                                         <  1:746971/51‑1 (MQ=255)                                                                                                                                     cAAGAGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCt                                                                                         >  2:746971/1‑51 (MQ=255)                                                                                                                                         |                                                                                                                                      GGTTTAGCGATTCGCGGCATTCGCCACGTAGCGGCTGGGCTGGATCTTTAAGGAAGCACGATGAAAGAAAAAACGATCATTATTGTCGGTGGCGGGCAAGCGGCGGCAATGGCTGCGGCCTCGCTACGCCAGCAAGGGTTCACCGGTGAGCTGCATCTGTTTTCCGATGAGCGACATCTTCCTTATGAACGACCTCCACTCTCGAAATCCATGTTGCTGGAAGATTCCCCGCAGTTACAGCAGGTGTTACCCGCTAACTGGTGGCAGGAA  >  NC_000913/2671987‑2672256

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