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%	116.8 / NA	33	P857S (CCT→TCT)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/17);  total (16/17)

ATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCC  >  NC_000913/4481202‑4481345                                                                    |                                                                             aTCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATa                                                                 >  2:857539/1‑81 (MQ=255)     ggATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACttt                                                             <  1:857494/81‑1 (MQ=255)      gATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTc                                                            <  1:857495/81‑1 (MQ=255)       aTCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTcc                                                           <  2:857495/81‑1 (MQ=255)            cagcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAcg                                                      <  2:857496/81‑1 (MQ=255)            cagcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAcg                                                      <  1:857498/81‑1 (MQ=255)             agcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAcgc                                                     >  1:857571/1‑81 (MQ=255)             agcTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGAcgc                                                     <  1:857499/81‑1 (MQ=255)                tctGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACgcgc                                                   >  1:857539/1‑80 (MQ=255)                      ccGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGtt                                            >  1:857561/1‑81 (MQ=255)                       cGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGttt                                           <  2:857500/81‑1 (MQ=255)                       cGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGttt                                           <  2:857499/81‑1 (MQ=255)                           gatgatCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCa                                        <  2:857501/80‑1 (MQ=255)                           gatgatCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAg                                       >  2:857547/1‑81 (MQ=255)                             tgatCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGa                                     <  2:857502/81‑1 (MQ=255)                                   tCGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCAc                               >  2:857542/1‑81 (MQ=255)                                    cGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAGGcc                                 <  2:857504/78‑1 (MQ=255)                                    cGTAACGGAAACCGGACCTTTGTCATCGGCAGACAGCACGATGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAg                                   <  2:857503/76‑1 (MQ=255)                                       aaCGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGtt                           >  1:857545/1‑81 (MQ=255)                                        aCGGAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGttt                          <  2:857505/81‑1 (MQ=255)                                          ggAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGtttt                         >  1:857547/1‑80 (MQ=255)                                          ggAAACCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGtttt                         >  1:857546/1‑80 (MQ=255)                                            aaaCCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCAt                      >  2:857506/1‑81 (MQ=255)                                              aCCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCAttt                    <  2:857507/81‑1 (MQ=255)                                              aCCGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCAtt                     >  2:857557/1‑80 (MQ=255)                                                cGGACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTAc                  >  2:857578/1‑81 (MQ=255)                                                  gACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTTATTTAcgc                <  1:857509/81‑1 (MQ=255)                                                  gACCTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTAcgc                >  1:857560/1‑81 (MQ=255)                                                     cTTTGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGAc             <  1:857510/81‑1 (MQ=255)                                                        tGTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGt           <  2:857510/80‑1 (MQ=255)                                                         gTCATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTc         >  1:857552/1‑81 (MQ=255)                                                           cATCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTg       >  2:857566/1‑81 (MQ=255)                                                            aTCGGCAGACAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTg       >  1:857563/1‑80 (MQ=255)                                                               ggcagacaGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATc   >  1:857565/1‑81 (MQ=255)                                                               ggcagacaGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATc   <  2:857511/81‑1 (MQ=255)                                                                gcagacagCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATcc  <  1:857513/81‑1 (MQ=255)                                                                    |                                                                             ATCGGGATCAGCAGCTCTGCACCGTCGATGATCTTCGTAACGGAAACCGGACCTTTGTCATCGGCAGGCAGCACGGTGATACTTTCCAGACGCGCCAGGGTTTGCAGGAAGCCACGGTTTTCATTTACGCGACGTTCTGCATCC  >  NC_000913/4481202‑4481345

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