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

CATATCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAACCAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCATCTGCCAGCAC  >  NC_000913/672822‑673079                                                                                                                              |                                                                                                                                     cATATCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCt                                                                                                                         <  1:567781/139‑1 (MQ=255)      cAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGc                                                                                                                    <  1:1089923/139‑1 (MQ=255)                  cAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCttt                                                                                                        >  1:1255686/1‑139 (MQ=255)                                 gACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCttt                                                                                         >  1:243415/1‑139 (MQ=255)                                  aCGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTtc                                                                                        >  2:784477/1‑139 (MQ=255)                                   cGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTtct                                                                                       >  2:766523/1‑139 (MQ=255)                                   cGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTtct                                                                                       >  2:368243/1‑139 (MQ=255)                                            tCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTc                                                                                                                                >  2:975296/1‑89 (MQ=255)                                            tCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTc                                                                                                                                <  1:975296/89‑1 (MQ=255)                                            tCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGtt                                                                              >  1:994772/1‑139 (MQ=255)                                            tCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGtt                                                                              <  1:368243/139‑1 (MQ=255)                                             ccGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTTATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGTCCTTTCTCGTCACGTTc                                                                             <  1:784477/139‑1 (MQ=255)                                                    gTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATa                                                                      >  2:335767/1‑139 (MQ=255)                                                      aaCGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGc                                                                    >  1:892107/1‑139 (MQ=255)                                                       aCGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCa                                                                   <  2:994772/139‑1 (MQ=255)                                                                 aTCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCgggg                                                         <  1:766523/139‑1 (MQ=255)                                                                 aTCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCgggg                                                         <  2:243415/139‑1 (MQ=255)                                                                  tCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCg                                                                                                                     >  1:1290493/1‑78 (MQ=255)                                                                  tCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCg                                                                                                                     <  2:1290493/78‑1 (MQ=255)                                                                   cACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGaa                                                       >  2:158471/1‑139 (MQ=255)                                                                            gTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGt                                                >  1:1178021/1‑137 (MQ=255)                                                                            gTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGt                                                <  2:1178021/137‑1 (MQ=255)                                                                                aTACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGtt                                          >  1:1047293/1‑139 (MQ=255)                                                                                      ttgttCTTCGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGt                                    >  1:185294/1‑139 (MQ=255)                                                                                         ttcttcGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCg                                                                                                                     <  1:274936/55‑1 (MQ=255)                                                                                         ttcttcGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCg                                                                                                                     >  2:274936/1‑55 (MQ=255)                                                                                         ttcttcGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTc                                                                                                       >  2:427277/1‑69 (MQ=255)                                                                                         ttcttcGACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTc                                                                                                       <  1:427277/69‑1 (MQ=255)                                                                                               gACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCaa                           <  1:158471/139‑1 (MQ=255)                                                                                               gACTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCaa                           <  1:1057683/139‑1 (MQ=255)                                                                                                 cTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGtt                                               >  2:360350/1‑117 (MQ=255)                                                                                                 cTTGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGtt                                               <  1:360350/117‑1 (MQ=255)                                                                                                   tGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTc                                                                                                       >  2:156300/1‑59 (MQ=255)                                                                                                   tGGACATTTTGCTCATGCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTc                                                                                                       <  1:156300/59‑1 (MQ=255)                                                                                                                   gCCGGTATAAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCATCTGcc       <  2:892107/139‑1 (MQ=255)                                                                                                                        tataAAACAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCATCTGCCAGCAc  >  1:855146/1‑139 (MQ=255)                                                                                                                              |                                                                                                                                     CATATCAGCCGGAGAAGCAAACATCATAAACAGACGAACGGTGTCCGCGCCGTAACGTTCAACCATCACCTGCGGGTCGATACCGTTGTTCTTCGACTTGGACATTTTGCTCATGCCGGTATAAACCAGTTCATGGCCTGCCGCATCTTTCGCTTTCACGATACGGCCTTTCTCGTCACGTTCAACGATAGCATCAACCGGGGAAACCCAGTTACGTTCGCCGTTTTCGCCAACATAGTAGAAGGCATCTGCCAGCAC  >  NC_000913/672822‑673079

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