Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,533,728	T→C	intergenic (+348/‑36)	cysK → / → ptsH	O‑acetylserine sulfhydrylase A/phosphocarrier protein HPr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,533,728	0	T	C	100.0%	177.1 / NA	50	intergenic (+348/‑36)	cysK/ptsH	O‑acetylserine sulfhydrylase A/phosphocarrier protein HPr
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/23);  total (27/23)

GTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGTTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCTAAGGGCTT  >  NC_000913/2533594‑2533849                                                                                                                                       |                                                                                                                          ttttAATGACGTACGAACCGTCAGCGGTCAACCCCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCcac                                                                                                                      <  2:159522/139‑1 (MQ=255)     aaTGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCacaaca                                                                                                                  >  1:1115271/1‑140 (MQ=255)        gACGTACGAAACGTCAGCGGTCAACCCCCGCCAGCAATCGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTa                                                                                                               <  2:1033616/140‑1 (MQ=255)          cGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCca                                                                                                                       >  2:749610/1‑130 (MQ=255)          cGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCca                                                                                                                       <  1:749610/130‑1 (MQ=255)          cGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAacacctgt                                                                                                             >  1:829724/1‑136 (MQ=255)            tACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAATCGCCAGGCTAGACTTTAGCTCCACAACACTAAAcc                                                                                                           >  2:538434/1‑140 (MQ=255)            tACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAAcc                                                                                                           >  2:1060779/1‑140 (MQ=255)            tACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAAcc                                                                                                           >  2:354611/1‑140 (MQ=255)            tACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAAcc                                                                                                           >  1:843592/1‑140 (MQ=255)            tACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAAcc                                                                                                           >  1:1164169/1‑140 (MQ=255)            tACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAAcc                                                                                                           >  1:719141/1‑140 (MQ=255)            tACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAAcc                                                                                                           >  1:636181/1‑140 (MQ=255)                   cGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGt                                                                                                    >  1:273702/1‑140 (MQ=255)                        gCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTgggg                                                                                               >  1:1215589/1‑140 (MQ=255)                        gCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTgggg                                                                                               >  2:339740/1‑140 (MQ=255)                        gCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTgggg                                                                                               >  2:576594/1‑140 (MQ=255)                          ggTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGaa                                                                                             <  1:538434/140‑1 (MQ=255)                             caccccccGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCTTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGAC‑TTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATAc                                                                                         <  1:719268/135‑1 (MQ=255)                               acacCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACa                                                                                        <  1:576594/140‑1 (MQ=255)                                      ccAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAAt                                                                                      >  1:103449/1‑135 (MQ=255)                                      ccAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAAt                                                                                      <  2:103449/135‑1 (MQ=255)                                          cAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCa                                                                             <  2:273702/140‑1 (MQ=255)                                          cAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCa                                                                             >  2:1200397/1‑140 (MQ=255)                                            aTGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCaag                                                                           <  2:814073/140‑1 (MQ=255)                                             tGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCaaga                                                                          >  1:912736/1‑140 (MQ=255)                                                    tATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTAcc                                                                   >  1:668334/1‑140 (MQ=255)                                                              ttCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTc                                                         >  2:139422/1‑140 (MQ=255)                                                              ttCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTc                                                         >  1:1060915/1‑140 (MQ=255)                                                              ttCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTc                                                         >  2:881146/1‑140 (MQ=255)                                                                  tGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGaa                                                     <  1:354611/140‑1 (MQ=255)                                                                   gcgtcgcgtcTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAAc                                                    <  1:1060779/140‑1 (MQ=255)                                                                     gtcgcgtcTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACgg                                                  <  2:912736/140‑1 (MQ=255)                                                                      tcgcgtcTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGt                                                 <  2:636181/140‑1 (MQ=255)                                                                        gcgtcTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCt                                               <  2:668334/140‑1 (MQ=255)                                                                          gtcTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGc                                             >  1:475092/1‑140 (MQ=255)                                                                              gTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCacac                                         <  2:542864/140‑1 (MQ=255)                                                                                       tGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCctgc                                <  1:139422/140‑1 (MQ=255)                                                                                         gcgcTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCctgctg                              >  1:1055805/1‑140 (MQ=255)                                                                                           gcTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGcc                            <  1:339740/140‑1 (MQ=255)                                                                                                       ggCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTaaaa                <  2:225021/140‑1 (MQ=255)                                                                                                        gCTAAAGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAaag               >  1:87860/1‑140 (MQ=255)                                                                                                             aGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGGTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCt          >  1:808380/1‑140 (MQ=255)                                                                                                             aGTCGAACCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCt          <  2:894814/140‑1 (MQ=255)                                                                                                                  aaCCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCTAAggg     >  2:1007020/1‑140 (MQ=255)                                                                                                                  aaCCGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCTAAggg     >  2:721059/1‑140 (MQ=255)                                                                                                                     cGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTTCCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCTAAGGGCtt  <  2:1164169/140‑1 (MQ=255)                                                                                                                     cGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCTAAGGGCtt  <  2:475092/140‑1 (MQ=255)                                                                                                                     cGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCTAAGGGCtt  <  2:1115271/140‑1 (MQ=255)                                                                                                                     cGCCAGGCTAGACTTTAGCTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCTAAGGGCtt  <  1:721059/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                          GTTTAATGACGTACGAAACGTCAGCGGTCAACACCCGCCAGCAATGGACTGTATTGCGCTCTTCGTGCGTCGCGTCTGTTAAAAACTGGCGCTAACAATACAGGCTAAAGTCGAACCGCCAGGCTAGACTTTAGTTCCACAACACTAAACCTATAAGTTGGGGAAATACAATGTTCCAGCAAGAAGTTACCATTACCGCTCCGAACGGTCTGCACACCCGCCCTGCTGCCCAGTTTGTAAAAGAAGCTAAGGGCTT  >  NC_000913/2533594‑2533849

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