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%	243.4 / NA	75	A26T (GCC→ACC)	hmp	fused nitric oxide dioxygenase/dihydropteridine reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/39);  total (36/39)

TATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAA  >  NC_000913/2685776‑2686034                                                                                                                                       |                                                                                                                             tATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCa                                                                                                                          <  1:174687/139‑1 (MQ=255)      tttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTAC‑GGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    <  1:2012710/139‑1 (MQ=255)           aaGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                                >  1:1489433/1‑139 (MQ=255)            aGATGCATTTGAGATACATCAATTAAGATGCAACAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAcc                                                                                                               <  1:1248732/139‑1 (MQ=255)            aGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAcc                                                                                                               <  2:1057541/139‑1 (MQ=255)               tGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           >  1:755132/1‑124 (MQ=255)               tGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           <  2:755132/124‑1 (MQ=255)                   tttGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGtt                                                                                                        <  1:1804671/139‑1 (MQ=255)                         aTACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  <  2:412551/139‑1 (MQ=255)                         aTACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  <  2:1489433/139‑1 (MQ=255)                                   tAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGaa                                                                                        >  2:1233783/1‑139 (MQ=255)                                       aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCa                                                                                    >  1:1183372/1‑139 (MQ=255)                                                       aCCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                              >  2:1128689/1‑97 (MQ=255)                                                       aCCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                              <  1:1128689/97‑1 (MQ=255)                                                         cATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGa                                                                  <  2:68406/139‑1 (MQ=255)                                                           tatGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                              >  1:778735/1‑93 (MQ=255)                                                           tatGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                              <  2:778735/93‑1 (MQ=255)                                                           tatGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  <  1:1136400/105‑1 (MQ=255)                                                           tatGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  >  2:1136400/1‑105 (MQ=255)                                                              gCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                             <  2:1183372/139‑1 (MQ=255)                                                                 tGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           <  2:1005953/74‑1 (MQ=255)                                                                 tGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           >  1:1005953/1‑74 (MQ=255)                                                                 tGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAAccc                                                                                            <  1:2038173/105‑1 (MQ=255)                                                                 tGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAAccc                                                                                            >  2:2038173/1‑105 (MQ=255)                                                                  gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                                >  1:1768677/1‑84 (MQ=255)                                                                  gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                                <  2:1768677/84‑1 (MQ=255)                                                                  gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                              <  1:1599851/86‑1 (MQ=255)                                                                  gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                              >  2:1599851/1‑86 (MQ=255)                                                                  gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGc                                                         >  1:1587349/1‑139 (MQ=255)                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    <  1:192460/79‑1 (MQ=255)                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    >  2:192460/1‑79 (MQ=255)                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                              <  2:1918036/85‑1 (MQ=255)                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                              >  1:1918036/1‑85 (MQ=255)                                                                   aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCg                                                        <  2:1118650/139‑1 (MQ=255)                                                                       tCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGACCAAAATTAACCACCCATTTCTACGACCCTATGTCTACTAATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAAt                                                    >  1:1557201/1‑139 (MQ=255)                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    >  2:1806367/1‑66 (MQ=255)                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    <  1:1806367/66‑1 (MQ=255)                                                                                cGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCaa                                           >  2:1364259/1‑139 (MQ=255)                                                                                     cAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           >  2:1083444/1‑54 (MQ=255)                                                                                     cAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                             <  1:259083/116‑1 (MQ=255)                                                                                     cAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                             >  2:259083/1‑116 (MQ=255)                                                                                     cAGTAAAAGCCACCATCCCTCTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           <  1:1083444/54‑1 (MQ=255)                                                                                                accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           >  2:1818375/1‑43 (MQ=255)                                                                                                accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           <  1:873954/43‑1 (MQ=255)                                                                                                accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           >  2:873954/1‑43 (MQ=255)                                                                                                accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                           <  1:1818375/43‑1 (MQ=255)                                                                                                accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTa                           >  2:2009626/1‑139 (MQ=255)                                                                                                 ccaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTaa                          >  2:1942800/1‑139 (MQ=255)                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTa                                                      <  2:218648/110‑1 (MQ=255)                                                                                                  caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTa                                                      >  1:218648/1‑110 (MQ=255)                                                                                                            cTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  >  1:1261450/1‑56 (MQ=255)                                                                                                            cTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  <  2:1261450/56‑1 (MQ=255)                                                                                                            cTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                 <  2:1270432/73‑1 (MQ=255)                                                                                                            cTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                 >  1:1270432/1‑73 (MQ=255)                                                                                                              ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTa                                                      <  1:632687/98‑1 (MQ=255)                                                                                                              ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTa                                                      >  2:632687/1‑98 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    <  1:1602307/35‑1 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    >  2:332003/1‑35 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    <  1:332003/35‑1 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                    >  2:1602307/1‑35 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTa                                                                                                            <  2:372342/43‑1 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTa                                                                                                            >  1:372342/1‑43 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  >  2:673865/1‑53 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  >  2:691519/1‑53 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  <  1:691519/53‑1 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  <  1:673865/53‑1 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGt                                                                <  2:1031626/87‑1 (MQ=255)                                                                                                               gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGt                                                                >  1:1031626/1‑87 (MQ=255)                                                                                                                tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  >  1:1556462/1‑52 (MQ=255)                                                                                                                tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                  <  2:1556462/52‑1 (MQ=255)                                                                                                                 ggAAACGGGGCCAAAGTTAGCCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACAGGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTGTAACGCTATTGCCGCCTAc          <  1:1364259/139‑1 (MQ=255)                                                                                                                 ggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTAc          <  1:2009626/139‑1 (MQ=255)                                                                                                                  gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                             >  1:162726/1‑87 (MQ=255)                                                                                                                  gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                             <  2:162726/87‑1 (MQ=255)                                                                                                                         ggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTaa  >  1:1677532/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                             TATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAA  >  NC_000913/2685776‑2686034

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: ‑