Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,298,294	G→A	T142I (ACC→ATC)	napB ←	nitrate reductase, small, cytochrome C550 subunit, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,298,294	0	G	A	100.0%	119.7 / NA	35	T142I (ACC→ATC)	napB	nitrate reductase, small, cytochrome C550 subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/18);  total (17/18)

CGATCAACAGCAGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGGTAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAGTAGGACTGATACGCGG  >  NC_000913/2298163‑2298427                                                                                                                                    |                                                                                                                                      cGATCAACAGCAGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAAGCTTACGGTCAGAAGTTCCCGTAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGt                                                                                                                                <  1:989902/139‑1 (MQ=255)                tCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTgg                                                                                                                 >  1:1035697/1‑139 (MQ=255)                tCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTgg                                                                                                                 >  2:335813/1‑139 (MQ=255)                                    ggggTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCGACGATTGGCGCGGTATCGGCCTGCGGTa                                                                                             <  2:730943/139‑1 (MQ=255)                                              accaTTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTg                                                                                   <  1:197252/139‑1 (MQ=255)                                                   ttCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACAATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGcagac                                                                              <  2:534645/139‑1 (MQ=255)                                                            cgcTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAac                                                                     >  2:436965/1‑139 (MQ=255)                                                              cTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAacga                                                                   >  2:131735/1‑139 (MQ=255)                                                                     aGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAg                                                            >  1:580359/1‑139 (MQ=255)                                                                           ggCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTg                                                                                                  <  1:224905/95‑1 (MQ=255)                                                                           ggCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTg                                                                                                  >  2:224905/1‑95 (MQ=255)                                                                                 cGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTGTCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCgcgc                                                >  2:943978/1‑139 (MQ=255)                                                                                      aGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGa                                                             >  1:729299/1‑121 (MQ=255)                                                                                      aGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGa                                                             <  2:729299/121‑1 (MQ=255)                                                                                      aGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCcgcg                                                 <  2:376901/133‑1 (MQ=255)                                                                                      aGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCcgcg                                                 >  1:376901/1‑133 (MQ=255)                                                                                                 taataaCCTCTTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCt                                <  1:131735/139‑1 (MQ=255)                                                                                                         tctTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGcagac                                                                              >  1:1116171/1‑85 (MQ=255)                                                                                                         tctTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGcagac                                                                              <  2:1116171/85‑1 (MQ=255)                                                                                                          ctTATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTTGCTGTCCATaaa                       <  1:436965/139‑1 (MQ=255)                                                                                                            tATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGgcg                                                                                                               <  1:1039360/49‑1 (MQ=255)                                                                                                            tATTTCCCGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGgcg                                                                                                               >  2:1039360/1‑49 (MQ=255)                                                                                                                 cccGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACg                                                                                           <  2:155946/64‑1 (MQ=255)                                                                                                                 cccGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACg                                                                                           >  1:155946/1‑64 (MQ=255)                                                                                                                   cGTAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAGTAgg              >  1:866611/1‑139 (MQ=255)                                                                                                                     tAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAg                  <  2:1109911/133‑1 (MQ=255)                                                                                                                     tAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAg                  <  2:48634/133‑1 (MQ=255)                                                                                                                     tAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAg                  >  1:1109911/1‑133 (MQ=255)                                                                                                                     tAACCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAg                  >  1:48634/1‑133 (MQ=255)                                                                                                                       aCCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGAcgc                                                                <  1:257778/85‑1 (MQ=255)                                                                                                                       aCCTTTTGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGAcgc                                                                >  2:257778/1‑85 (MQ=255)                                                                                                                            ttGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGcag                                                                                >  1:808548/1‑64 (MQ=255)                                                                                                                            ttGATGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGcag                                                                                <  2:808548/64‑1 (MQ=255)                                                                                                                               aTGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATGACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAGTAGGACTGATACgcgg  <  1:943978/139‑1 (MQ=255)                                                                                                                               aTGGGATAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAGTAGGACTGATACgcgg  <  2:1035697/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      CGATCAACAGCAGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCTTATTTCCCGTAACCTTTTGATGGGGTAAAGGTATTCCCCACGATTGGCGCGGTATCGGCCTGCGGTACGTGACATTGCAGACAGAAATAACGACGCGGAGCCACTTCCGCGCCCACTTTGCCGTCGCTGTCCATAAAGTGAGTAGGACTGATACGCGG  >  NC_000913/2298163‑2298427

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