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

CTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCGCCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGA  >  NC_000913/2685775‑2686040                                                                                                                                        |                                                                                                                                   cTATTTTTTATAAGATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                                 >  1:233661/1‑139 (MQ=255)              gATGCATTTGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCg                                                                                                                    <  2:233661/139‑1 (MQ=255)                    tttGAGATACATCAATTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGtt                                                                                                              <  1:1346747/139‑1 (MQ=255)                                  aTTAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAg                                                                                                >  1:1038140/1‑139 (MQ=255)                                   ttAAGATGCAAAAAAAGGAAGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGa                                                                                               <  2:989278/139‑1 (MQ=255)                                                     aaGACCATATGCTTGACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                             <  2:692251/139‑1 (MQ=255)                                                                   gACGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGc                                                               >  2:1083713/1‑139 (MQ=255)                                                                    aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCg                                                              <  2:1038140/139‑1 (MQ=255)                                                                    aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCg                                                              >  1:503012/1‑139 (MQ=255)                                                                    aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGATAACCACCCATTTCTa                                                                                                                          >  1:893421/1‑79 (MQ=255)                                                                    aCGCTCAAACCATCGCTACAGTAAAAGCCACCATCCATTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTa                                                                                                                          <  2:893421/79‑1 (MQ=255)                                                                          aaaCCATCGCTACAGTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATgg                                                        >  1:1073415/1‑139 (MQ=255)                                                                                        gTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                        >  2:639319/1‑77 (MQ=255)                                                                                        gTAAAAGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCa                                                                                                        <  1:639319/77‑1 (MQ=255)                                                                                          aaaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGcc                                        <  1:1083713/139‑1 (MQ=255)                                                                                          aaaaGCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGcc                                        <  2:1073415/139‑1 (MQ=255)                                                                                              gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                                 >  2:851268/1‑46 (MQ=255)                                                                                              gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCAccc                                                                                                                                 <  1:851268/46‑1 (MQ=255)                                                                                              gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACt                                                                                                          <  1:568936/69‑1 (MQ=255)                                                                                              gCCACCATCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACt                                                                                                          >  2:568936/1‑69 (MQ=255)                                                                                                   caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTg                                             >  2:1150019/1‑125 (MQ=255)                                                                                                   caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTg                                             <  1:1150019/125‑1 (MQ=255)                                                                                                   caTCCCTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAAc                               >  2:352468/1‑139 (MQ=255)                                                                                                       ccTTTACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTa                           <  1:352468/139‑1 (MQ=255)                                                                                                          ttACTGGTGGAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTg                        <  1:810759/139‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAg                                                                                                <  2:257274/61‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAg                                                                                                >  1:257274/1‑61 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                             <  1:300344/80‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTaa                                                                             >  2:300344/1‑80 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCg                                                              >  1:648510/1‑95 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCg                                                              <  2:648510/95‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGa                                                     <  2:1215420/104‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGa                                                     >  1:1215420/1‑104 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGGTATTgccgcc                   >  2:891699/1‑138 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTgccgcc                   <  1:891699/138‑1 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt                  >  2:247804/1‑139 (MQ=255)                                                                                                                gtggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCt                  >  2:1109800/1‑139 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATc                                                   >  1:865135/1‑105 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATc                                                   <  2:865135/105‑1 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTa                 >  2:1025779/1‑139 (MQ=255)                                                                                                                 tggAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTa                 >  1:363956/1‑139 (MQ=255)                                                                                                                   gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGt                                                                                                                   <  2:22860/39‑1 (MQ=255)                                                                                                                   gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGt                                                                                                                   >  1:22860/1‑39 (MQ=255)                                                                                                                   gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTa                                                                     >  1:637280/1‑85 (MQ=255)                                                                                                                   gAAACGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTa                                                                     <  2:637280/85‑1 (MQ=255)                                                                                                                    aaaCGGGGCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGc              <  2:363956/139‑1 (MQ=255)                                                                                                                          ggCCAAAGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTaa        >  2:668017/1‑139 (MQ=255)                                                                                                                                aGTTAACCACCCATTTCTACGACCGTATGTTTACTCATAACCCAGAACTCAAAGAAATTTTTAACATGAGTAACCAGCGTAATGGCGATCAACGTGAAGCCCTGTTTAACGCTATTGCCGCCTACGCCAGTAATATTGa  <  1:1109800/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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑