Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,157,712	(C)6→7	coding (63/927 nt)	yqhG →	DUF3828 family putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,157,706	1	.	C	100.0%	194.3 / NA	54	Q19H (CAG→CAC)	yqhG	DUF3828 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (25/29);  total (25/29)

ATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑CCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGG  >  NC_000913/3157570‑3157802                                                                                                                                          |                                                                                                 atcatcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATCATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑cc                                                                                                <  1:161509/139‑1 (MQ=255) atcatcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑cc                                                                                                <  2:147295/139‑1 (MQ=255)   catcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑cccc                                                                                              >  2:275882/1‑139 (MQ=255)     tcatGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAA‑cccccc                                                                                            >  1:214300/1‑139 (MQ=255)        tGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGAGGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑ccccccctc                                                                                         >  2:548077/1‑136 (MQ=255)           agagCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTcaca                                                                                       >  2:39301/1‑139 (MQ=255)           agagCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTcaca                                                                                       <  1:705577/139‑1 (MQ=255)           agagCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTcaca                                                                                       <  1:43444/139‑1 (MQ=255)             agCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAg                                                                                                   >  2:445770/1‑125 (MQ=255)             agCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAg                                                                                                   <  1:445770/125‑1 (MQ=255)             agCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGa                                                                                     <  1:13686/139‑1 (MQ=255)             agCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGa                                                                                     >  1:980021/1‑139 (MQ=255)                 ggAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGt                                                                                 <  2:541358/139‑1 (MQ=255)                         ttACAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaacaaa                                                                         <  2:651752/139‑1 (MQ=255)                          tAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaacaaac                                                                        <  2:214300/139‑1 (MQ=255)                                aacaacTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGGCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCg                                                                  <  1:548077/139‑1 (MQ=255)                                aacaacTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGAc                                                                                    >  1:583459/1‑121 (MQ=255)                                aacaacTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGAc                                                                                    <  2:583459/121‑1 (MQ=255)                                  caacTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGc                                                                <  2:980021/139‑1 (MQ=255)                                     cTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaacaaaca                                                                       >  1:213520/1‑129 (MQ=255)                                     cTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGaacaaaca                                                                       <  2:213520/129‑1 (MQ=255)                                     cTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGc                                                                <  1:580089/136‑1 (MQ=255)                                     cTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGc                                                                >  2:580089/1‑136 (MQ=255)                                             tatCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                     >  2:38730/1‑139 (MQ=255)                                             tatCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                     <  1:39301/139‑1 (MQ=255)                                                  gggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGAcc                                                                                   <  2:461593/104‑1 (MQ=255)                                                  gggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGAcc                                                                                   >  1:461593/1‑104 (MQ=255)                                                  gggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACta                                                >  1:670486/1‑139 (MQ=255)                                                   ggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCg                                                                  >  2:442395/1‑120 (MQ=255)                                                   ggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCg                                                                  <  1:442395/120‑1 (MQ=255)                                                   ggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                     <  1:1000868/133‑1 (MQ=255)                                                   ggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                     >  2:1000868/1‑133 (MQ=255)                                                   ggATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtat                                               >  2:922862/1‑139 (MQ=255)                                                    gATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCg                                                                  <  2:458207/119‑1 (MQ=255)                                                    gATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCg                                                                  >  1:458207/1‑119 (MQ=255)                                                    gATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTAt                                                        >  2:949952/1‑129 (MQ=255)                                                    gATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTAt                                                        <  1:949952/129‑1 (MQ=255)                                                    gATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtata                                              >  1:426433/1‑139 (MQ=255)                                                        aaCGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCa                                                               <  2:327940/118‑1 (MQ=255)                                                        aaCGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCa                                                               >  1:327940/1‑118 (MQ=255)                                                        aaCGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATc                                          <  2:6679/139‑1 (MQ=255)                                                                  agagGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACt                                <  2:670486/139‑1 (MQ=255)                                                                     ggTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGcc                             <  1:38730/139‑1 (MQ=255)                                                                            aaCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtat                                               >  2:90425/1‑114 (MQ=255)                                                                            aaCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtat                                               <  1:90425/114‑1 (MQ=255)                                                                                gATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGa                                                                                     <  1:188529/72‑1 (MQ=255)                                                                                gATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGa                                                                                     >  2:188529/1‑72 (MQ=255)                                                                                 aTGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGt                 >  2:829281/1‑139 (MQ=255)                                                                                           taCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGc                                                                <  1:983397/82‑1 (MQ=255)                                                                                           taCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGc                                                                >  2:983397/1‑82 (MQ=255)                                                                                                cTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCgg  >  2:174200/1‑139 (MQ=255)                                                                                                        agcagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTAt                                                        <  2:835614/77‑1 (MQ=255)                                                                                                        agcagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTAt                                                        >  1:835614/1‑77 (MQ=255)                                                                                                                       ttttACCGTACACGCACAGCCCCCCCTCACAGAcc                                                                                   >  2:881324/1‑35 (MQ=255)                                                                                                                       ttttACCGTACACGCACAGCCCCCCCTCACAGAcc                                                                                   <  1:881324/35‑1 (MQ=255)                                                                                                                       ttttACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtata                                              >  1:910235/1‑72 (MQ=255)                                                                                                                       ttttACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACtata                                              <  2:910235/72‑1 (MQ=255)                                                                                                                         ttACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATc                                                       <  2:793312/61‑1 (MQ=255)                                                                                                                         ttACCGTACACGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATc                                                       >  1:793312/1‑61 (MQ=255)                                                                                                                                 cacGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                     <  2:185541/55‑1 (MQ=255)                                                                                                                                 cacGCACAGCCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAg                                                     >  1:185541/1‑55 (MQ=255)                                                                                                                                          |                                                                                                 ATCATCATGCAGAGCCGGAAGCTCTTAAAAGAACAACTCATCTATATCCGGGATAAACGCAACGGAGAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAG‑CCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGG  >  NC_000913/3157570‑3157802

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