Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,280,561	Δ1 bp	coding (698/3744 nt)	narG →	nitrate reductase A subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,280,558	0	G	.	100.0%	102.4 / NA	26	coding (695/3744 nt)	narG	nitrate reductase A subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (12/14);  total (12/14)

GGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAA  >  NC_000913/1280425‑1280685                                                                                                                                      |                                                                                                                                ggTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTgggcgag                                                                                                                           >  2:105781/1‑136 (MQ=255)                  ggCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGa                                                                                                         <  1:59580/140‑1 (MQ=255)                    cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAACTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAAt                                                                                                       <  1:231882/140‑1 (MQ=255)                       tGTCGATGGTTTCTTACGCATCGGGTGCATGCTATCTGTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATctg                                                                                                    <  1:268303/140‑1 (MQ=255)                                   cTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa                                                                                        >  2:541258/1‑140 (MQ=255)                                   cTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACa                                                                                        >  2:471755/1‑140 (MQ=255)                                      aCGCATCGGGTGCACGCTATCTCTCGCTGATGGGCTGTACTTGCTTAATCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACt                                                                                     <  2:306452/140‑1 (MQ=255)                                          aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGAACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTGCCGGAATCTGCTGGCTGGTACAACTCCAg                                                                                 >  1:181672/1‑140 (MQ=255)                                                               gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GTGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTc                                                            >  2:605740/1‑140 (MQ=255)                                                                 tGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGGCaa                                                          >  2:569026/1‑140 (MQ=255)                                                                 tGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGCAATCTGCTCACTTGTACAACTCCAGGTACATCATCGCCTGGGGGTCaa                                                          >  2:239276/1‑140 (MQ=25)                                                                   aTTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTAGTACAACGCCAACTACATCATAGCCTGGGGGTCAAAc                                                        >  1:215593/1‑140 (MQ=255)                                                                    ttGGCGGTACTTGCTTAATCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACg                                                       <  2:127054/140‑1 (MQ=255)                                                                    ttGGCGGTACTTGCTTAAGCTTCTACGACTGTTACTGCGACTTGCCCCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACg                                                       <  1:372011/140‑1 (MQ=255)                                                                           tACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCa                                                <  2:345092/140‑1 (MQ=255)                                                                                 cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAcgcg                                          <  1:471755/140‑1 (MQ=255)                                                                                  ttAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGt                                         >  1:2859/1‑140 (MQ=255)                                                                                         tCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCgg                                  <  1:541258/140‑1 (MQ=255)                                                                                          cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGATCATACGTGCCGCTGACGCGTACCCCGGa                                 >  1:312173/1‑140 (MQ=21)                                                                                              gACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GTGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCt                             <  1:605740/140‑1 (MQ=255)                                                                                              gACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCt                             <  2:555944/140‑1 (MQ=255)                                                                                                   gTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActt                        <  2:129854/140‑1 (MQ=255)                                                                                                   gTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActt                        >  1:608195/1‑140 (MQ=255)                                                                                                   gTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActt                        <  1:105781/140‑1 (MQ=255)                                                                                                   gTACTGCGACTTGCCACCTGCGTCTTCGCTAACCT‑GGTCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActt                        <  1:224115/140‑1 (MQ=21)                                                                                                       tGCGACTTACCTCCTGTCTCTCGGCAAACC‑CCGGCCAGCTAACTGACGCACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCttt                    <  2:312173/140‑1 (MQ=18)                                                                                                           aCTTGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTg                >  2:123595/1‑140 (MQ=255)                                                                                                              tGCCTCCTGCGTCTCCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAg             >  2:85/1‑140 (MQ=255)                                                                                                                         tctcCGCAAACCT‑GGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACaaa  <  2:608195/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                GGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTTACTGAAGTGCGTTACAAA  >  NC_000913/1280425‑1280685

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