Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,686,857	T→C	G341G (GGT→GGC)	hmp →	nitric oxide dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,686,857	0	T	C	100.0%	173.6 / NA	48	G341G (GGT→GGC)	hmp	nitric oxide dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/28);  total (20/28)

TGGCGATGTTCACGCCTTTGCCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGTCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCAT  >  NC_000913/2686728‑2686992                                                                                                                                  |                                                                                                                                        tGGCGATGTTCACGCCTTTGCCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGAtt                                                                                                                               <  2:320976/140‑1 (MQ=255)      aTGTTCACGCCTTTACCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGc                                                                                                                          <  2:615271/140‑1 (MQ=255)       tGTTCACGCCTTTGCCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCa                                                                                                                         <  1:572722/140‑1 (MQ=255)       tGTTCACGCCTTTGCCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCa                                                                                                                         >  1:516777/1‑140 (MQ=255)           cACGCCTTTGCCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACt                                                                                                                     >  2:686572/1‑140 (MQ=255)               ccTTTGCCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGaa                                                                                                                 <  1:686572/140‑1 (MQ=255)                  ttGCCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGt                                                                                                              >  1:46033/1‑140 (MQ=255)                        aTGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTc                                                                                                        >  1:468041/1‑140 (MQ=255)                             gTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGa                                                                                                   <  1:624377/140‑1 (MQ=255)                             gTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGa                                                                                                   <  2:110296/140‑1 (MQ=255)                             gTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGa                                                                                                   <  2:46033/140‑1 (MQ=255)                               tAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATc                                                                                                 >  1:728693/1‑140 (MQ=255)                                  ggAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGa                                                                                              <  1:227533/140‑1 (MQ=255)                                     aCTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACaa                                                                                           >  2:539963/1‑140 (MQ=255)                                      cTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAAt                                                                                          <  2:516777/140‑1 (MQ=255)                                       tGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTc                                                                                                        >  2:710072/1‑125 (MQ=255)                                       tGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTc                                                                                                        <  1:710072/125‑1 (MQ=255)                                       tGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATg                                                                                         <  1:39948/140‑1 (MQ=255)                                       tGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATg                                                                                         <  2:275465/140‑1 (MQ=255)                                       tGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATg                                                                                         <  1:666543/140‑1 (MQ=255)                                            cAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGtt                                                                                    <  1:189800/140‑1 (MQ=255)                                                         gcTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGc                                                                       >  1:648864/1‑140 (MQ=255)                                                         acTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGc                                                                       <  2:833124/139‑1 (MQ=255)                                                                          tGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCa                                                      >  1:527524/1‑140 (MQ=255)                                                                                 gTCAGCCGATCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTAcc                                               >  2:143832/1‑140 (MQ=255)                                                                                 gTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTAcc                                               >  2:30572/1‑140 (MQ=255)                                                                                 gTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTAcc                                               >  2:618335/1‑140 (MQ=255)                                                                                 gTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTAcc                                               >  2:423915/1‑140 (MQ=255)                                                                                   cAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACcgc                                             >  1:585470/1‑140 (MQ=255)                                                                                     gCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGa                                           >  1:619359/1‑140 (MQ=255)                                                                                       cGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGaaa                                         <  1:893515/140‑1 (MQ=255)                                                                                       cGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGaaa                                         <  2:761223/140‑1 (MQ=255)                                                                                       cGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGaaa                                         <  2:728693/140‑1 (MQ=255)                                                                                       cGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGaaa                                         <  1:30572/140‑1 (MQ=255)                                                                                         atcGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACa                                       <  1:143832/138‑1 (MQ=255)                                                                                         aGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACa                                       <  2:468041/140‑1 (MQ=255)                                                                                         aGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACa                                       <  1:423915/140‑1 (MQ=255)                                                                                         aGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACa                                       <  2:648864/140‑1 (MQ=255)                                                                                           cGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGt                                     <  2:585470/140‑1 (MQ=255)                                                                                                     cgcgCTAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCt                           <  1:877261/140‑1 (MQ=255)                                                                                                          tAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTa                                   <  1:30335/127‑1 (MQ=255)                                                                                                          tAAAGGTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTa                                   >  2:30335/1‑127 (MQ=255)                                                                                                              ggTCAGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGaa                  <  1:539963/140‑1 (MQ=255)                                                                                                                  aGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAg              >  2:359070/1‑140 (MQ=255)                                                                                                                  aGTTTGATAGCGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAg              >  2:765599/1‑140 (MQ=255)                                                                                                                            cGAAGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAAcattc    <  1:359070/140‑1 (MQ=255)                                                                                                                              aaGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAAcattcat  >  1:382462/1‑140 (MQ=255)                                                                                                                              aaGGCCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAAcattcat  >  2:221421/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                        TGGCGATGTTCACGCCTTTGCCGATGAAGTTAAGGAACTGGGGCAGTCACTGCCGCGCTTTACCGCGCACACCTGGTATCGTCAGCCGAGCGAAGCCGATCGCGCTAAAGGTCAGTTTGATAGCGAAGGTCTGATGGATTTGAGCAAACTGGAAGGTGCGTTCAGCGATCCGACAATGCAGTTCTATCTCTGCGGCCCGGTTGGCTTCATGCAGTTTACCGCGAAACAGTTAGTGGATCTGGGCGTGAAGCAGGAAAACATTCAT  >  NC_000913/2686728‑2686992

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