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%	138.3 / NA	43	*141* (TAG→TAA)	yphA	DoxX family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/18);  total (25/18)

TGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGA  >  NC_000913/2673643‑2673900                                                                                                                              |                                                                                                                                     tGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGa                                                                                                                         >  2:796993/1‑139 (MQ=255)     atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcg                                                                                                                     >  1:212404/1‑139 (MQ=255)        tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcggcg                                                                                                                  >  2:802276/1‑139 (MQ=255)        tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcggcg                                                                                                                  >  2:22069/1‑139 (MQ=255)          aCCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTg                                                                                                                <  1:796993/139‑1 (MQ=255)            cGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGaa                                                                                                              <  1:284878/139‑1 (MQ=255)               cGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGc                                                                                                           >  1:465993/1‑139 (MQ=255)                                       tAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAc                                                                                   <  2:37540/139‑1 (MQ=255)                                        aaTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACa                                                                                  <  1:558810/139‑1 (MQ=255)                                         aTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAg                                                                                 >  1:172122/1‑139 (MQ=255)                                          tttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                >  2:545759/1‑139 (MQ=255)                                                gaagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActct                                                                          <  2:212404/139‑1 (MQ=255)                                                gaagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActct                                                                          <  2:465993/139‑1 (MQ=255)                                                           gTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                               <  1:744385/139‑1 (MQ=255)                                                            tATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGAc                                                              <  2:172122/139‑1 (MQ=255)                                                               cGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTCCAGCCACTCTGCCCGCAACGACAtc                                                           <  1:586140/139‑1 (MQ=255)                                                                 cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                                         >  2:597585/1‑139 (MQ=255)                                                                         ttCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                                                 >  2:38702/1‑139 (MQ=255)                                                                           cTTGCTATTGTCACTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                                               >  1:536752/1‑139 (MQ=255)                                                                              gCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGAt                                            >  2:386665/1‑139 (MQ=255)                                                                                    ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                <  2:323138/97‑1 (MQ=255)                                                                                    ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                >  1:323138/1‑97 (MQ=255)                                                                                    ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                                      >  2:6261/1‑139 (MQ=255)                                                                                    ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                                      >  1:417151/1‑139 (MQ=255)                                                                                    ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                                      >  1:46543/1‑139 (MQ=255)                                                                                    ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                                      >  2:26302/1‑139 (MQ=255)                                                                                         ttACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGc                                 >  2:767356/1‑139 (MQ=255)                                                                                         ttACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGc                                 >  1:650047/1‑139 (MQ=255)                                                                                            ccgggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg                              <  2:46543/139‑1 (MQ=255)                                                                                              gggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                                         <  2:15266/110‑1 (MQ=255)                                                                                              gggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                                         >  1:15266/1‑110 (MQ=255)                                                                                               ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcgc                           >  1:752732/1‑139 (MQ=255)                                                                                                 ccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCaa                         <  1:802276/139‑1 (MQ=255)                                                                                                    ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                                             <  1:469939/116‑1 (MQ=255)                                                                                                    ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                                             >  2:469939/1‑116 (MQ=255)                                                                                                    ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATga                      >  2:757610/1‑139 (MQ=255)                                                                                                    ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATga                      >  2:764945/1‑139 (MQ=255)                                                                                                    ggggCAATTTCTCTCGATCGCCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAAtca                      >  1:471103/1‑137 (MQ=255)                                                                                                                  cGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTcc        <  1:386665/139‑1 (MQ=255)                                                                                                                    aTCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCgg      >  1:167243/1‑139 (MQ=255)                                                                                                                     tCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCggg     <  1:6261/139‑1 (MQ=255)                                                                                                                     tCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCggg     <  1:38702/139‑1 (MQ=255)                                                                                                                        gCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCga  <  1:26302/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                     TGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGA  >  NC_000913/2673643‑2673900

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: ‑