Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,821,522	G→A	P164S (CCA→TCA)	alaS ←	alanine‑‑tRNA ligase/DNA‑binding transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,821,522	0	G	A	100.0%	171.8 / NA	48	P164S (CCA→TCA)	alaS	alanine‑‑tRNA ligase/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/30);  total (18/30)

ATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGATTTCGTAGGCT  >  NC_000913/2821448‑2821598                                                                           |                                                                             aTCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCAcc                                                                        >  2:508924/1‑81 (MQ=255) aTCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCAc                                                                         <  2:508859/80‑1 (MQ=255)  tCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCt                                                                       <  1:508860/81‑1 (MQ=255)    gTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCttt                                                                     >  2:508903/1‑81 (MQ=255)       gaagaTTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGt                                                                   <  1:508861/80‑1 (MQ=255)         agaTTTCGGTGCACGGGCAGCACGGACCAGTGTAACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTAt                                                                <  1:508862/81‑1 (MQ=255)           aTTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCg                                                              <  1:508863/81‑1 (MQ=255)            tttCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGc                                                             <  1:508864/81‑1 (MQ=255)            tttCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGc                                                             >  1:508919/1‑81 (MQ=255)              tCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGcc                                                            <  1:508865/80‑1 (MQ=255)                 gTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATg                                                        <  1:508866/81‑1 (MQ=255)                 gTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATg                                                        <  2:508867/81‑1 (MQ=255)                  tGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGc                                                       <  1:508868/81‑1 (MQ=255)                   gCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGc                                                       >  1:508933/1‑80 (MQ=255)                    cACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCg                                                      <  1:508869/80‑1 (MQ=255)                    cACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGa                                                     <  2:508870/81‑1 (MQ=255)                       gggCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGaat                                                   >  1:508917/1‑80 (MQ=255)                         gCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGaataa                                                 >  2:508945/1‑80 (MQ=255)                          ccGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAataata                                               <  1:508871/81‑1 (MQ=255)                           cGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATAc                                              <  2:508872/81‑1 (MQ=255)                             aaCGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGt                                            <  1:508873/80‑1 (MQ=255)                              aCGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGtt                                           <  1:508874/81‑1 (MQ=255)                                  aCCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTcgcg                                       >  2:508938/1‑81 (MQ=255)                                      gtgtCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCggg                                   <  2:508875/81‑1 (MQ=255)                                         tCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGAt                                 <  2:508876/80‑1 (MQ=255)                                         tCACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATc                                >  1:508940/1‑81 (MQ=255)                                          cACCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATcc                               >  2:508926/1‑81 (MQ=255)                                           aCCCATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATccc                              <  1:508877/81‑1 (MQ=255)                                            cccATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCt                             <  1:508878/81‑1 (MQ=255)                                              cATCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTAc                           >  1:508925/1‑81 (MQ=255)                                               aTCTGCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTAct                          >  1:508937/1‑81 (MQ=255)                                                   gCCAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTActtctt                      >  2:508939/1‑81 (MQ=255)                                                     cAGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCtttt                    >  2:508944/1‑81 (MQ=255)                                                      aGAAGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCttttt                   <  2:508879/81‑1 (MQ=255)                                                        aaGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTcc                 <  2:508880/81‑1 (MQ=255)                                                         aGTTGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTccc                <  1:508881/81‑1 (MQ=255)                                                           ttGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTTCTTCTTTTTCCCAg              <  2:508882/81‑1 (MQ=255)                                                            tGTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTcc                 >  2:508920/1‑77 (MQ=255)                                                             gTCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGAt            >  2:508927/1‑81 (MQ=255)                                                              tCAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGAtt           >  2:508928/1‑81 (MQ=255)                                                               cAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGAttt          >  2:508934/1‑81 (MQ=255)                                                               cAGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACATCTTTTTCCCAGAttt          >  2:508946/1‑81 (MQ=255)                                                                aGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGAttt          <  1:508883/80‑1 (MQ=255)                                                                aGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGATTTc         <  1:508885/81‑1 (MQ=255)                                                                aGATGCGTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGATTTc         <  1:508884/81‑1 (MQ=255)                                                                      gTATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGATTTCGTAGGc   <  1:508886/81‑1 (MQ=255)                                                                       tATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGATTTCGTAGGCt  <  1:508888/81‑1 (MQ=255)                                                                       tATGACGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGATTTCGTAGGCt  <  1:508887/81‑1 (MQ=255)                                                                           |                                                                             ATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACGTTCGCGCGGGATCCCTACTTCTTTTTCCCAGATTTCGTAGGCT  >  NC_000913/2821448‑2821598

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: ‑