Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	967,118	A→G	Q166Q (CAA→CAG)	msbA →	lipid ABC transporter permease/ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	967,118	0	A	G	100.0%	239.4 / NA	67	Q166Q (CAA→CAG)	msbA	lipid ABC transporter permease/ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/35);  total (32/35)

TACGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCT  >  NC_000913/966986‑967252                                                                                                                                     |                                                                                                                                       tACGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCg                                                                                                                                  <  1:928554/139‑1 (MQ=255) tACGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCg                                                                                                                                  >  1:357963/1‑139 (MQ=255)   cGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                                <  2:357963/139‑1 (MQ=255)                  ttACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                                 >  2:277176/1‑139 (MQ=255)                  ttACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                                 >  1:769721/1‑139 (MQ=255)                  ttACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                                 >  1:127000/1‑139 (MQ=255)                  ttACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                                 >  2:470516/1‑139 (MQ=255)                     ccTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGc                                                                                                              >  1:419279/1‑139 (MQ=255)                          gATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATtgtg                                                                                                                   <  2:825607/129‑1 (MQ=255)                          gATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATtgtg                                                                                                                   >  1:825607/1‑129 (MQ=255)                               cGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGtt                                                                                                    >  2:385930/1‑139 (MQ=255)                               cGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTA‑TACAGTTGGCAGCTGTCg                                                                                                                                  <  1:1049588/108‑1 (MQ=255)                               cGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTA‑TACAGTTGGCAGCTGTCg                                                                                                                                  >  2:1049588/1‑108 (MQ=255)                                      gTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCAt                                                                                                                             >  2:154866/1‑107 (MQ=255)                                      gTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCAt                                                                                                                             <  1:154866/107‑1 (MQ=255)                                          cttcttcttcttCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                                         <  2:448638/139‑1 (MQ=255)                                               tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    <  1:757105/139‑1 (MQ=255)                                               tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    <  1:662292/139‑1 (MQ=255)                                               tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    <  1:470516/139‑1 (MQ=255)                                               tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    <  1:385930/139‑1 (MQ=255)                                               tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    <  2:379519/139‑1 (MQ=255)                                                 ttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGtt                                                                                  <  1:1009294/139‑1 (MQ=255)                                                          cgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAg                                                                         >  2:851023/1‑139 (MQ=255)                                                                       gttgtGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGGTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACa                                                            <  1:851023/139‑1 (MQ=255)                                                                           tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCg                                                                                                                                  >  2:1210838/1‑65 (MQ=255)                                                                           tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCg                                                                                                                                  <  1:1210838/65‑1 (MQ=255)                                                                           tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAg                                                        >  2:451657/1‑139 (MQ=255)                                                                           tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAg                                                        >  1:1210615/1‑139 (MQ=255)                                                                           tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAg                                                        >  1:846362/1‑139 (MQ=255)                                                                             cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTa                                                      <  1:451657/139‑1 (MQ=255)                                                                             cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTa                                                      >  2:865253/1‑139 (MQ=255)                                                                                  aGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                                 >  1:214238/1‑139 (MQ=255)                                                                                  aGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                                 <  1:277176/139‑1 (MQ=255)                                                                                       cGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCa                                            <  2:846362/139‑1 (MQ=255)                                                                                          cGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                         <  2:1210615/139‑1 (MQ=255)                                                                                            atcatcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    >  1:920006/1‑94 (MQ=255)                                                                                            atcatcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    <  2:920006/94‑1 (MQ=255)                                                                                               atcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCa                                    <  2:419279/139‑1 (MQ=255)                                                                                                   gCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATggg                                <  1:1155142/139‑1 (MQ=255)                                                                                                     cTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGc                              <  2:769721/139‑1 (MQ=255)                                                                                                     cTGTTCATCATGATGTTCTATTACAGTTGGCAGATGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGc                              <  2:214238/139‑1 (MQ=255)                                                                                                        ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                         >  2:460967/1‑125 (MQ=255)                                                                                                        ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                         <  1:460967/125‑1 (MQ=255)                                                                                                                  tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGa                                                                                                                        <  2:56918/36‑1 (MQ=255)                                                                                                                  tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGa                                                                                                                        >  1:56918/1‑36 (MQ=255)                                                                                                                  tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    >  2:111967/1‑72 (MQ=255)                                                                                                                  tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                    <  1:111967/72‑1 (MQ=255)                                                                                                                  tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                         <  2:253209/115‑1 (MQ=255)                                                                                                                  tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                         >  1:253209/1‑115 (MQ=255)                                                                                                                  tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGCCCACCAg                 >  1:229731/1‑139 (MQ=255)                                                                                                                  tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAg                 >  2:624250/1‑139 (MQ=255)                                                                                                                       tATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                                         >  2:571698/1‑62 (MQ=255)                                                                                                                       tATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                                         <  1:571698/62‑1 (MQ=255)                                                                                                                         ttACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGc                              <  1:1047444/119‑1 (MQ=255)                                                                                                                         ttACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGc                              >  2:1047444/1‑119 (MQ=255)                                                                                                                          tACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                                 >  2:1051832/1‑99 (MQ=255)                                                                                                                          tACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                                 <  1:1051832/99‑1 (MQ=255)                                                                                                                           aCAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACa        <  2:127000/139‑1 (MQ=255)                                                                                                                              gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                         >  2:430811/1‑103 (MQ=255)                                                                                                                              gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                         <  1:430811/103‑1 (MQ=255)                                                                                                                               ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAtt                                                                                        >  1:128222/1‑55 (MQ=255)                                                                                                                               ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAtt                                                                                        <  2:128222/55‑1 (MQ=255)                                                                                                                                 ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGt                                                                      <  1:1101431/71‑1 (MQ=255)                                                                                                                                 ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGt                                                                      >  2:1101431/1‑71 (MQ=255)                                                                                                                                 ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                                 <  2:1041715/92‑1 (MQ=255)                                                                                                                                 ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                                 >  1:1041715/1‑92 (MQ=255)                                                                                                                                 ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCt  >  1:506246/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                       TACGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCT  >  NC_000913/966986‑967252

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