Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,991,148	A→G	intergenic (‑382/‑5)	hemC ← / → cyaA	hydroxymethylbilane synthase/adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,148	0	A	G	100.0%	260.6 / NA	73	intergenic (‑382/‑5)	hemC/cyaA	hydroxymethylbilane synthase/adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/42);  total (31/42)

TTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCA  >  NC_000913/3991012‑3991280                                                                                                                                         |                                                                                                                                     ttttCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCg                                                                                                                                    <  2:733716/139‑1 (MQ=255)      gtcgtGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTg                                                                                                                               <  2:220322/139‑1 (MQ=255)      gtcgtGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTg                                                                                                                               <  2:373627/139‑1 (MQ=255)            aaaCTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACctc                                                                                                                         <  2:108569/139‑1 (MQ=255)            aaaCTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACctc                                                                                                                         <  1:1314540/139‑1 (MQ=255)                    aaaaCCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTga                                                                                                                 <  2:562701/139‑1 (MQ=255)                      aaCCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTgaga                                                                                                               <  1:345958/139‑1 (MQ=255)                          aGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGActct                                                                                                           <  1:287212/139‑1 (MQ=255)                          aGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGActct                                                                                                           >  1:753493/1‑139 (MQ=255)                                      tGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTAcc                                                                                                                           >  1:502592/1‑111 (MQ=255)                                      tGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTAcc                                                                                                                           <  2:502592/111‑1 (MQ=255)                                              gTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCa                                                                                       <  1:224854/139‑1 (MQ=255)                                              gTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCa                                                                                       <  2:419246/139‑1 (MQ=255)                                                 aCCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCg                                                                                                                                    >  2:317722/1‑91 (MQ=255)                                                       gACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAAt                                                                              <  2:753493/139‑1 (MQ=255)                                                       gACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAAt                                                                              >  2:1775774/1‑139 (MQ=255)                                                       gACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCATCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAAt                                                                              <  2:1391543/139‑1 (MQ=255)                                                            atttcgtatCCCGAATGCCGCGTGTTACCGTTGATGTTGGCTGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGaaa                                                                                                       <  2:355837/103‑1 (MQ=255)                                                            atttcgtatCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGaaa                                                                                                       >  1:355837/7‑109 (MQ=255)                                                             cGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCgtg                                                                        >  1:842679/1‑139 (MQ=255)                                                                    tCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATcg                                                                 >  2:745376/1‑139 (MQ=255)                                                                     cccGAATGCCGCGTGTTACCGTTGATGTTGGCGGAGTCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTg                                                                                                                               <  2:1118220/76‑1 (MQ=255)                                                                     cccGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTg                                                                                                                               >  1:1118220/1‑76 (MQ=255)                                                                     cccGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATcgc                                                                >  1:1036058/1‑139 (MQ=255)                                                                          aTGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCt                                                           >  1:428591/1‑139 (MQ=255)                                                                              cgcgTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTtgct                                                       >  2:680746/1‑139 (MQ=255)                                                                              cacgTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACct                                                                                                                          <  2:387819/70‑1 (MQ=255)                                                                              cacgTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACct                                                                                                                          >  1:387819/3‑72 (MQ=255)                                                                                cgTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTgaga                                                                                                               <  1:201039/81‑1 (MQ=255)                                                                                cgTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTgaga                                                                                                               >  2:201039/1‑81 (MQ=255)                                                                                   gtTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTAt                                                  <  2:428591/139‑1 (MQ=255)                                                                                       ccGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTga                                                                                                                 >  2:1793104/1‑72 (MQ=255)                                                                                       ccGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTga                                                                                                                 <  1:1793104/72‑1 (MQ=255)                                                                                         gTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGcc                                            <  1:837845/139‑1 (MQ=255)                                                                                         gTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGcc                                            <  2:1036058/139‑1 (MQ=255)                                                                                           tGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTg                                          <  2:600943/139‑1 (MQ=255)                                                                                            gATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATcgcg                                                               <  1:153479/117‑1 (MQ=255)                                                                                            gATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATcgcg                                                               >  2:153479/1‑117 (MQ=255)                                                                                              tGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCAt                                       <  1:745376/139‑1 (MQ=255)                                                                                                    cGGAATCACAGTGATGACGAGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGGGACTCTGAAACAGAGACTGGATGCCACAAATCAATTTCGTGTGGATCGCGGGCGTGCTGCTATGGGGCCTGTATTCCAAc                                 <  1:1494363/139‑1 (MQ=255)                                                                                                    cGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAAc                                 <  1:1775774/139‑1 (MQ=255)                                                                                                        aTCACAGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGt                             <  2:1268625/139‑1 (MQ=255)                                                                                                             aGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTgaga                                                                                                               <  2:1054525/52‑1 (MQ=255)                                                                                                             aGTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTgaga                                                                                                               >  1:1054525/1‑52 (MQ=255)                                                                                                              gTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGTGCCTGCATTCCAACAGGTCTACAg                       >  2:3427/1‑139 (MQ=255)                                                                                                              gTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAg                       >  2:16069/1‑139 (MQ=255)                                                                                                              gTCATGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAg                       >  2:661804/1‑139 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTAc                                                                                                                            <  1:1002669/35‑1 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTAc                                                                                                                            >  2:1002669/1‑35 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAAt                                                                              >  2:1398096/1‑81 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAAt                                                                              <  2:199221/81‑1 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAAt                                                                              <  1:1398096/81‑1 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAAt                                                                              >  1:199221/1‑81 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTg                                                                            <  2:933108/83‑1 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTg                                                                            >  1:933108/1‑83 (MQ=255)                                                                                                                 aTGACGGGTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCt                    >  2:1544889/1‑139 (MQ=255)                                                                                                                       ggTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATcg                                                                 >  1:678520/1‑88 (MQ=255)                                                                                                                       ggTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATcg                                                                 <  2:678520/88‑1 (MQ=255)                                                                                                                        gTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGa                                                                                            <  1:1307565/60‑1 (MQ=255)                                                                                                                        gTAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGa                                                                                            >  2:1307565/1‑60 (MQ=255)                                                                                                                         tAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGcc                                            <  2:1354769/107‑1 (MQ=255)                                                                                                                         tAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGcc                                            >  1:1354769/1‑107 (MQ=255)                                                                                                                         tAGCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGa            <  1:680746/139‑1 (MQ=255)                                                                                                                           gCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACa          >  1:29240/1‑139 (MQ=255)                                                                                                                           gCAAATCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACa          <  1:1544889/139‑1 (MQ=255)                                                                                                                               aTCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACAttgt      >  2:1262608/1‑139 (MQ=255)                                                                                                                                tCAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACAttgtt     <  1:1262608/139‑1 (MQ=255)                                                                                                                                 cAGGCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACAttgttg    <  2:29240/139‑1 (MQ=255)                                                                                                                                   ggCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCgtg                                                                        >  2:1249706/1‑69 (MQ=255)                                                                                                                                   ggCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCa  <  1:16069/139‑1 (MQ=255)                                                                                                                                   ggCAATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCgtg                                                                        <  1:1249706/69‑1 (MQ=255)                                                                                                                                    gCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGa            >  1:1763925/1‑128 (MQ=255)                                                                                                                                    gCGATGCGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGa            <  2:1763925/128‑1 (MQ=255)                                                                                                                                         |                                                                                                                                     TTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCA  >  NC_000913/3991012‑3991280

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