Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,685,910	G→A	A26T (GCC→ACC)	hmp →	fused nitric oxide dioxygenase/dihydropteridine reductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,685,910	0	G	A	100.0%	224.7 / NA	69	A26T (GCC→ACC)	hmp	fused nitric oxide dioxygenase/dihydropteridine reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/35);  total (34/35)

CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAA  >  NC_000913/2685775‑2686043                                                                                                                                        |                                                                                                                                      cTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                                    <  2:732065/139‑1 (MQ=255)    ttttttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCAtt                                                                                                                                 >  1:263954/1‑139 (MQ=255)        ttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                             <  2:1165977/139‑1 (MQ=255)            aaGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                                         <  1:85869/139‑1 (MQ=255)              gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                            >  2:124457/1‑134 (MQ=255)              gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                            <  1:124457/134‑1 (MQ=255)              gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                                       <  2:263954/139‑1 (MQ=255)                         gaTACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTc                                                                                                            >  1:1129975/1‑139 (MQ=255)                             cATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATaa                                                                                                        >  1:99645/1‑139 (MQ=255)                               tCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAAcc                                                                                                      >  1:1104816/1‑139 (MQ=255)                                cAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAAccc                                                                                                     <  2:1104816/139‑1 (MQ=255)                                          gCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGt                                                                                                                      <  1:414596/112‑1 (MQ=255)                                          gCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGt                                                                                                                      >  2:414596/1‑112 (MQ=255)                                                     aaGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                                <  2:1125712/139‑1 (MQ=255)                                                              tGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATa                                                                                                         >  2:294886/1‑105 (MQ=255)                                                              tGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATa                                                                                                         <  1:294886/105‑1 (MQ=255)                                                                  tGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAg                                                                   >  2:827889/1‑139 (MQ=255)                                                                   gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGc                                                                  >  1:79039/1‑139 (MQ=255)                                                                    aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                           <  1:762964/97‑1 (MQ=255)                                                                    aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                           >  2:762964/1‑97 (MQ=255)                                                                    aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCg                                                                 >  2:774169/1‑139 (MQ=255)                                                                      gCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGt                                                                                                                  >  2:32739/1‑88 (MQ=255)                                                                      gCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGt                                                                                                                  <  1:32739/88‑1 (MQ=255)                                                                       cTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTaa                                                              >  1:707311/1‑139 (MQ=255)                                                                             ccATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGa                                                        <  2:99645/139‑1 (MQ=255)                                                                                 cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                                         <  1:482909/70‑1 (MQ=255)                                                                                 cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                                         >  2:482909/1‑70 (MQ=255)                                                                                     aCAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTg                                                <  1:827889/139‑1 (MQ=255)                                                                                          aaaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGcc                                           <  2:1129975/139‑1 (MQ=255)                                                                                           aaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGccc                                          >  1:1123582/1‑139 (MQ=255)                                                                                           aaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGccc                                          >  1:652605/1‑139 (MQ=255)                                                                                              gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                                    <  1:1110502/46‑1 (MQ=255)                                                                                              gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                                    >  2:1110502/1‑46 (MQ=255)                                                                                               ccaccaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGtt                                      <  2:692119/139‑1 (MQ=255)                                                                                                     tCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                                      <  1:152978/101‑1 (MQ=255)                                                                                                     tCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                                      >  2:152978/1‑101 (MQ=255)                                                                                                          ttACTGGTGGAAACGGGGCCAAAGTCAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTg                           <  1:1275211/139‑1 (MQ=255)                                                                                                               ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGAt                                                       >  1:1288361/1‑106 (MQ=255)                                                                                                               ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGAt                                                       <  2:1288361/106‑1 (MQ=255)                                                                                                               ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc                      <  2:652605/139‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                             <  1:1199744/35‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                             >  2:1199744/1‑35 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGttt                                                                                                                >  1:819664/1‑48 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGttt                                                                                                                <  2:819664/48‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                                >  1:135006/1‑80 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                                <  2:135006/80‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGt                                                                         <  1:1000456/87‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGt                                                                         >  2:1000456/1‑87 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                                      <  2:1077227/90‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                                      <  1:1231357/90‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                                      >  2:1231357/1‑90 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                                      >  1:1077227/1‑90 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTg                                                <  2:1048986/112‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTg                                                >  1:1048986/1‑112 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGcc                                           <  2:656197/117‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGcc                                           >  1:656197/1‑117 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc                      <  1:971363/138‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc                      >  2:971363/1‑138 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTGAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt                     >  2:821371/1‑139 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTa                    >  2:337545/1‑139 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTa                    >  1:717983/1‑139 (MQ=255)                                                                                                                       cGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGTTATTGCCGCCTACGCCAg              <  2:717983/139‑1 (MQ=255)                                                                                                                         gggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTa            <  2:1123582/139‑1 (MQ=255)                                                                                                                          ggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                                <  2:1192494/70‑1 (MQ=255)                                                                                                                          ggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                                >  1:1192494/1‑70 (MQ=255)                                                                                                                          ggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgc                       >  2:422538/1‑127 (MQ=255)                                                                                                                          ggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgc                       <  1:422538/127‑1 (MQ=255)                                                                                                                           ggccaaGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAat          >  2:1096191/5‑139 (MQ=255)                                                                                                                                   tAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGaaaa  <  1:337545/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGAAAA  >  NC_000913/2685775‑2686043

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