Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,958,496	G→A	intergenic (‑364/+24)	torY ← / ← cutC	TMAO reductase III (TorYZ), cytochrome c‑type subunit/copper homeostasis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,958,496	0	G	A	100.0%	158.5 / NA	49	intergenic (‑364/+24)	torY/cutC	TMAO reductase III (TorYZ), cytochrome c‑type subunit/copper homeostasis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/21);  total (28/21)

TTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGA  >  NC_000913/1958363‑1958624                                                                                                                                      |                                                                                                                                 ttttAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGAc                                                                                                                             >  2:665386/1‑139 (MQ=255)                  gcgcACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTaaaa                                                                                                            <  1:665386/139‑1 (MQ=255)                      aCTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCa                                                                                                        >  2:110483/1‑139 (MQ=255)                        tGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCAtt                                                                                                      >  1:1729967/1‑139 (MQ=255)                                gCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGa                                                                                              >  2:9761/1‑139 (MQ=255)                                        gATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTc                                                                                      >  1:1876776/1‑139 (MQ=255)                                        gATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTc                                                                                      >  1:1821552/1‑139 (MQ=255)                                                          ttttttATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCtttcatttc                                                                    >  2:401402/1‑139 (MQ=255)                                                             tttATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGc                                                                 <  2:1729967/139‑1 (MQ=255)                                                                   cTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAAccgc                                                           >  2:1625606/1‑139 (MQ=255)                                                                     gTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAAccgccg                                                         >  2:1735638/1‑139 (MQ=255)                                                                      ttCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCGTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACcgccgc                                                        >  1:1579679/1‑139 (MQ=255)                                                                      ttCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACcgccgc                                                        >  1:1294077/1‑139 (MQ=255)                                                                         aGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGcccc                                                     <  2:1550892/139‑1 (MQ=255)                                                                         aGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGcccc                                                     <  1:9761/139‑1 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACa                                                                                                                            >  1:439357/1‑61 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACa                                                                                                                            <  2:439357/61‑1 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAAc                                                                                                                   >  2:1581064/1‑70 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAAc                                                                                                                   >  1:1833807/1‑70 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAAc                                                                                                                   <  1:1581064/70‑1 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAAc                                                                                                                   <  2:1833807/70‑1 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCaa                                                                                    <  1:78978/101‑1 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCaa                                                                                    >  2:78978/1‑101 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCCATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACg                                              >  1:1673268/1‑139 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCgccgcc                                                       >  1:853/1‑130 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCgccgcc                                                       <  2:853/130‑1 (MQ=255)                                                                                aaaaaaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACg                                              >  2:1708183/1‑139 (MQ=255)                                                                                    aaGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGatat                                          <  1:401402/139‑1 (MQ=255)                                                                                     aGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGAtata                                         >  2:1664405/1‑139 (MQ=255)                                                                                            tGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCgccgcc                                                       <  2:382119/118‑1 (MQ=255)                                                                                            tGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCgccgcc                                                       >  1:382119/1‑118 (MQ=255)                                                                                                      gCAGGGGCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCcg                        >  2:1757348/1‑139 (MQ=255)                                                                                                            gCCTGGTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGtt                  <  1:1735638/139‑1 (MQ=255)                                                                                                                 gTACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCt             >  2:1797319/1‑139 (MQ=255)                                                                                                                  tACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCtg            <  1:1757348/139‑1 (MQ=255)                                                                                                                  tACGAGCAAGCATCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCtg            <  1:110483/139‑1 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAGTCgccgcc                                                       >  2:1428158/1‑92 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAGTCgccgcc                                                       <  1:1428158/92‑1 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGcccc                                                     <  2:1654972/94‑1 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGcccc                                                     >  1:1654972/1‑94 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGt                                                   <  2:1376267/96‑1 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGt                                                   >  1:1376267/1‑96 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGAtata                                         >  2:575426/1‑106 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGAtata                                         <  1:575426/106‑1 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCtgatga        >  1:782509/1‑139 (MQ=255)                                                                                                                      agcaagcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCtgatga        >  1:477297/1‑139 (MQ=255)                                                                                                                          agcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATg    <  2:1294077/139‑1 (MQ=255)                                                                                                                           gcaTCATATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATgg   <  2:477297/139‑1 (MQ=255)                                                                                                                            catcatATTGAGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGa  <  2:1673268/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 TTTTAGCACAACGCTATGCGCACTGTCACCCGCTCAACAGATCACGGTTATCGTTCGTTTTTTATACTGTTCAGGGATAAAAAAAGGCCCCTGTTGAAATTGCAGGGGCCTGGTACGAGCAAGCATCATATTGGGCGACATGATGCAACGGTAAAAATCATTTGGCCTGATGGCGTTCAATGATTCCTTTCATTTCAGCAACCGCCGCCCCGTCTACGATATAACGCGAATACTCGTCCGCGTGTTCATCTGATGACATGGA  >  NC_000913/1958363‑1958624

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