Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,286,859	T→C	intergenic (+333/+207)	narI → / ← rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,286,859	0	T	C	100.0%	106.3 / NA	36	intergenic (+333/+207)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/15);  total (21/15)

AAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTC  >  NC_000913/1286810‑1286994                                                  |                                                                                                                                          aaGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                <  1:1397881/139‑1 (MQ=17)  aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                                                                      <  1:106013/116‑1 (MQ=255)  aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                                                                      >  2:106013/1‑116 (MQ=12)   gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCa                                              >  1:1573105/1‑139 (MQ=9)          aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                                       <  1:993174/139‑1 (MQ=17)           tCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                      <  2:1856116/139‑1 (MQ=18)               tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                  >  2:891293/1‑139 (MQ=12)                  ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                 >  2:492535/1‑121 (MQ=2)                  ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                 <  1:492535/121‑1 (MQ=17)                    gcatGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                <  2:321170/117‑1 (MQ=255)                    gcatGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                >  1:321170/4‑120 (MQ=19)                      atGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTggggg                                                                                                 <  2:276418/68‑1 (MQ=255)                      atGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTggggg                                                                                                 >  1:276418/2‑69 (MQ=255)                      atGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                           <  2:1712897/138‑1 (MQ=255)                        gggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                         >  2:966724/1‑139 (MQ=255)                          gCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                       >  1:1874554/1‑139 (MQ=38)                                   ttGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCtt              >  1:634563/1‑139 (MQ=255)                                   ttGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCtt              >  2:899737/1‑139 (MQ=255)                                        tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCg         >  2:128343/1‑139 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGCAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                                       <  1:1344340/107‑1 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAt                                                                                                              <  1:1051176/38‑1 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAt                                                                                                              >  2:1051176/1‑38 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCg         >  1:1353470/1‑139 (MQ=14)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTc                                                           >  2:1708771/1‑87 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTc                                                           <  1:1708771/87‑1 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                                       >  2:1344340/1‑107 (MQ=38)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGaa       >  2:1375077/1‑139 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGaa       >  1:483342/1‑139 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGaa       >  1:453929/1‑139 (MQ=255)                                          cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGaa       >  2:932810/1‑139 (MQ=255)                                           ggATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGa        <  2:1353470/139‑1 (MQ=14)                                            gATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                     <  1:989089/123‑1 (MQ=255)                                            gATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                     >  2:989089/1‑123 (MQ=255)                                            gATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGc     >  1:619309/1‑139 (MQ=255)                                              ttCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCtt   <  1:899737/139‑1 (MQ=255)                                               tCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGAAGCTtc  <  1:1375077/139‑1 (MQ=255)                                                  |                                                                                                                                          AAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATG‑‑GTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTC  >  NC_000913/1286810‑1286994

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