Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,159,997	A→T	L265Q (CTG→CAG)	sthA ←	soluble pyridine nucleotide transhydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,159,997	0	A	T	100.0%	112.2 / NA	33	L265Q (CTG→CAG)	sthA	soluble pyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (16/17);  total (16/17)

GGCTTCGCCTTTTACCAGCGCCTGCGCGGCAATGCGCCCCTGGTCATAGGCCGCCGACGCCAGGCTCGGATAACCAATCACGTCGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCA  >  NC_000913/4159777‑4160239                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   ggCTTCGCCTTTTACCAGCGCCTGCGCGGCAATGCGCCCCTGGTCATAGGCCGCCGACGCCAGGCTCGGATAACCAATCACGTCGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACt                                                                                                                                                                                                                      <  1:229233/251‑1 (MQ=255)    ttCGCCTTTTACCAGCGCCTGCGCGGCAATGCGCCCCTGGTCATAGGCCGCCGACGCCAGGCTCGGATAACCAATCACGTCGCCCACCGCGTAAACGAGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCCGCTTTCAGTTTTTTACCCGACTTc                                                                                                                                                                                                                    >  1:51070/1‑250 (MQ=255)                     ccTGCGCGGCAATGCGCCCCTGGTCATAGGCCGCCGACGCCAGGCTCGGATAACCAATCACGTCGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCa                                                                                                                                                                                                  <  1:180004/251‑1 (MQ=255)                     ccTGCGCGGCAATGCGCCCCTGGTCATAGGCCGCCGACGCCAGGCTCGGATAACCAATCACGTCGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGAATTCAGATGCATGATCACAcc                                                                                                                                                                                                   >  1:262615/1‑250 (MQ=255)                      cTGCGCGGCAATGCGCCCCTGGTCATAGGCCGCCGACGCCAGGCTCGGATAACCAATCACGTCGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCa                                                                                                                                                                                                  >  2:180004/1‑250 (MQ=255)                                        ctgatCATAGGCCGCCGACGCCAGGCTCGGATCACCAATCACGTAGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTCACCGGCTTCAGATGCAAGATCACACCATCGTCACAGCCTTCGAtct                                                                                                                                                                               >  2:133038/5‑251 (MQ=255)                                        cTGGTCATAGGCCGCCGACGCCAGGCTCGGATAACCAATCACGTCGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGAtct                                                                                                                                                                               >  1:14488/1‑251 (MQ=255)                                                       cgACGCCAGGCTCGGATAACCAATCACGTCGCCCCCCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGt                                                                                                                                                                >  1:78803/1‑251 (MQ=255)                                                                       tAACCAATCACGTCGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGCATCACTAc                                                                                                                                                >  1:77453/1‑251 (MQ=255)                                                                           caagcactacGCCCACCGCGTAAACGTGTGGCTGTGCGGACTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGTCTGTAAGTCTCTAGCTCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTAGCAGTCAGCTGTCAGTTTTTTGCCCGACTTCAGATGCATGATCTCACCAACGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAAGCACTACGcta                                                                                                                                            <  2:189055/242‑3 (MQ=255)                                                                                   ccGCCCACCGCGTAAACGGGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTGCGTTGTGACGAATCACTACGCCACTGTTCCa                                                                                                                                    <  2:14488/250‑1 (MQ=255)                                                                                                    gtgtGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGGGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTAACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGCATCACTACGCCACTGTTCCAGGAGTTATAGgagag                                                                                                                     >  2:145525/1‑249 (MQ=255)                                                                                                        ggCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGGGAGAGAATc                                                                                                               >  2:3577/1‑251 (MQ=255)                                                                                                                ggTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGTTAGGAGAGAGAATCTGACAtct                                                                                                       >  1:95743/1‑251 (MQ=255)                                                                                                                               gTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAAt                                                                                         <  2:77453/250‑1 (MQ=255)                                                                                                                                  gACCTTCAGCTTGTC‑GCGGCTGTCAGTTGCTAGCCCAATGTTCTGTAACGCCTGCTTATCGGTCTTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTGTCAGTGTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTCCTCTTCGTTGTGACGAATCACTCCGCTACTGTGCAAGAAGTGATAGGAGATAGAATCTGACATCTCCTGATCGAGAACTGCca                                                                                     <  2:43760/251‑1 (MQ=255)                                                                                                                                     cTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCcagc                                                                                   >  2:259850/1‑250 (MQ=255)                                                                                                                                       tCAGCT‑GGCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCcagca                                                                                  <  1:133038/249‑1 (MQ=255)                                                                                                                                                   gggCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTGTCAGTTTGTTACCCGACTTCAGAGGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgt                                                                     <  1:175273/250‑1 (MQ=255)                                                                                                                                                        gTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGtt                                                                  <  1:137948/249‑1 (MQ=255)                                                                                                                                                                       aaTGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                  >  1:240819/1‑250 (MQ=255)                                                                                                                                                                               ggaaCGCCAGCGAATCGGTATTGCCGGTGCGGCCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCCCAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTTATCGAGATATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCCCTTTTACATCCAt                                          <  2:17855/248‑1 (MQ=255)                                                                                                                                                                                 aaCGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAc                                        <  2:51070/250‑1 (MQ=255)                                                                                                                                                                                             tCGGTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCg                            >  1:56684/1‑250 (MQ=255)                                                                                                                                                                                              cGGTATTACCTGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGAGCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGCAGATCGAc                          <  2:53266/251‑1 (MQ=255)                                                                                                                                                                                               ggTATTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACg                         >  2:182110/1‑251 (MQ=255)                                                                                                                                                                                                gTATTACCGGTTCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGc                        <  2:232909/251‑1 (MQ=255)                                                                                                                                                                                                  aTTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                      <  1:128933/251‑1 (MQ=255)                                                                                                                                                                                                  aTTACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCa                       <  1:224656/250‑1 (MQ=255)                                                                                                                                                                                                   ttACCGGTGCGACCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCCGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                      >  1:169022/1‑250 (MQ=255)                                                                                                                                                                                                              aCCGTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCGTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGa          <  2:29051/251‑1 (MQ=255)                                                                                                                                                                                                                 gTTGGCATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc        >  2:147676/1‑250 (MQ=255)                                                                                                                                                                                                                      cATAGAGCTGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCa  <  1:182110/251‑1 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   GGCTTCGCCTTTTACCAGCGCCTGCGCGGCAATGCGCCCCTGGTCATAGGCCGCCGACGCCAGGCTCGGATAACCAATCACGTCGCCCACCGCGTAAACGTGTGGCTGTGCGGTCTGATACATGCTGTTGACCTTCAGCT‑GTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCA  >  NC_000913/4159777‑4160239

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