Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,673,768	G→A	*141* (TAG→TAA)	yphA →	DoxX family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,673,768	0	G	A	100.0%	183.1 / NA	57	*141* (TAG→TAA)	yphA	DoxX family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/27);  total (30/27)

GTCACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGAT  >  NC_000913/2673632‑2673901                                                                                                                                         |                                                                                                                                      gTCACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGa                                                                                                                                     >  2:1298864/1‑139 (MQ=255)        ttACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGc                                                                                                                              <  1:1433495/139‑1 (MQ=255)                   tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcggcg                                                                                                                   >  2:1362802/1‑139 (MQ=255)                     aCCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTg                                                                                                                 <  1:1298864/139‑1 (MQ=255)                                       ccAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACTGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                                <  2:1053786/122‑1 (MQ=255)                                       ccAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACTGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                                >  1:1053786/1‑122 (MQ=255)                                                aTTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTta                                                                                             <  1:3450/132‑1 (MQ=255)                                                aTTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTta                                                                                             >  2:3450/1‑132 (MQ=255)                                                aTTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAtt                                                                                      >  1:1464648/1‑139 (MQ=255)                                                     tttCTGGAAGAATGTCAGTATCGCTGGTGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATCTAATAAATTACAGc                                                                                 >  2:899385/1‑139 (MQ=255)                                                     tttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAg                                                                                                                                  >  2:303259/1‑90 (MQ=255)                                                     tttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAg                                                                                                                                  <  1:303259/90‑1 (MQ=255)                                                     tttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                 <  2:1464648/139‑1 (MQ=255)                                                               aaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGc                                                                                                                              >  2:395705/1‑84 (MQ=255)                                                               aaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGc                                                                                                                              <  1:395705/84‑1 (MQ=255)                                                               aaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAccc                                                                       >  1:605725/1‑139 (MQ=255)                                                                 tGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcgg                                                                                                                     <  2:1357917/91‑1 (MQ=255)                                                                 tGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcgg                                                                                                                     >  1:1357917/1‑91 (MQ=255)                                                                   tCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTg                                                                                                                 >  1:309943/1‑93 (MQ=255)                                                                   tCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTg                                                                                                                 <  2:309943/93‑1 (MQ=255)                                                                        aTCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGAt                                                                                                                                >  2:536895/1‑73 (MQ=255)                                                                        aTCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGAt                                                                                                                                <  1:536895/73‑1 (MQ=255)                                                                            cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                                          >  1:1303310/1‑139 (MQ=255)                                                                               gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGa                                                       >  2:1205173/1‑139 (MQ=255)                                                                                 gcgTTCTTGCTATTGGCAATTACCGGGCC‑GGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAAc                                                                                                              >  1:446608/1‑81 (MQ=39)                                                                                 gcgTTCTTGCTATTGGCAATTACCGGGCC‑GGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAAc                                                                                                              <  2:446608/81‑1 (MQ=39)                                                                                      cTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGc                                                                                                            <  1:850922/79‑1 (MQ=255)                                                                                      cTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGc                                                                                                            >  2:850922/1‑79 (MQ=255)                                                                                            aTTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTGAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAAc                                                                  >  2:1065473/1‑115 (MQ=255)                                                                                            aTTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAAc                                                                  <  1:1065473/115‑1 (MQ=255)                                                                                               ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGGAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                                       >  2:1068943/1‑139 (MQ=255)                                                                                               ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                    <  2:407809/78‑1 (MQ=255)                                                                                               ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                    >  1:407809/1‑78 (MQ=255)                                                                                               ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                                       <  1:1205173/139‑1 (MQ=255)                                                                                               ggCAATTACCGGGCC‑GGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTg                                      >  1:159161/1‑139 (MQ=255)                                                                                                    ttACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGc                                  >  1:286768/1‑139 (MQ=255)                                                                                                      aCCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGTATGCGGCGTGAACGCCTTATCCGGCATCTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTa                                <  1:899385/139‑1 (MQ=255)                                                                                                       ccgggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg                                   <  2:901952/135‑1 (MQ=255)                                                                                                       ccgggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg                                   >  1:901952/1‑135 (MQ=255)                                                                                                          ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                                >  1:155850/1‑55 (MQ=255)                                                                                                          ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                                >  1:1085255/1‑55 (MQ=255)                                                                                                          ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                                <  2:155850/55‑1 (MQ=255)                                                                                                          ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                                <  2:1085255/55‑1 (MQ=255)                                                                                                          ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                 <  1:1412041/86‑1 (MQ=255)                                                                                                          ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                 >  2:1412041/1‑86 (MQ=255)                                                                                                               ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                           <  1:562666/55‑1 (MQ=255)                                                                                                               ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                           >  2:562666/1‑55 (MQ=255)                                                                                                                   cAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATc                   <  1:1068943/139‑1 (MQ=255)                                                                                                                      tttCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                                <  1:580963/43‑1 (MQ=255)                                                                                                                      tttCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                                >  2:580963/1‑43 (MQ=255)                                                                                                                      tttCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCAcc                <  2:1303310/139‑1 (MQ=255)                                                                                                                           ctcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                                >  1:618692/1‑86 (MQ=255)                                                                                                                           ctcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                                <  2:618692/86‑1 (MQ=255)                                                                                                                            tcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                      >  1:486301/1‑47 (MQ=255)                                                                                                                            tcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                      <  2:486301/47‑1 (MQ=255)                                                                                                                                  ggCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCg    <  1:1362802/139‑1 (MQ=255)                                                                                                                                    cGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgat  >  2:341676/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                      GTCACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGAT  >  NC_000913/2673632‑2673901

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