Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,481,269	G→A	P857S (CCT→TCT)	valS ←	valine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,481,269	0	G	A	100.0%	137.4 / NA	38	P857S (CCT→TCT)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/23);  total (15/23)

GAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCG  >  NC_000913/4481192‑4481346                                                                              |                                                                              gaGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGgcagacag                                                                            >  2:927423/1‑81 (MQ=255)   ggCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGacagca                                                                          >  1:927419/1‑81 (MQ=255)    gccagccaTCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCAc                                                                         <  2:927365/81‑1 (MQ=255)          caTCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGa                                                                    >  1:927418/1‑80 (MQ=255)           aTCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATa                                                                  >  1:927427/1‑81 (MQ=255)            tCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATAc                                                                 <  2:927366/81‑1 (MQ=255)                gATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTc                                                             <  1:927367/81‑1 (MQ=255)                    agcagcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAg                                                          <  1:927368/80‑1 (MQ=255)                    agcagcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGa                                                         >  2:927424/1‑81 (MQ=255)                     gcagcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAc                                                        <  1:927369/81‑1 (MQ=255)                       agcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAcgc                                                      <  1:927370/81‑1 (MQ=255)                       agcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAcgc                                                      <  1:927372/81‑1 (MQ=255)                       agcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAcgc                                                      <  2:927371/81‑1 (MQ=255)                        gcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAcgcg                                                     <  2:927373/81‑1 (MQ=255)                            tGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGcc                                                   <  2:927374/79‑1 (MQ=255)                            tGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGcc                                                   <  2:927375/79‑1 (MQ=255)                              cACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAggg                                               <  2:927376/81‑1 (MQ=255)                                    cGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTg                                           >  2:927443/1‑79 (MQ=255)                                     gatgatCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAg                                        <  1:927377/81‑1 (MQ=255)                                      atgatCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAgg                                       <  2:927378/81‑1 (MQ=255)                                       tgatCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGa                                      <  2:927379/81‑1 (MQ=255)                                        gatCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAgg                                       >  2:927433/1‑79 (MQ=255)                                         atCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGaa                                     <  2:927381/80‑1 (MQ=255)                                         atCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGaa                                     <  2:927380/80‑1 (MQ=255)                                         atCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAg                                    >  2:927435/1‑81 (MQ=255)                                              cGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAg                                    <  1:927382/76‑1 (MQ=255)                                               gTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACgg                              >  1:927445/1‑81 (MQ=255)                                               gTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACgg                              >  2:927430/1‑81 (MQ=255)                                                tAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGt                             <  1:927383/81‑1 (MQ=255)                                                tAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGt                             <  1:927384/81‑1 (MQ=255)                                                  aCGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGttt                           >  2:927436/1‑81 (MQ=255)                                                  aCGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGtt                            <  2:927385/80‑1 (MQ=255)                                                   cGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGtttt                          >  2:927449/1‑81 (MQ=255)                                                     gAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCa                        >  2:927439/1‑81 (MQ=255)                                                      aaaCCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCAt                       >  2:927453/1‑81 (MQ=255)                                                       aaCCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCAtt                      <  1:927386/81‑1 (MQ=255)                                                           ggACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTAcg                  >  1:927437/1‑81 (MQ=255)                                                           ggACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTAcg                  >  1:927447/1‑81 (MQ=255)                                                                  tGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGtt           >  1:927440/1‑81 (MQ=255)                                                                    tCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCt         <  2:927387/81‑1 (MQ=255)                                                                          gcaggcagTACGCTGCTATTTTCCAGCCGCGTCAGGCTTTGCCTGAATCCACGCTTTTCATTTCCGCGACGTTCTGCATcc   >  2:1051799/1‑81 (MQ=255)                                                                           cagacagCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCcg  <  2:927388/81‑1 (MQ=255)                                                                              |                                                                              GAGGCCAGCCATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCCG  >  NC_000913/4481192‑4481346

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