Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	720,946	C→T	intergenic (+106/+110)	ybfK → / ← kdpE	uncharacterized protein/response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	720,946	0	C	T	100.0%	198.8 / NA	57	intergenic (+106/+110)	ybfK/kdpE	uncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/38);  total (19/38)

AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGC  >  NC_000913/720849‑721038                                                                                                  |                                                                                              aaGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAg                                                                                             <  2:931389/100‑1 (MQ=255)  aGGTGAAGATATTCTAACCCGGAAATAAGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTCAATTTATTTAATGATGAGCTTTTTACTTAGt                                                                                            <  1:317025/100‑1 (MQ=255)   ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTa                                                                                           >  2:1466667/1‑100 (MQ=255)      gAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAata                                                                                        <  1:271296/100‑1 (MQ=255)       aaGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAatat                                                                                       <  2:1134887/100‑1 (MQ=255)            atTCTATACAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                  <  2:842430/100‑1 (MQ=255)               cTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑Att                                                                              >  1:293548/1‑100 (MQ=255)                tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTg                                                                             <  1:120307/100‑1 (MQ=255)                    cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAACATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAAtt                                                                         >  2:978884/1‑100 (MQ=255)                    cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAAtt                                                                         >  2:177469/1‑100 (MQ=255)                    cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAAt                                                                          <  2:253113/99‑1 (MQ=255)                     aGGAAATACTGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTg                                                                        <  1:672882/100‑1 (MQ=255)                      ggAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGt                                                                       <  1:718408/100‑1 (MQ=255)                       gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTg                                                                        <  2:970791/98‑1 (MQ=255)                       gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGtt                                                                      <  2:1073861/100‑1 (MQ=255)                            aCGGTTGAACTGTGAAGGTTATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTAtgtt                                                                 <  1:1264214/100‑1 (MQ=255)                              ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTtg                                                               <  2:482782/100‑1 (MQ=255)                                ttGAACTGTGAAGGTGATTAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTtgtg                                                             <  2:443147/100‑1 (MQ=255)                                     cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGtgttg                                                        <  2:1300197/100‑1 (MQ=255)                                      tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGtgttgt                                                       <  2:107769/100‑1 (MQ=255)                                       gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAATAT‑ATTGAATTGTTATTTTTGTGttttg                                                        >  2:1027876/1‑98 (MQ=255)                                        tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGttt                                                     >  1:1178794/1‑100 (MQ=255)                                           aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                  <  1:1226143/100‑1 (MQ=255)                                            ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                 <  1:298378/100‑1 (MQ=255)                                              taataaaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTa                                                    <  2:760068/90‑1 (MQ=255)                                               gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                              >  2:531261/1‑100 (MQ=255)                                                 tcggtataaacTGGTCAAGCTTAATTTATTTAATGCTGAGCTCTTTTTATAGTAATATAAAATGATTG‑ATTGTACATTTTGTGCTATGTTGAAAATaaaa                                            >  2:602507/8‑100 (MQ=255)                                                  cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                           >  1:697601/1‑100 (MQ=255)                                                  cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                           >  1:985527/1‑100 (MQ=255)                                                  cAATAAAAAATGATCAATCTTAATTTATTTAATGATCAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                           >  1:1114927/1‑100 (MQ=255)                                                   aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                          >  1:595367/1‑100 (MQ=255)                                                   aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCATTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                           >  2:608200/1‑99 (MQ=255)                                                    aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAgn                                         <  1:964096/100‑2 (MQ=255)                                                    aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                         <  1:281427/100‑1 (MQ=255)                                                     tAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGcc                                        >  1:1032826/1‑100 (MQ=255)                                                       aaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGt                                      <  1:985099/100‑1 (MQ=255)                                                        aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                     >  1:929947/1‑100 (MQ=255)                                                             gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                <  1:1030814/100‑1 (MQ=255)                                                             gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                <  1:1405047/100‑1 (MQ=255)                                                               tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                              <  2:1290681/100‑1 (MQ=255)                                                                cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                             <  2:33913/100‑1 (MQ=255)                                                                cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                             <  1:188093/100‑1 (MQ=255)                                                                   tCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                          <  2:293548/100‑1 (MQ=255)                                                                   tCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                          >  1:96970/1‑100 (MQ=255)                                                                         tttatttaATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                    <  2:1267829/100‑1 (MQ=255)                                                                          ttatttaATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                   <  2:824694/100‑2 (MQ=38)                                                                           tatttaATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                  >  1:725420/1‑100 (MQ=255)                                                                           tatttaATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                  <  1:1262593/100‑1 (MQ=255)                                                                           tatttaATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttt                  <  2:730365/100‑3 (MQ=37)                                                                            atttaATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTCTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                 <  1:1492198/100‑1 (MQ=255)                                                                             tttaATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGTTttaatt                <  1:829722/100‑1 (MQ=255)                                                                               taATGATGAGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa              >  2:684637/1‑100 (MQ=255)                                                                                  tgatgaGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaataa           <  1:740851/100‑1 (MQ=255)                                                                                    atgaGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa         <  1:706405/100‑1 (MQ=255)                                                                                      gaGCTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg       <  1:1250796/100‑1 (MQ=255)                                                                                         cTTTTTACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag    >  1:384149/1‑100 (MQ=255)                                                                                           ttttACTTAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAATaggcaggt  <  2:885253/100‑2 (MQ=255)                                                                                                  |                                                                                              AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAAT‑ATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGC  >  NC_000913/720849‑721038

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