Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,494,514	T→C	G168G (GGT→GGC)	nirB →	nitrite reductase catalytic subunit NirB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,494,514	0	T	C	100.0%	103.8 / NA	31	G168G (GGT→GGC)	nirB	nitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/16);  total (15/16)

TGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGCGTGCGCGTTCACACCAGCAAAAACACCCTT  >  NC_000913/3494386‑3494643                                                                                                                                 |                                                                                                                                  tGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTg                                                                                                                               >  1:485467/1‑133 (MQ=255) tGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTGAGGTCTGGAAGCAGCACGCGCGCTGGAAAAATTAGGCATTg                                                                                                                               >  2:213453/1‑133 (MQ=255)        tCTACCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCa                                                                                                                        >  2:60152/1‑133 (MQ=255)          tATCGCACTCTTGAACACCTCAATACTATTGAATCCTGCGCTCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACg                                                                                                                      <  2:250117/133‑1 (MQ=255)             cGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTa                                                                                                                   >  1:301005/1‑133 (MQ=255)                  tATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGaa                                                                                                              >  1:363039/1‑133 (MQ=255)                    ttGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAAtt                                                                                                            <  2:363039/133‑1 (MQ=255)                            cTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTa                                                                                                    >  2:292851/1‑133 (MQ=255)                              cAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCt                                                                                                     >  2:53742/1‑130 (MQ=255)                              cAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCt                                                                                                     <  1:53742/130‑1 (MQ=255)                              cAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATg                                                                                                  >  1:23320/1‑133 (MQ=255)                              cAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATg                                                                                                  >  1:283360/1‑133 (MQ=255)                                aCGCCATTTAATCCTGCGACCGTCGCAGCACACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCt                                                                                                <  1:248403/133‑1 (MQ=255)                                        gAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAg                                                                                        <  2:485467/133‑1 (MQ=255)                                         aaTCCTGCTCCCGTCCCAACAAACGCGGTGACGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGa                                                                                       <  1:199007/133‑1 (MQ=255)                                          aTCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGaa                                                                                      <  2:301005/133‑1 (MQ=255)                                             cTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAg                                                                                   >  1:59206/1‑133 (MQ=255)                                              tGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGc                                                                                  >  1:399007/1‑133 (MQ=255)                                              tGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGc                                                                                  >  1:112043/1‑133 (MQ=255)                                                 gcCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTg                                                                               <  2:112043/133‑1 (MQ=255)                                                          agcaAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATgg                                                                      <  2:498432/133‑1 (MQ=255)                                                                cgcgGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGgcggcg                                                                >  2:290903/1‑133 (MQ=255)                                                                                      gcCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGa                                          <  1:292851/133‑1 (MQ=255)                                                                                                tCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGc                                >  1:298256/1‑133 (MQ=255)                                                                                                  tGGAAGCCGCAGGCGCGCTGAAAAACTTAGGCATGGAAACCCACGTTATCTAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGGTGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGCGt                              <  2:307453/133‑1 (MQ=255)                                                                                                         cGCAGGCGCGCTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGCGTGCGCGtt                       >  1:441242/1‑133 (MQ=255)                                                                                                               cgcgcTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGCGTGCGCGTTCACAcc                 <  2:409630/133‑1 (MQ=255)                                                                                                                 cgcTGAAAAACTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGCGTGCGCGTTCACACCAg               <  1:437420/133‑1 (MQ=255)                                                                                                                       aaaaCTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGCGTGCGCGTTCACACCAGCaaaaa         <  1:60152/133‑1 (MQ=255)                                                                                                                       aaaaCTTAGGCATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGCGTGCGCGTTCACACCAGCaaaaa         <  2:59206/133‑1 (MQ=255)                                                                                                                              aGGCATTGAAACCCAAGTGATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGTGTCGCAAAATCGAAAGTATGGGCGTGCGCGTTCACACCAGCAAAAACACCCtt  <  1:213453/133‑1 (MQ=255)                                                                                                                                 |                                                                                                                                  TGCTTTGTCTATCGCACTATTGAAGACCTCAACGCCATTGAATCCTGCGCCCGTCGCAGCAAACGCGGTGCCGTTGTTGGTGGCGGCCTGTTAGGTCTGGAAGCCGCAGGCGCGCTGAAAAACTTAGGTATTGAAACCCACGTTATCGAATTTGCCCCTATGCTGATGGCAGAACAGCTTGATCAGATGGGCGGCGAGCAGCTGCGTCGCAAAATCGAAAGTATGGGCGTGCGCGTTCACACCAGCAAAAACACCCTT  >  NC_000913/3494386‑3494643

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