Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,287,743	C→T	intergenic (‑372/‑21)	acs ← / → nrfA	acetyl‑CoA synthetase (AMP‑forming)/cytochrome c552 nitrite reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,287,743	0	C	T	100.0%	175.1 / NA	54	intergenic (‑372/‑21)	acs/nrfA	acetyl‑CoA synthetase (AMP‑forming)/cytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/20);  total (34/20)

GCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGG  >  NC_000913/4287613‑4287878                                                                                                                                   |                                                                                                                                        gCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCgg                                                                                                                                >  2:524218/1‑140 (MQ=255)   aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGAt                                                                                                                              >  1:705616/1‑140 (MQ=255)   aCTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGCCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGAt                                                                                                                              <  2:819483/140‑1 (MQ=255)      tACAATTGATTAAAGACAACATTCACATTGTGTTTTTTTGTTACACATGGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAAAAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGaa                                                                                                                           <  2:816420/140‑1 (MQ=255)      tACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGaa                                                                                                                           <  2:712417/140‑1 (MQ=255)       aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATTTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAg                                                                                                                          <  1:788100/140‑1 (MQ=255)       aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGaa                                                                                                                           >  2:1182657/1‑139 (MQ=255)       aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGaa                                                                                                                           <  1:1182657/139‑1 (MQ=255)           ttGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGcccc                                                                                                                      <  2:705616/140‑1 (MQ=255)             gATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGaa                                                                                                                           >  1:169105/1‑133 (MQ=255)             gATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTa                                                                                                                    >  1:811103/1‑140 (MQ=255)             gATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGaa                                                                                                                           <  2:169105/133‑1 (MQ=255)              aTTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTAt                                                                                                                   >  2:885644/1‑140 (MQ=255)              aTTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTAt                                                                                                                   >  1:1006865/1‑140 (MQ=255)              aTTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTAt                                                                                                                   >  1:1213318/1‑140 (MQ=255)                 aaaGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGAc                                                                                                                <  2:811103/140‑1 (MQ=255)                       aacaTTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGAt                                                                                                          >  1:178660/1‑140 (MQ=255)                                   gtgGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCTGATGAAGCCCCTATGACAAGGATAAAAATAAacgc                                                                                              >  1:925366/1‑140 (MQ=255)                                    tgtttttttGTTACACATAGGGGCGAGAATTGTAATGACAGTGTAGGTGCGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca                                                                                             <  2:855271/134‑1 (MQ=255)                                    ggttttttttttACAAATGGGGGCGAGAAATGTCATGACAGTGTAGGTGGGTTTACTGTGGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca                                                                                             <  2:4606/134‑1 (MQ=255)                                      gTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcacg                                                                                           >  2:357139/1‑140 (MQ=255)                                        tATTTGTTACACATGGGGGCGAGAAATGCAATGACAGTGTGGGGGGGTTTATGGCTGAGAAAAAAAAAAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGcc                                                                                         <  2:178660/140‑1 (MQ=255)                                            tGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTAt                                                                                     >  1:542626/1‑140 (MQ=255)                                                acacATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTc                                                                                 >  2:197680/1‑140 (MQ=255)                                                      aGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTa                                                                           >  2:610047/1‑140 (MQ=255)                                                      aGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTa                                                                           >  1:364419/1‑140 (MQ=255)                                                        gggCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                                         >  2:555274/1‑140 (MQ=255)                                                             aGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTc                                                                    >  2:328739/1‑140 (MQ=255)                                                             aGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTc                                                                    >  2:1154616/1‑140 (MQ=255)                                                              gCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                                                   >  1:709608/1‑140 (MQ=255)                                                                     gagaagggtgtGGGTGGGGTTACTGTGGTAAAAAACAAAAAAACCCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttc                                                            <  2:43563/133‑1 (MQ=255)                                                                       tGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGGAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctt                                                          >  2:1129125/1‑140 (MQ=255)                                                                          cggtgtGGGTGCGTTTCCTGTCGAGAAAAAAAAAAAAACCCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                                       <  2:709608/138‑1 (MQ=255)                                                                          cAGTGAAGGGGCGGTTACTGTCGAAAAAAAAAGAAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                                       <  2:1116800/140‑1 (MQ=255)                                                                          aaGTGGAGTTGCGTTTACTGTCGTAAAAAACAAAAAAACCCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                                       <  2:364419/139‑1 (MQ=255)                                                                             tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCAc                                                    >  2:812581/1‑140 (MQ=255)                                                                              gtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACt                                                   >  2:538783/1‑140 (MQ=255)                                                                               tAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                                       >  2:198119/1‑135 (MQ=255)                                                                               tAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                                       <  1:198119/135‑1 (MQ=255)                                                                                  gTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg                                               >  1:35973/1‑140 (MQ=255)                                                                                    ggggTTACTGTCGAAACAAAAAAAAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGtt                                             <  2:1135899/138‑1 (MQ=255)                                                                                     cGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc                                            >  2:1218505/1‑140 (MQ=255)                                                                                     cGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc                                            >  1:159034/1‑140 (MQ=255)                                                                                     cGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc                                            >  2:597412/1‑140 (MQ=255)                                                                                      ggTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCa                                           >  1:1135899/1‑140 (MQ=255)                                                                                       gTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTcctgtctc                                          >  1:43563/1‑133 (MQ=255)                                                                                       gTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCAc                                          >  1:646434/1‑140 (MQ=255)                                                                                         tACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGc                                        >  1:1235973/1‑140 (MQ=255)                                                                                                        cAATAAAAACCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc                         >  1:942632/1‑140 (MQ=255)                                                                                                         aaaaaaacCCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgct                        <  2:1190474/137‑1 (MQ=255)                                                                                                               acccGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCa                  <  1:322112/138‑1 (MQ=255)                                                                                                               aaCCGCCATTGCAACAATGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCa                  <  1:1154616/140‑1 (MQ=255)                                                                                                                               aTGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTgg  >  1:713363/1‑140 (MQ=255)                                                                                                                               aTGGTGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTgg  >  1:934876/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                        GCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTTACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACCTGTAACTGTGG  >  NC_000913/4287613‑4287878

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