Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,296,984	G→A	*192* (TAG→TAA)	yraP →	outer membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,296,984	0	G	A	100.0%	144.7 / NA	45	*192* (TAG→TAA)	yraP	outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/23);  total (22/23)

TGAAAGTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCAGAACAACAGCGCC  >  NC_000913/3296848‑3297113                                                                                                                                         |                                                                                                                                  tGAAAGTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaaca                                                                                                                                 <  1:1000022/139‑1 (MQ=255) tGAAAGTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaaca                                                                                                                                 <  2:1743985/139‑1 (MQ=255)    aaGTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaacagc                                                                                                                               <  1:1205428/138‑1 (MQ=255)    aaGTGACCAACGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaacagc                                                                                                                               >  2:1205428/1‑138 (MQ=255)          ccaccGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaaca                                                                                                                                 <  2:1364291/130‑1 (MQ=255)          ccaccGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaaca                                                                                                                                 >  1:1364291/1‑130 (MQ=255)                                    aTGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGt                                                                                              <  2:1687244/139‑1 (MQ=255)                                      ggggCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGcc                                                                                                                              <  1:1316865/105‑1 (MQ=255)                                      ggggCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGcc                                                                                                                              >  2:1316865/1‑105 (MQ=255)                                        ggCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                          >  2:1555788/1‑139 (MQ=255)                                        ggCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                          >  1:1344504/1‑139 (MQ=255)                                         gCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                                   >  1:836565/1‑129 (MQ=255)                                         gCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                                   <  2:836565/129‑1 (MQ=255)                                         gCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         <  2:1151713/139‑1 (MQ=255)                                         gCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         <  2:386920/139‑1 (MQ=255)                                         gCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         <  2:264032/139‑1 (MQ=255)                                         gCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         >  1:1317424/1‑139 (MQ=255)                                                         gAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATa                                                                         <  2:291824/139‑1 (MQ=255)                                                          aaGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaa                                                                              <  2:122412/133‑1 (MQ=255)                                                          aaGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaa                                                                              >  1:122412/1‑133 (MQ=255)                                                           aGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACa                                                                       <  2:267838/139‑1 (MQ=255)                                                           aGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACa                                                                       <  2:1344504/139‑1 (MQ=255)                                                           aGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACa                                                                       >  2:710238/1‑139 (MQ=255)                                                                aaGCGGCGGCAGATATTGCCAGTCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGcccc                                                                                                                            <  1:644185/81‑1 (MQ=255)                                                                aaGCGGCGGCAGATATTGCCAGTCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGcccc                                                                                                                            >  2:644185/1‑81 (MQ=255)                                                                 aGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGc                                                                                                                   >  1:1032117/1‑89 (MQ=255)                                                                 aGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGc                                                                                                                   <  2:1032117/89‑1 (MQ=255)                                                                 aGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTACAACGATACAAACCGt                                                                 >  2:1057943/1‑139 (MQ=255)                                                                                      ccGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                                   <  1:1640672/84‑1 (MQ=255)                                                                                      ccGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                                   >  2:1640672/1‑84 (MQ=255)                                                                                      ccGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTAcg                                            >  1:1604566/1‑139 (MQ=255)                                                                                         ggTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                                   <  1:1594899/81‑1 (MQ=255)                                                                                         ggTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGc                                                                                                   >  2:1594899/1‑81 (MQ=255)                                                                                         ggTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         <  2:351151/91‑1 (MQ=255)                                                                                         ggTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         >  1:351151/1‑91 (MQ=255)                                                                                          gTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaacag                                                                                                                                >  1:1776638/1‑51 (MQ=255)                                                                                          gTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATaacag                                                                                                                                <  2:1776638/51‑1 (MQ=255)                                                                                           tGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCa                                       <  2:1317424/139‑1 (MQ=255)                                                                                           tGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCGACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCa                                       <  2:1604566/139‑1 (MQ=255)                                                                                            gAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCAt                                      <  2:1660690/139‑1 (MQ=255)                                                                                               cggcgTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCg                                   >  2:1408803/1‑139 (MQ=255)                                                                                                  cgTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACa                                <  1:710238/139‑1 (MQ=255)                                                                                                            gTAACTACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCt                      >  2:873722/1‑139 (MQ=255)                                                                                                                 tACGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAAt                 <  1:1408803/139‑1 (MQ=255)                                                                                                                   cGGCGTTTACGTTTATTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCa               >  2:827774/1‑139 (MQ=255)                                                                                                                       gtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         <  2:1457744/61‑1 (MQ=255)                                                                                                                       gtttacgtttaTTAAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                         >  1:1457744/1‑61 (MQ=255)                                                                                                                                tattaAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                        <  2:1393601/53‑1 (MQ=255)                                                                                                                                tattaAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                        >  1:1393601/1‑53 (MQ=255)                                                                                                                                tattaAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACa                                                                       <  1:749182/70‑1 (MQ=255)                                                                                                                                tattaAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACa                                                                       >  2:749182/1‑70 (MQ=255)                                                                                                                                tattaAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCAGAACAACAGCGcc  >  2:1651151/1‑139 (MQ=255)                                                                                                                                tattaAATAACAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCAGAACAACAGCGcc  >  1:20393/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                  TGAAAGTGACCACCGAAAACGGTGAAGTGTTCCTGATGGGGCTGGTGACTGAACGTGAAGCGAAAGCGGCGGCAGATATTGCCAGCCGGGTGAGCGGCGTGAAGCGGGTAACTACGGCGTTTACGTTTATTAAATAGCAGCCCCTTGTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAACGATACAAACCGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATGTTTTCCCTCAAATCAGAACAACAGCGCC  >  NC_000913/3296848‑3297113

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