Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,710,250	G→A	intergenic (‑86/‑170)	rseD ← / → nadB	rpoE leader peptide/L‑aspartate oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,710,250	0	G	A	100.0%	117.7 / NA	35	intergenic (‑86/‑170)	rseD/nadB	rpoE leader peptide/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/19);  total (16/19)

ATCCGAGGTAAAGTCTCCCCAAACCAAATTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGATTATCGGTAGCGGCGCAGCCGGACTTTCACTGGCGCTACG  >  NC_000913/2710011‑2710496                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       aTCCGAGGTAAAGTCTCCCCAAACCAAATTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTAtt                                                                                                                                                                                                                                             >  1:502977/1‑251 (MQ=255)    cGAGGTAAAGTCTCCCCAAACCAAATTTCCACGCGCTATCGAAACGCCACGCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTGTTTGTAAAGTTCCGTAACGCATCTGTTTTTGGGTGTCATGCTTTAGACGGATCATTATCTATCAGTATAATTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAAGGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTAATTTca                                                                                                                                                                                                                                          <  1:345230/251‑1 (MQ=255)                 ccccAAACCAAATTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAAtcctc                                                                                                                                                                                                                             >  1:66844/1‑251 (MQ=255)                        ccaaaTTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCAccagc                                                                                                                                                                                                                      <  1:533747/251‑1 (MQ=255)                          aaaTTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCAccagcc                                                                                                                                                                                                                     >  2:148455/1‑250 (MQ=255)                              ttCCACGCGCTATCGAAACTCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATGATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACGGCCTGAAGAGTAACCCAACAGCCTTTTTATTGCACCACCTAATCCTCCACCAGCCAGTa                                                                                                                                                                                                                 <  2:473718/250‑1 (MQ=255)                                ccACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACt                                                                                                                                                                                                              >  1:385992/1‑251 (MQ=255)                                   cgcgCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct                                                                                                                                                                                                           >  1:130942/1‑251 (MQ=255)                                      gcTATCGAAACGCCACTCCATTAGCGGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTATAAATGCCCATTTAACTGCCTGAAGAGGAACCCAACAGCCTTTTGATTTCACCACCTAATCCCCCACCAGCCAGTAACTTCTCtt                                                                                                                                                                                                        <  2:252154/251‑1 (MQ=255)                                                    aCTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTATCTTCTCTTTTTCTCGCCGCCCt                                                                                                                                                                                          <  1:410266/251‑1 (MQ=255)                                                            aGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGc                                                                                                                                                                                  >  2:404597/1‑251 (MQ=255)                                                                gCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCGCGCCGCCCTGCGTCAGCgtg                                                                                                                                                                               <  2:448214/250‑1 (MQ=255)                                                                           aagcaaAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGACATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAAAATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGTCCTGCGTCAGCGTGTTTAGCAACTGt                                                                                                                                                                   <  2:469572/251‑1 (MQ=255)                                                                              caaAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAAc                                                                                                                                                                <  2:294688/251‑1 (MQ=255)                                                                               aaaGGGTTAGAGTGTCTAGTTTTTGTAAAGTTCCGTAACGCAGCTGATTTTGTTTGTCATGCTGTAGACGGCTCATTATCTATCATTATAAGTATACAGGATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGCAGCGTAACCAAACAGCCTTTTTATTTCACCACCTCATCATCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGAGTTTAGCAACTGtcac                                                                                                                                                                >  2:490621/1‑247 (MQ=255)                                                                                          gtgtCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAAt                                                                                                                                                    <  1:138105/251‑1 (MQ=255)                                                                                                       gTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTTTCTAGCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGGACCCCAACAGCCTTTTTATGTCCCCGCCTAATCCTCCGCCAGCCAGTAACGTCTCTTTTTCTCTCCGCCCTGCGTCAGCGTGTTTAGGAGCTGTAACAAATATTAAAATAGCAGGTGTTTAt                                                                                                                                       <  2:347335/251‑1 (MQ=255)                                                                                                                         tCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGatgct                                                                                                                      >  2:173255/1‑250 (MQ=255)                                                                                                                            gTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCCCCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGatgct                                                                                                                      <  2:436317/247‑1 (MQ=255)                                                                                                                              ttttGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGatgctatgct                                                                                                                 <  2:299541/250‑1 (MQ=255)                                                                                                                                    ttgtCATGCTGTAGACGGATCATTATCTTTCATTATAATTTTAAAGATTCTGAACATCCCATTATCTTGGTAGAAATGCCCATTTACCTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCCCCTAATCCCCCTCCAGCCAGTAACTTCTCTTTT‑CTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCaaa                                                                                                          <  1:372662/250‑1 (MQ=255)                                                                                                                                             tGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGttt                                                                                                  >  2:115545/1‑250 (MQ=255)                                                                                                                                                                      ataAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAAt                                                                         <  1:336704/250‑1 (MQ=255)                                                                                                                                                                        aaGTCTTCAGAATCTGAACCTCGCATTAACTGGGTAGAAATTCCCCTTTAACTGCCTAAAGAGTAACCCAACAGCCTTTTTATTTCACCCCCTAATCCCCCACCAGCCAGTAACTTCGCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGATTATACGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATAAAACAAAGAAAATGAAtccc                                                                      <  2:49887/251‑4 (MQ=255)                                                                                                                                                                                             cGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGAc                                                 >  1:358429/1‑251 (MQ=255)                                                                                                                                                                                                      tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCTTCAGCGTGTTGAGCAACTGTAACAAATATTAAAATATCAGgtg                                                                                                                                            <  1:413263/151‑1 (MQ=255)                                                                                                                                                                                                      tgtgTAGAAATGCCCATTTAACTGACTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAgggg                                                                                                                                            >  2:413263/1‑149 (MQ=255)                                                                                                                                                                                                       gtgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGAtt                                        >  1:444789/1‑250 (MQ=255)                                                                                                                                                                                                              aaTGCCCATTTAACTGCCTGACGAGTAACCCAACAGCCTTTTTATTTCACCACATAATCATCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGATTATCGGt                                  >  2:537400/1‑249 (MQ=255)                                                                                                                                                                                                              aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCGGAACATTCATGTGACGTGTTGATTATCGGTAg                                <  2:385992/251‑1 (MQ=255)                                                                                                                                                                                                                         aaCTGCCTGAAGAGTAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGATTATCGGTAGCGGCGCAGCCg                     >  2:462123/1‑251 (MQ=255)                                                                                                                                                                                                                                       tAACCCAACAGCCTTTTTATTTCACCACCTACTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAAGTAATCAAAGAAAATGAATACTCTCCCTGTACATTCAt                                                      <  2:327212/204‑1 (MQ=255)                                                                                                                                                                                                                                       tAACCCAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGACAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTACACAAAGAAAATGAATACTCTCCCTGAACATTCAt                                                      >  1:327212/1‑204 (MQ=255)                                                                                                                                                                                                                                          cccAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGATTATCGGTAGCGGCGCAGCCGGACTTTCACTGGCGCTa    <  2:444789/251‑1 (MQ=255)                                                                                                                                                                                                                                            cAACAGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGATTATCGGTAGCGGCGCAGCCGGACTTTCACTGGCGCTACg  >  2:140601/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       ATCCGAGGTAAAGTCTCCCCAAACCAAATTTCCACGCGCTATCGAAACGCCACTCCATTAGCTGCAATTTGAGCAAGCAAAGGGTTAGAGTGTCTCGTTTTTGTAAAGTTCCGTAACGCATCTGTTTTTGTTTGTCATGCTGTAGACGGATCATTATCTATCATTATAAGTCTACAGAATCTGAACATCGCATTATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTAAAATAGCAGGTGTTTATCCGCACAACATGATGCTATGCTGACCAAACCATGTTTAGTAAATTAAACAAAGAAAATGAATACTCTCCCTGAACATTCATGTGACGTGTTGATTATCGGTAGCGGCGCAGCCGGACTTTCACTGGCGCTACG  >  NC_000913/2710011‑2710496

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