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

TCACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTT  >  NC_000913/2673633‑2673892                                                                                                                                        |                                                                                                                             tCACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGaa                                                                                                                           >  1:469734/1‑139 (MQ=255)     cATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGAt                                                                                                                       >  2:119188/1‑139 (MQ=255)     cATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGAt                                                                                                                       >  2:571375/1‑139 (MQ=255)               atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcg                                                                                                             >  2:197134/1‑139 (MQ=255)               atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcg                                                                                                             >  1:1573156/1‑139 (MQ=255)                tatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcgg                                                                                                            <  1:14299/139‑1 (MQ=255)                  tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCCATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATgcggcg                                                                                                          >  2:161091/1‑139 (MQ=255)                    aCCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTg                                                                                                        <  1:571375/139‑1 (MQ=255)                    aCCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTg                                                                                                        <  1:197134/139‑1 (MQ=255)                    aCCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTg                                                                                                        <  1:119188/139‑1 (MQ=255)                     ccGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                       <  1:111334/139‑1 (MQ=255)                     ccGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                       <  2:1573156/139‑1 (MQ=255)                                   gggCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                         <  1:1479537/139‑1 (MQ=255)                                         aaTATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                                   <  1:985063/139‑1 (MQ=255)                                          atatGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                       <  1:243165/118‑1 (MQ=255)                                          atatGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                       >  2:243165/1‑118 (MQ=255)                                                  aaTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACa                                                                          <  2:469734/139‑1 (MQ=255)                                                                        tCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAt                                                    >  2:802161/1‑139 (MQ=255)                                                                         cGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                       <  1:893652/135‑1 (MQ=255)                                                                         cGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                       >  2:893652/1‑135 (MQ=255)                                                                           cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                                 >  2:255596/1‑139 (MQ=255)                                                                              gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                                   <  2:854087/102‑1 (MQ=255)                                                                              gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                                   >  1:854087/1‑102 (MQ=255)                                                                              gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                                     <  1:1181616/116‑1 (MQ=255)                                                                              gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                                     >  2:1181616/1‑116 (MQ=255)                                                                              gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGa                                              >  2:1371722/1‑139 (MQ=255)                                                                                   ttCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                                         >  1:977476/1‑139 (MQ=255)                                                                                         cTATTGGCAATTCCCGGTCCGTGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAGCGCCTTATCCGGCATTTAATAAATTACAGCCACTCTTCCCGCAACGACATCTCTGAAGTTGTTGAtt                                   <  1:161091/139‑1 (MQ=255)                                                                                             tGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAGCGCCTTATCCGGCATTtaat                                                                                  >  2:1005102/1‑88 (MQ=255)                                                                                             tGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaac                                                                                  <  1:1005102/88‑2 (MQ=255)                                                                                             tGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTcccc                               <  1:1371722/139‑1 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGTCGGATGCGGCGTg                                                                                                        >  1:369755/1‑65 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGTCGGATGCGGCGTg                                                                                                        <  2:369755/65‑1 (MQ=38)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                                  <  1:1201781/87‑1 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                                  >  2:1201781/1‑87 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                                   >  1:169070/1‑86 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                                   <  2:169070/86‑1 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                       >  1:143839/1‑114 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTg                                               >  1:1114723/1‑122 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTg                                               <  2:1114723/122‑1 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                              >  1:1091694/1‑139 (MQ=255)                                                                                              ggCAATTACCGGGCCGGGCGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                       <  2:143839/114‑1 (MQ=255)                                                                                                  aTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg                          >  2:472037/1‑139 (MQ=255)                                                                                                   ttACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGc                         >  2:1057192/1‑139 (MQ=255)                                                                                                      ccgggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg                      <  2:977476/139‑1 (MQ=255)                                                                                                      ccgggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg                      <  1:255596/139‑1 (MQ=255)                                                                                                         ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                       >  2:1264094/1‑55 (MQ=255)                                                                                                         ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                       <  1:1264094/55‑1 (MQ=255)                                                                                                         ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                             >  1:1524537/1‑65 (MQ=255)                                                                                                         ggccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                             <  2:1524537/65‑1 (MQ=255)                                                                                                          gccgggGGCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCa                  <  2:78415/139‑1 (MQ=255)                                                                                                              ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATg                                                                                                               <  2:971369/42‑1 (MQ=255)                                                                                                              ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATg                                                                                                               >  1:971369/1‑42 (MQ=255)                                                                                                              ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                             <  2:1080103/60‑1 (MQ=255)                                                                                                              ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                             >  1:1080103/1‑60 (MQ=255)                                                                                                              ggggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATga              >  1:1190942/1‑139 (MQ=255)                                                                                                               gggCAATTTCTCTCGATCGGCGTTCAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                        >  1:1091593/1‑64 (MQ=255)                                                                                                               gggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                        <  2:1091593/64‑1 (MQ=255)                                                                                                                ggCAATTTCTCTCGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATgaga            <  2:1091694/139‑1 (MQ=255)                                                                                                                          ctcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACa                                                                          >  1:1388739/1‑67 (MQ=255)                                                                                                                          ctcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACa                                                                          <  2:1388739/67‑1 (MQ=255)                                                                                                                          ctcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                       >  1:960210/1‑86 (MQ=255)                                                                                                                          ctcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGa                                                       <  2:960210/86‑1 (MQ=255)                                                                                                                          ctcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCGCCCCCtt  <  1:1057192/139‑1 (MQ=255)                                                                                                                           tcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                       >  1:107126/1‑37 (MQ=25)                                                                                                                           tcGATCGGCGTTAAGAAAGATGCCGGATGCGGCGTGa                                                                                                       <  2:107126/37‑1 (MQ=25)                                                                                                                                   cGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAccc                                                              <  2:1130105/70‑1 (MQ=255)                                                                                                                                   cGTTAAGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAccc                                                              >  1:1130105/1‑70 (MQ=255)                                                                                                                                        |                                                                                                                             TCACCATTACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTT  >  NC_000913/2673633‑2673892

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