Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,169,127	G→A	A19T (GCC→ACC)	bhsA →	DUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,169,127	0	G	A	100.0%	103.0 / NA	32	A19T (GCC→ACC)	bhsA	DUF1471 domain‑containing multiple stress resistance outer membrane protein BhsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/15);  total (17/15)

GTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTT  >  NC_000913/1168890‑1169366                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                attCACGAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGAATGTGATACATATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTGGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCATCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTAGCGGc                                                                                                                                                                                                                                    >  2:69321/2‑251 (MQ=255)                                gCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTGACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCcaa                                                                                                                                                                                                     <  1:283065/251‑1 (MQ=255)                                 cTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCcaa                                                                                                                                                                                                     >  1:145143/1‑250 (MQ=255)                                                             tttAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATATATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAAcgc                                                                                                                                                                        >  1:21706/1‑251 (MQ=255)                                                                      attaatGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACaaa                                                                                                                                                               <  2:75672/251‑1 (MQ=255)                                                                                        aaaTATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGa                                                                                                                                                >  2:315185/1‑248 (MQ=255)                                                                                        aaaTATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGa                                                                                                                                                >  2:122694/1‑248 (MQ=255)                                                                                              cTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCt                                                                                                                                       >  2:167888/1‑251 (MQ=255)                                                                                                            aGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTAATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCgg                                                                                                                         >  1:282312/1‑251 (MQ=255)                                                                                                                cacaTCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATg                                                                                                                      <  2:87828/250‑1 (MQ=255)                                                                                                                 acaTCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATgag                                                                                                                    <  1:102852/251‑1 (MQ=255)                                                                                                                 acaTCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATg                                                                                                                      >  1:49758/1‑249 (MQ=255)                                                                                                                               gTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATc                                                                                                      >  2:455993/1‑251 (MQ=255)                                                                                                                                         gTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTAtt                                                                                            <  1:167888/251‑1 (MQ=255)                                                                                                                                          ttAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTAATATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGt                                                                                               <  2:421772/247‑1 (MQ=255)                                                                                                                                              gTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTc                                                                                       >  2:206628/1‑251 (MQ=255)                                                                                                                                                aaTTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCACTGCTAACGCGGGGCCAACTCTGGGATCGCTGGACGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCACAATCTTTCCGTATTACTTCTg                                                                                     >  1:423534/1‑251 (MQ=255)                                                                                                                                                          gaagTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAAGCAGTGCTAACGCGGGGACAAATCTGGTATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTcc                                                                           <  2:282312/248‑1 (MQ=255)                                                                                                                                                          gACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTcc                                                                           >  2:119176/1‑251 (MQ=255)                                                                                                                                                              tcattcatCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGATGTTCAGCCAACGCCAGAAGGCCAACAAAAAGTCGGCACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGATGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAACATATTTCCGTATTACTTCTGTAACCGGTCCGAAt                                                                       >  2:23496/1‑251 (MQ=255)                                                                                                                                                                            aaGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAAc                                                         <  1:122694/251‑1 (MQ=255)                                                                                                                                                                                ccACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAAcagca                                                     <  2:145143/251‑1 (MQ=255)                                                                                                                                                                                     attatGAAAAACGTAAAAACCCTCAGCGCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTaa                                                 <  1:19125/250‑1 (MQ=255)                                                                                                                                                                                                  tAAAAACCCTCATCGCTGCGGAGATTTTAAGCTCCATGTAATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAGAAGTCGGTACAATCAGTGCTAACGCGGGGAAAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCATCCATGGCACCACAGTAATTTATAAATAAga                                    >  2:190048/1‑249 (MQ=255)                                                                                                                                                                                                                gCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGc                                                                                                                                                                              >  1:456519/1‑98 (MQ=255)                                                                                                                                                                                                                gCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGc                                                                                                                                                                              <  2:456519/98‑1 (MQ=255)                                                                                                                                                                                                                gCTGCGGCGATTTTAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTTGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCCCATTaa                     <  2:423534/251‑1 (MQ=255)                                                                                                                                                                                                                          ttttAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTAc           >  1:235949/1‑251 (MQ=255)                                                                                                                                                                                                                           tttAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTAc           <  2:235949/250‑1 (MQ=255)                                                                                                                                                                                                                           tttAAGCTCCATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTAc           <  1:354491/250‑1 (MQ=255)                                                                                                                                                                                                                                   ccATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGAttt  >  1:385178/1‑251 (MQ=255)                                                                                                                                                                                                                                   ccATGTCATTTACCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGAt    <  2:275704/249‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                GTTCACCAGTCCAGATCCCATAAAAATAATTGCTTTCTATTTAACTGAAATTTAAAGATTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAAAACCCTCATCGCTGCGGCGATTTTAAGCTCCATGTCATTTGCCAGCTTTGCGGCTGTCGAAGTTCAGTCAACGCCAGAAGGCCAACAAAAAGTCGGTACAATCAGTGCTAACGCGGGGACAAATCTGGGATCGCTGGAAGAGCAGCTGGCGCAAAAAGCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTT  >  NC_000913/1168890‑1169366

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