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%	114.7 / NA	37	V613F (GTT→TTT)	leuS	leucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/17);  total (20/17)

TCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAACCAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCATCTGCCAGCAC  >  NC_000913/672826‑673079                                                                                                                          |                                                                                                                                     tCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATg                                                                                                                             >  1:455799/1‑131 (MQ=255) tCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATg                                                                                                                             <  2:455799/131‑1 (MQ=255)   aGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCa                                                                                                                   <  1:498409/139‑1 (MQ=255)         agaAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTc                                                                                                             >  2:813936/1‑139 (MQ=255)         agaAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTc                                                                                                             >  2:20682/1‑139 (MQ=255)                             gACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCttt                                                                                         >  2:423146/1‑139 (MQ=255)                               cGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTtct                                                                                       >  2:36784/1‑139 (MQ=255)                               cGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTtct                                                                                       >  1:436711/1‑139 (MQ=255)                               cGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTACTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTtct                                                                                       >  2:452069/1‑139 (MQ=255)                                      tgtCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACg                                                                                <  1:423146/139‑1 (MQ=255)                                          cgcgCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGt                                                                                    >  2:784600/1‑131 (MQ=255)                                          cgcgCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGt                                                                                    <  1:784600/131‑1 (MQ=255)                                             gcCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACg                                                                         >  2:757640/1‑139 (MQ=255)                                              ccGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGa                                                                        <  1:20682/139‑1 (MQ=255)                                                gTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATa                                                                      >  2:337223/1‑139 (MQ=255)                                                             aTCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTTCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCgggg                                                         <  2:792121/139‑1 (MQ=255)                                                             aTCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCgggg                                                         <  2:436711/139‑1 (MQ=255)                                                             aTCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCgggg                                                         <  1:337223/139‑1 (MQ=255)                                                              tCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGa                                                        >  2:410152/1‑139 (MQ=255)                                                                     cGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTAGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAg                                                 >  2:569877/1‑139 (MQ=255)                                                                           gATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCt                                                                                                          >  1:301533/1‑76 (MQ=255)                                                                           gATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCt                                                                                                          <  2:301533/76‑1 (MQ=255)                                                                                  ttgttCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGt                                    >  2:253609/1‑139 (MQ=255)                                                                                  ttgttCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGt                                    >  1:491490/1‑139 (MQ=255)                                                                                  ttgttCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAACATCAACCGGGGAAACCCAGTTACGTTCGCCGt                                    >  2:791847/1‑139 (MQ=255)                                                                                        ttcGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGc                              <  1:253609/139‑1 (MQ=255)                                                                                        ttcGACTTGGACATTGTGCTCATGCCGGTATAGAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGAGGAAACCCAGTTACGTTCGCCGTTTTCGc                              <  2:552759/139‑1 (MQ=255)                                                                                        ttcGACTTGGACATATTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCGTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGc                              <  1:569877/139‑1 (MQ=255)                                                                                          cGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCa                            <  1:883592/139‑1 (MQ=255)                                                                                                 gACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACAtagt                     >  1:109942/1‑139 (MQ=255)                                                                                                          cTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAACATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCAt            <  1:791847/139‑1 (MQ=255)                                                                                                             aTGCCGGTATAAAACAGATCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCa                                                                   >  2:51622/1‑81 (MQ=255)                                                                                                             aTGCCGGTATAAAACAGATCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCa                                                                   <  1:51622/81‑1 (MQ=255)                                                                                                              tGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCgggg                                                         >  1:881404/1‑90 (MQ=255)                                                                                                              tGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCgggg                                                         <  2:881404/90‑1 (MQ=255)                                                                                                                    tataAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCATCTGCCAGCAc  <  1:410152/139‑1 (MQ=255)                                                                                                                    tataAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCATCTGCCAGCAc  >  2:804962/1‑139 (MQ=255)                                                                                                                          |                                                                                                                                     TCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAACCAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCATCTGCCAGCAC  >  NC_000913/672826‑673079

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