Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,279,364	A→G	intergenic (+16/‑500)	narK → / → narG	nitrate/nitrite transporter/nitrate reductase 1, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,279,364	0	A	G	100.0%	67.7 / NA	19	intergenic (+16/‑500)	narK/narG	nitrate/nitrite transporter/nitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (3/16);  total (3/16)

GGTAGCTCGCTGGCATTAACGGGTTCGCCAGTCGGCGCAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAA  >  NC_000913/1279232‑1279476                                                                                                                                     |                                                                                                                 ggTAGCTCGCTGGCATTAACGGGTTCGCCAGTCGGCGCAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGtgct                                                                                                <  1:611122/151‑1 (MQ=255)                          cGCCAGTCGGCGCAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGt                                                                        <  1:540639/150‑1 (MQ=255)                                 cGGCGCAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCacttac                                                                 <  2:593912/150‑1 (MQ=255)                                   gcgcAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACttt                                                              <  1:363491/151‑1 (MQ=255)                                   gcgcAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACttt                                                              <  1:695837/151‑1 (MQ=255)                                             ggTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTcc                                                    <  1:391095/151‑1 (MQ=255)                                                             tataTCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACtt                                    >  1:694374/1‑151 (MQ=255)                                                             tataTCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACtt                                    >  2:432162/1‑151 (MQ=255)                                                             tataTCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACtt                                    >  2:477543/1‑151 (MQ=255)                                                              atatCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTg                                   <  1:432183/151‑1 (MQ=255)                                                              atatCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTg                                   <  2:694332/151‑1 (MQ=255)                                                                atCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTgcc                                 <  2:29549/151‑1 (MQ=255)                                                                   gCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTgccgcc                              <  1:724991/151‑1 (MQ=255)                                                                       gCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCtt                            <  1:725407/149‑1 (MQ=255)                                                                         gTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCAGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTg                        <  2:3150/151‑1 (MQ=255)                                                                              gATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAAc                   <  1:470904/151‑1 (MQ=255)                                                                                  aCCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGaa               <  1:532531/151‑1 (MQ=255)                                                                                      gggCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATta            <  1:586654/150‑1 (MQ=255)                                                                                               atGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATGCTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTaa  <  2:59200/151‑1 (MQ=255)                                                                                                                                     |                                                                                                                 GGTAGCTCGCTGGCATTAACGGGTTCGCCAGTCGGCGCAATGAAGGTATTTTTGATTTTCTATATCGCCTGCGTAGTGATTACCTGGGCGGTATATGGTCGGCATTCTAAAAAATAAACCGTTACTCGTCATACTTCGGGTTACATGTGCTGCGGCTGCGTTCATTCACCCCAGTCACTTACTTTAGTAAGCTCCTGGGATTCATTCACTTGCCGCCTTCCTGTAAACCGAATTATATAGAGTAA  >  NC_000913/1279232‑1279476

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