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

CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGA  >  NC_000913/2685775‑2686040                                                                                                                                        |                                                                                                                                   cTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                                 >  2:1297717/1‑139 (MQ=255)   aTTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCAt                                                                                                                               <  2:478695/139‑1 (MQ=255)        ttATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                          <  2:565178/139‑1 (MQ=255)            aaGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                                      <  1:1297717/139‑1 (MQ=255)            aaGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGAc                                                                                                                      <  2:447601/139‑1 (MQ=255)              gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                         >  2:48553/1‑134 (MQ=255)              gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTAc                                                                                                                         <  1:48553/134‑1 (MQ=255)              gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                                    <  1:928430/139‑1 (MQ=255)                tGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTa                                                                                                                  <  2:674036/139‑1 (MQ=255)                  cATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATg                                                                                                                <  1:848965/139‑1 (MQ=255)                    tttGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGtt                                                                                                              <  1:240795/139‑1 (MQ=255)                            aCATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATa                                                                                                      >  1:1206652/1‑139 (MQ=255)                             cATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATaa                                                                                                     >  1:306426/1‑139 (MQ=255)                             cATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATa                                                                                                      <  2:586231/138‑1 (MQ=255)                             cATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATa                                                                                                      >  1:586231/1‑138 (MQ=255)                                  aTTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAg                                                                                                >  2:993067/1‑139 (MQ=255)                                        ttGCTAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATGGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTAGTCATAAccc                                                                                                  >  1:808462/2‑131 (MQ=255)                                        aTGCTAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAAccc                                                                                                  <  2:808462/131‑1 (MQ=255)                                        aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCATCCATTTCTACGACCGt                                                                                                                   >  2:262617/1‑114 (MQ=255)                                        aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCATCCATTTCTACGACCGt                                                                                                                   <  1:262617/114‑1 (MQ=255)                                        aTGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCa                                                                                          >  1:465387/1‑139 (MQ=255)                                           cAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                       >  1:91517/1‑139 (MQ=255)                                             aaaaaaGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCTTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                                    >  2:988727/1‑108 (MQ=255)                                             aaaaaaGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCTTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                                    <  1:988727/108‑1 (MQ=255)                                                     aaGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                             <  2:43438/139‑1 (MQ=255)                                                          cATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGa                                                                        <  2:306426/139‑1 (MQ=255)                                                             atGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCAt                                                                                                       <  2:114822/105‑1 (MQ=255)                                                             atGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCAt                                                                                                       >  1:114822/1‑105 (MQ=255)                                                                  tGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAg                                                                >  2:276942/1‑139 (MQ=255)                                                                   gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                       >  1:874793/1‑115 (MQ=255)                                                                   gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                       <  2:874793/115‑1 (MQ=255)                                                                          aaaCCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                          >  1:468763/1‑73 (MQ=255)                                                                          aaaCCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                          <  2:468763/73‑1 (MQ=255)                                                                          aaaCCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATgg                                                        <  1:276942/139‑1 (MQ=255)                                                                                     aCAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTg                                             <  2:602299/139‑1 (MQ=255)                                                                                        gTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                        >  1:700214/1‑77 (MQ=255)                                                                                        gTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                        <  2:700214/77‑1 (MQ=255)                                                                                            aaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                       >  1:653423/1‑90 (MQ=255)                                                                                            aaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAg                                                                                       <  2:653423/90‑1 (MQ=255)                                                                                                 accaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTa                                 >  2:9882/1‑139 (MQ=255)                                                                                                  ccaTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTaa                                <  2:465387/139‑1 (MQ=255)                                                                                                               ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc                   <  2:91517/139‑1 (MQ=255)                                                                                                               ggtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc                   >  2:814484/1‑139 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                          <  1:1215857/35‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                          >  2:1215857/1‑35 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACt                                                                                                          >  2:448076/1‑51 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACt                                                                                                          <  1:448076/51‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCa                                                                                                 <  2:639432/60‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCa                                                                                                 >  1:639432/1‑60 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                             <  1:1241963/80‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                             >  2:1241963/1‑80 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                                   >  2:7716/1‑90 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAAc                                                                   <  1:7716/90‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt                  >  2:523897/1‑139 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt                  >  2:978868/1‑139 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATAAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt                  >  2:260241/1‑139 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGttt                                                                                                             >  1:904390/1‑47 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGttt                                                                                                             <  2:904390/47‑1 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                        >  1:133117/1‑52 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                        <  2:133117/52‑1 (MQ=255)                                                                                                                   gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACg               >  2:114572/1‑139 (MQ=255)                                                                                                                     aaCGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGcc             >  1:263098/1‑139 (MQ=255)                                                                                                                                aGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGa  <  1:978868/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGA  >  NC_000913/2685775‑2686040

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