Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	672,947	C→A	V613F (GTT→TTT)	leuS ←	leucyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	672,947	0	C	A	100.0%	109.2 / NA	35	V613F (GTT→TTT)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/19);  total (16/19)

TATCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAACCAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGG  >  NC_000913/672824‑673066                                                                                                                            |                                                                                                                        tatCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGc                                                                                                          <  2:361717/139‑1 (MQ=255)   tCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCg                                                                                                        <  2:280150/139‑1 (MQ=255)         ggAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCtt                                                                                                  >  1:436872/1‑139 (MQ=255)           agaAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCtt                                                                                                  <  1:363801/137‑1 (MQ=255)           agaAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCtt                                                                                                  >  2:363801/1‑137 (MQ=255)             aaGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGc                                                                                              <  1:478051/139‑1 (MQ=255)                       catAAACAGACGAACGGTGTCCGCGCCGTAACGTTCATCCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATa                                                                                    <  1:681015/139‑1 (MQ=255)                          aaaCAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACgg                                                                                 <  2:436872/139‑1 (MQ=255)                               gACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCttt                                                                            >  2:582192/1‑139 (MQ=255)                                 cGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTtct                                                                          >  2:205499/1‑139 (MQ=255)                                 cGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTtct                                                                          >  2:638764/1‑139 (MQ=255)                                      ggTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCa                                                                     <  2:106042/139‑1 (MQ=255)                                                ccGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGa                                                           <  1:247974/139‑1 (MQ=255)                                                 cGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGAt                                                          <  1:638764/139‑1 (MQ=255)                                                    aaCGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGc                                                       <  2:600887/139‑1 (MQ=255)                                                         tCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGc                                                                                              <  1:656539/95‑1 (MQ=255)                                                         tCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGc                                                                                              >  2:656539/1‑95 (MQ=255)                                                           aaCCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCa                                                                     >  2:795996/1‑118 (MQ=255)                                                           aaCCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCa                                                                     <  1:795996/118‑1 (MQ=255)                                                                  aCCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGaaa                                         >  2:43073/1‑139 (MQ=255)                                                                          gTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGa                                                           <  2:397316/113‑1 (MQ=255)                                                                          gTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGa                                                           >  1:397316/1‑113 (MQ=255)                                                                                    ttgttCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCttt                                                                                           >  2:573905/1‑71 (MQ=255)                                                                                    ttgttCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCttt                                                                                           <  1:573905/71‑1 (MQ=255)                                                                                    ttgttCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGt                       >  2:700421/1‑139 (MQ=255)                                                                                     tgttCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGtt                      <  1:498952/139‑1 (MQ=255)                                                                                       ttcttcGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGtttt                    >  2:648024/1‑139 (MQ=255)                                                                                            cGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGc                                                       <  1:697934/99‑1 (MQ=255)                                                                                            cGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGc                                                       >  2:697934/1‑99 (MQ=255)                                                                                                      aTTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGa     <  1:43073/139‑1 (MQ=255)                                                                                                         ttGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAgg  >  2:192431/1‑139 (MQ=255)                                                                                                                tGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCa                                                               <  2:307642/71‑1 (MQ=255)                                                                                                                tGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCa                                                               >  1:307642/1‑71 (MQ=255)                                                                                                                  ccGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGtttt                    <  2:580790/112‑1 (MQ=255)                                                                                                                  ccGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGtttt                    >  1:580790/1‑112 (MQ=255)                                                                                                                            |                                                                                                                        TATCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAACCAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGG  >  NC_000913/672824‑673066

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