Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,822,103	A→C	intergenic (‑92/+36)	alaS ← / ← recX	alanyl‑tRNA synthetase/regulatory protein for RecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,822,103	0	A	C	100.0%	164.1 / NA	45	intergenic (‑92/+36)	alaS/recX	alanyl‑tRNA synthetase/regulatory protein for RecA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/27);  total (18/27)

CTATGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAGGTTCG  >  NC_000913/2821959‑2822250                                                                                                                                                 |                                                                                                                                                    cTATGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTc                                                                                                                                               >  2:555414/1‑151 (MQ=255)   aTGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCCCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTc                                                                                                                                               >  1:782467/1‑149 (MQ=255)     ggAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTg                                                                                                                                            <  2:63817/150‑1 (MQ=255)        aaaaGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTg                                                                                                                                        <  1:8808/151‑1 (MQ=255)        aaaaGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTg                                                                                                                                        <  2:865396/151‑1 (MQ=255)         aaaGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTgg                                                                                                                                       <  2:188488/151‑1 (MQ=255)                     cGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATAc                                                                                                                           <  2:584485/151‑1 (MQ=255)                                   cAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGtt                                                                                                                 <  2:686150/147‑1 (MQ=255)                                   cAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGt                                                                                                             >  2:772075/1‑151 (MQ=255)                                    aGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTc                                                                                                            >  2:896416/1‑151 (MQ=255)                                     gCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCg                                                                                                           <  1:478259/151‑1 (MQ=255)                                       ggTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGc                                                                                                         <  1:93168/151‑1 (MQ=255)                                               tGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAttt                                                                                                  <  1:648515/150‑1 (MQ=255)                                               tGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAAttt                                                                                                  <  1:747825/150‑1 (MQ=255)                                                  tCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGcc                                                                                              >  2:497745/1‑151 (MQ=255)                                                   cATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCa                                                                                             <  2:266412/151‑1 (MQ=255)                                                      aaTTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAAt                                                                                          >  2:880720/1‑151 (MQ=255)                                                           tCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTAACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCt                                                                                     >  1:779622/1‑151 (MQ=255)                                                                aaaTCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAGGATACGTTTTCTTGACTGGGAAGAAAAATACCGAATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                                                                                 >  2:645684/1‑150 (MQ=255)                                                                aaaTCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGat                                                                                 >  1:492105/1‑150 (MQ=255)                                                                aaaTCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGCCGGCAAAATTTCGCCAAATCTCCTGGat                                                                                 >  1:254993/1‑150 (MQ=255)                                                                 aaTCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGatat                                                                               <  2:112949/151‑1 (MQ=255)                                                                     aaGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTc                                                                           <  1:449730/151‑1 (MQ=255)                                                                      aGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTcc                                                                          <  1:112209/151‑1 (MQ=255)                                                                      aGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTcc                                                                          <  1:779386/151‑1 (MQ=255)                                                                          aaCGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCAt                                                                        <  2:257221/149‑1 (MQ=255)                                                                          aaCGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATc                                                                       >  2:262094/1‑150 (MQ=255)                                                                              aaaTATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGAt                                                                   <  1:838643/150‑1 (MQ=255)                                                                               aaTATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGAt                                                                   >  2:249729/1‑149 (MQ=255)                                                                                atatCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGc                                                                <  1:107179/151‑1 (MQ=255)                                                                                atatCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGc                                                                >  2:366950/1‑151 (MQ=255)                                                                                  atCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCa                                                              <  2:840328/151‑1 (MQ=255)                                                                                                    gtgtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACg                                             <  1:863805/150‑1 (MQ=255)                                                                                                      gtCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTg                                          <  1:728450/151‑1 (MQ=255)                                                                                                       tCTAAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCt                                           <  2:782487/149‑1 (MQ=255)                                                                                                        cTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTgg                                         >  1:641576/1‑150 (MQ=255)                                                                                                          tAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGAt                                       <  1:136223/150‑1 (MQ=255)                                                                                                               gccgGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAAc                                 <  1:853639/151‑1 (MQ=255)                                                                                                                ccgGGCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACt                                <  2:417605/151‑1 (MQ=255)                                                                                                                    gCTGGTAACTGAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACttttt                            >  1:936025/1‑151 (MQ=255)                                                                                                                      tGGTAACTGAAAAGCGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCt                          >  1:687989/1‑151 (MQ=255)                                                                                                                              gAAAAGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGAc                  <  1:682745/151‑1 (MQ=255)                                                                                                                               aaaaGTGGGAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGAc                  >  2:681712/1‑150 (MQ=255)                                                                                                                                     gggAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGc           >  1:795788/1‑151 (MQ=255)                                                                                                                                       gAATAAGATACGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCaa         <  1:549220/151‑1 (MQ=255)                                                                                                                                           aagatacgTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAGGt     >  1:159952/1‑151 (MQ=255)                                                                                                                                            agatacgTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAGGtt    <  1:656016/151‑1 (MQ=255)                                                                                                                                              atacgTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAGGTTCg  >  2:918400/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CTATGGAAAAAGTCGAGAAACGCCTGACGGATCTCAGCGGTGCTCTTGCTCATAATTATCCTGAAATCAAGCTAACGAAATATCGCCACCAGCTCCAGCGTGTCTTAACCGCCGGGCTGGTAACTGAAAAGTGGGAATAAGATAAGTTTTCTTGACTGGGAAGTAAAATACCGTATGCGTTCAGTCGGCAAAATTTCGCCAAATCTCCTGGATATCTTCCATCAGATAGCCACGATAGAGCAGAAAACGCTGGATCTTAACTTTTTCTGAAAAGACAGTTGGCAAAGGTTCG  >  NC_000913/2821959‑2822250

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