Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,957,020	G→A	intergenic (‑13/+12)	torZ ← / ← torY	trimethylamine N‑oxide reductase 2/cytochrome c quinol dehydrogenase TorY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,957,020	0	G	A	100.0%	151.7 / NA	44	intergenic (‑13/+12)	torZ/torY	trimethylamine N‑oxide reductase 2/cytochrome c quinol dehydrogenase TorY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/23);  total (21/23)

CCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGG  >  NC_000913/1956884‑1957156                                                                                                                                         |                                                                                                                                         ccGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATc                                                                                                                                       >  1:619342/1‑140 (MQ=255)  cgcgGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCa                                                                                                                                      >  1:365643/1‑140 (MQ=255)      gTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAAc                                                                                                                                    >  1:528900/1‑138 (MQ=255)      gTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAAc                                                                                                                                    <  2:528900/138‑1 (MQ=255)          ggATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGAt                                                                                                                              <  2:461330/140‑1 (MQ=255)                   cGCCCTTCTCTTCTGCCGATGCCGGGAGCGGAGCTGGCGTTGTCTCCTCCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTtc                                                                                                                     <  2:258215/140‑1 (MQ=255)                   cGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTtc                                                                                                                     <  1:543658/140‑1 (MQ=255)                   cGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTtc                                                                                                                     <  1:332587/140‑1 (MQ=255)                   cGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTtc                                                                                                                     <  1:401605/140‑1 (MQ=255)                         tctcttctGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTa                                                                                                               <  1:588683/140‑1 (MQ=255)                                tGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCtt                                                                                                        >  2:663263/1‑140 (MQ=255)                                tGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCtt                                                                                                        <  1:9120/140‑1 (MQ=255)                                    cATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGc                                                                                                    >  2:157365/1‑140 (MQ=255)                                     aTGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATTCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCg                                                                                                   <  2:630664/140‑1 (MQ=255)                                     aTGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCg                                                                                                   <  2:504716/140‑1 (MQ=255)                                     aTGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCg                                                                                                   <  1:553822/140‑1 (MQ=255)                                       gCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGtt                                                                                                 >  1:425555/1‑140 (MQ=255)                                                       ccGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACcgc                                                                                 >  2:596879/1‑140 (MQ=255)                                                        cGATGTCACCACCAGTGCTCCATCGACAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACcgcg                                                                                <  1:278339/140‑1 (MQ=255)                                                        cGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACcgcg                                                                                <  1:160349/140‑1 (MQ=255)                                                            gTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTaa                                                                            <  2:433809/140‑1 (MQ=255)                                                                  accaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATgac                                                                      <  1:663263/140‑3 (MQ=255)                                                                         cTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCg                                                               >  2:652694/1‑140 (MQ=255)                                                                                               tttAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTg                                         >  1:19679/1‑140 (MQ=255)                                                                                               tttAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTGAAAATGTACAGATCATTTTCGATCATGCTGGTTCGTg                                         >  1:312643/1‑140 (MQ=255)                                                                                                       aTTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAtt                                 >  2:9144/1‑140 (MQ=255)                                                                                                       aTTCACGTCTTGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAtt                                 >  1:339605/1‑140 (MQ=255)                                                                                                                ttGTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATCTCCAGTTCTTTTCCGCGCATGCTGGGTCGTTCCCCCATTCCTTTGGCa                        >  2:204174/1‑140 (MQ=255)                                                                                                                  gTTAATGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAAt                      >  2:549052/1‑140 (MQ=255)                                                                                                                     aaTGTCATAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGa                   <  2:425555/140‑1 (MQ=255)                                                                                                                           aTAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCa             <  2:525939/140‑1 (MQ=255)                                                                                                                            tAACTTCCTCCCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAATCCTTTGGCAATGGACGGCCAc            >  1:300713/1‑140 (MQ=255)                                                                                                                                 tcctccCTAATCAACGATGATCGCTGTTTCTCGGTAATAACTTTGGCGATGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCAtt       <  2:258163/140‑1 (MQ=255)                                                                                                                                 tcctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCAtt       <  2:19679/140‑1 (MQ=255)                                                                                                                                  cctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg      <  1:652694/140‑1 (MQ=255)                                                                                                                                  cctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg      >  1:435725/1‑140 (MQ=255)                                                                                                                                  cctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg      >  2:118193/1‑140 (MQ=255)                                                                                                                                  cctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg      >  1:334840/1‑140 (MQ=255)                                                                                                                                   ctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGa     <  2:339605/140‑1 (MQ=255)                                                                                                                                   ctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGa     >  1:649582/1‑140 (MQ=255)                                                                                                                                   ctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGa     <  1:118193/140‑1 (MQ=255)                                                                                                                                   ctccCTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAGTCCGTTGGCAATGGACGGCCACGCATTGa     >  1:258163/1‑140 (MQ=255)                                                                                                                                     cccTAATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg   <  2:619342/140‑1 (MQ=255)                                                                                                                                      ccTAATCAACTAGGATCACTGTTTCTCGGTAATATCTTTCGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACgg  >  1:14730/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                         CCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGG  >  NC_000913/1956884‑1957156

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