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

GGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGG  >  NC_000913/2673645‑2673897                                                                                                                            |                                                                                                                                  ggATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAGGATGCCGGATg                                                                                                                    >  1:738960/1‑139 (MQ=255)      tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcggcg                                                                                                               >  1:173042/1‑139 (MQ=255)        aCCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTg                                                                                                             >  2:80372/1‑139 (MQ=255)              gATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGcc                                                                                                       <  2:738960/139‑1 (MQ=255)                       gggCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                              >  2:539287/1‑139 (MQ=255)                            aaaTATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTta                                                                                         <  1:539287/139‑1 (MQ=255)                              atatGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                                       <  1:80372/139‑1 (MQ=255)                                                          tATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGAc                                                           <  2:173042/139‑1 (MQ=255)                                                               cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                            >  1:1022929/1‑85 (MQ=255)                                                               cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                            <  2:1022929/85‑1 (MQ=255)                                                               cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                                      >  1:491892/1‑139 (MQ=255)                                                               cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                                      >  1:919808/1‑139 (MQ=255)                                                                  gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACg                                                                                                         >  2:718006/1‑85 (MQ=255)                                                                  gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACg                                                                                                         <  1:718006/85‑1 (MQ=255)                                                                     cgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgt                                                <  2:919808/139‑1 (MQ=255)                                                                      gTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgtt                                               >  1:586456/1‑139 (MQ=255)                                                                             cTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGAtt                                        <  2:491892/139‑1 (MQ=255)                                                                                ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTccc                                     <  2:586456/139‑1 (MQ=255)                                                                                 tGGCAATTACCGGGCCGGGGGGAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTcccc                                    <  2:647500/139‑1 (MQ=255)                                                                                  ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAt                                                                                   >  1:432899/1‑91 (MQ=255)                                                                                  ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAt                                                                                   <  2:432899/91‑1 (MQ=255)                                                                                  ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTa                                                                      <  1:812720/104‑1 (MQ=255)                                                                                  ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTa                                                                      >  2:812720/1‑104 (MQ=255)                                                                                  ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAt                                                         >  1:850860/1‑117 (MQ=255)                                                                                  ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAt                                                         <  2:850860/117‑1 (MQ=255)                                                                                          ccgggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAAc                                                                                                          <  2:31013/60‑1 (MQ=255)                                                                                          ccgggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAAc                                                                                                          >  1:31013/1‑60 (MQ=255)                                                                                             ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcgc                        >  1:195/1‑139 (MQ=255)                                                                                                  ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                                       >  2:482526/1‑71 (MQ=255)                                                                                                  ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                                       <  1:482526/71‑1 (MQ=255)                                                                                                  ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                                          >  2:103840/1‑116 (MQ=255)                                                                                                  ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                                          <  1:103840/116‑1 (MQ=255)                                                                                                            ctctcGATCGGCGTTAAGAAAGATGCCGGATGTGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg                               <  2:282049/117‑1 (MQ=255)                                                                                                            ctctcGATCGGCGTTAAGAAAGATGCCGGATGTGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg                               >  1:282049/1‑117 (MQ=255)                                                                                                                   tCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCggg  <  2:244225/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                  GGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGG  >  NC_000913/2673645‑2673897

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