Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,060,660	G→A	G469S (GGT→AGT)	torA →	trimethylamine N‑oxide reductase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,060,660	0	G	A	100.0%	131.3 / NA	38	G469S (GGT→AGT)	torA	trimethylamine N‑oxide reductase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/17);  total (21/17)

TTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGATAACCA  >  NC_000913/1060533‑1060793                                                                                                                                |                                                                                                                                      tttCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTaa                                                                                                                           >  2:292178/1‑140 (MQ=255) tttCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTaa                                                                                                                           >  1:503212/1‑140 (MQ=255)        ggCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTgcc                                                                                                                    >  1:328442/1‑140 (MQ=255)              aCGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCt                                                                                                              >  2:534830/1‑140 (MQ=255)                  cGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGaaa                                                                                                          <  1:534830/140‑1 (MQ=255)                          ccTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTAt                                                                                                  >  2:8150/1‑140 (MQ=255)                           cTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTAtt                                                                                                 <  1:321967/140‑1 (MQ=255)                                cACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGc                                                                                            <  1:293619/140‑1 (MQ=255)                                cACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGc                                                                                            <  2:503212/140‑1 (MQ=255)                                cACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGc                                                                                            <  2:370386/140‑1 (MQ=255)                                        caGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTa                                                                                    >  1:307593/1‑140 (MQ=255)                                        caGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTa                                                                                    >  1:24155/1‑140 (MQ=255)                                        caGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTa                                                                                    >  1:79982/1‑140 (MQ=255)                                                     ggCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCg                                                                       >  1:385945/1‑140 (MQ=255)                                                      gCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGc                                                                      >  2:481386/1‑140 (MQ=255)                                                         aCAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCa                                                                          <  2:516406/133‑1 (MQ=255)                                                         aCAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCa                                                                          >  1:516406/1‑133 (MQ=255)                                                         aCAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCAt                                                                   >  2:760411/1‑140 (MQ=255)                                                           agcagcACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATca                                                                 >  2:760301/1‑140 (MQ=255)                                                                caCTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGa                                                            >  1:194933/1‑140 (MQ=255)                                                                          aTTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCat                                                  <  2:24155/140‑1 (MQ=255)                                                                          aTTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCat                                                  >  1:402772/1‑140 (MQ=255)                                                                          aTTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCat                                                  >  2:469379/1‑140 (MQ=255)                                                                             gCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCattat                                               >  2:531016/1‑140 (MQ=255)                                                                             gCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCattat                                               >  1:503174/1‑140 (MQ=255)                                                                                         gATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTgcg                                   <  2:426846/140‑1 (MQ=255)                                                                                         gATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTgcg                                   <  2:67782/140‑1 (MQ=255)                                                                                          aTGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTgcgc                                  <  1:760301/140‑1 (MQ=255)                                                                                           tGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCa                                 <  2:194933/140‑1 (MQ=255)                                                                                            gCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCaa                                <  2:385945/140‑1 (MQ=255)                                                                                                    cGAACCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGaaa                        >  1:260072/1‑140 (MQ=255)                                                                                                       aCCGGGGAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACgg                     >  1:171928/1‑140 (MQ=255)                                                                                                          ggggAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTa                  <  2:171928/140‑1 (MQ=255)                                                                                                            ggAAAGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATc                >  2:187350/1‑140 (MQ=255)                                                                                                              aaaGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGc              <  2:503174/140‑1 (MQ=255)                                                                                                              aaaGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGc              <  1:8150/140‑1 (MQ=255)                                                                                                              aaaGTGATCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGc              <  2:79982/140‑1 (MQ=255)                                                                                                                     tCAACTGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGAt       <  2:328442/140‑1 (MQ=255)                                                                                                                          tGGAACAGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGATAACCa  <  1:469379/140‑1 (MQ=255)                                                                                                                                |                                                                                                                                      TTTCTCCGGCTCTACGTCGATTCCGCCTGTTCACGACAACAGTGACTACAAAGGCTACAGCAGCACTATTCCGATTGCCCGTTTTATCGATGCGATCCTCGAACCGGGGAAAGTGATCAACTGGAACGGTAAATCGGTAAAACTGCCGCCGCTGAAAATGTGTATTTTTGCCGGAACTAACCCATTCCATCGCCATCAGCAGATCAACCGCATTATTGAAGGCTTGCGCAAGCTGGAAACGGTTATCGCCATAGATAACCA  >  NC_000913/1060533‑1060793

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