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%	176.9 / NA	59	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 (32/27);  total (32/27)

CCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGA  >  NC_000913/1286801‑1286988                                                           |                                                                                                                                  ccttaatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCtttt                                                                                                               <  1:399878/73‑1 (MQ=11) ccttaatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCtttt                                                                                                               >  2:399878/7‑79 (MQ=11)     aatCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                               >  1:1209825/3‑139 (MQ=12)          aaGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                          <  1:1573270/139‑1 (MQ=12)             cTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                       <  2:1209825/139‑1 (MQ=18)               tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                      <  1:292879/90‑1 (MQ=255)               tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                      >  2:292879/1‑90 (MQ=255)                       ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                             >  2:1800207/1‑139 (MQ=12)                        tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                                                                 <  1:773170/86‑1 (MQ=255)                        tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                                                                 >  2:773170/1‑86 (MQ=255)                        tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACg                                     >  2:1785049/1‑130 (MQ=16)                        tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACg                                     <  1:1785049/130‑1 (MQ=16)                         ccccGCATGGGCAGCATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                      >  1:181351/1‑80 (MQ=11)                         ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtggtg                                                                                               >  1:1040347/1‑71 (MQ=255)                         ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAtggtggtg                                                                                               <  2:1040347/71‑1 (MQ=255)                         ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                      <  2:181351/80‑1 (MQ=255)                         ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt                           >  2:1135122/1‑139 (MQ=12)                          cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGa                                                                                                             <  2:1473899/56‑1 (MQ=11)                          cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGa                                                                                                             >  1:1473899/1‑56 (MQ=14)                          cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACg                                     >  2:1478983/1‑128 (MQ=9)                          cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACg                                     <  1:1478983/128‑1 (MQ=18)                           ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTATGAAAGTGATGGTGGTGGGGGAAGGATCACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                         >  1:1246118/1‑139 (MQ=255)                             gcatGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt                                                  <  2:83552/109‑1 (MQ=16)                             gcatGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt                                                  >  1:83552/4‑112 (MQ=255)                               atGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGaa                                                                                                            <  1:1829205/51‑1 (MQ=14)                               atGGGCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGaa                                                                                                            >  2:1829205/2‑52 (MQ=14)                                 gggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGaa                                                                                                            >  1:1270684/1‑50 (MQ=255)                                 gggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGaa                                                                                                            <  2:1270684/50‑1 (MQ=255)                                 gggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                           <  1:1418186/67‑1 (MQ=255)                                 gggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                           >  2:1418186/1‑67 (MQ=255)                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                                                                                                                     >  1:1203517/1‑40 (MQ=14)                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCggg                                                                                                                     <  2:1203517/40‑1 (MQ=11)                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGa                                                                                                             <  1:1703934/48‑1 (MQ=255)                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGa                                                                                                             >  2:1703934/1‑48 (MQ=255)                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                           >  1:1532732/1‑66 (MQ=255)                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                           <  2:1532732/66‑1 (MQ=255)                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                  >  2:1198541/1‑139 (MQ=38)                                  ggCAGAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                  >  2:904666/1‑139 (MQ=255)                                      gAATATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg              <  1:1198541/139‑1 (MQ=255)                                       aaTATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                           >  2:207209/1‑61 (MQ=255)                                       aaTATTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                           <  1:207209/61‑1 (MQ=255)                                         tatTTGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc           >  2:1129380/1‑139 (MQ=255)                                            ttGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                           >  1:853700/1‑56 (MQ=255)                                            ttGATTGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggggg                                                                                           <  2:853700/56‑1 (MQ=255)                                                 tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt                                                                  <  2:920825/76‑1 (MQ=255)                                                 tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt                                                                  >  1:920825/1‑76 (MQ=255)                                                 tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCg   >  1:1457620/1‑139 (MQ=255)                                                 tGCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAAGGTGGTGGGGGAAGTATTAGTCAGCGCTGCGCGCTTCGCCCTTCTGGTCGTTGCCTGCGGCAACGCTCTCTCGCTAACGCTCGAGTCGAACCTTAATCg   >  1:1356049/1‑139 (MQ=11)                                                  gCGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGa  <  2:1457620/139‑1 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAgtggtg                                                                                                     >  2:1967720/1‑39 (MQ=17)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGAtg                                                                                                     <  1:1967720/39‑1 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTg                                                                         <  2:2039379/67‑1 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTg                                                                         >  1:2039379/1‑67 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                                                                <  1:326744/76‑1 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                                                                >  2:326744/1‑76 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                               >  2:787667/1‑93 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                               <  1:787667/93‑1 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGa  >  1:1747440/1‑138 (MQ=255)                                                   cGGATTCGCTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAATCGa  <  2:1747440/138‑1 (MQ=255)                                                           |                                                                                                                                  CCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGA  >  NC_000913/1286801‑1286988

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: ‑