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%	226.9 / NA	65	Q166Q (CAA→CAG)	msbA	lipid ABC transporter permease/ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/31);  total (34/31)

GTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTG  >  NC_000913/966995‑967253                                                                                                                            |                                                                                                                                         gtCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCg                                                                                                                                    >  1:785080/1‑130 (MQ=255) gtCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCg                                                                                                                                    <  2:785080/130‑1 (MQ=255) gtCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTg                                                                                                                           >  1:702953/1‑139 (MQ=255)         ttACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                                   >  2:1556110/1‑139 (MQ=255)         ttACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                                   >  1:536739/1‑139 (MQ=255)           aCCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                                  >  2:557620/1‑122 (MQ=255)           aCCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                                  <  1:557620/122‑1 (MQ=255)              tACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCAc                                                                                                              <  1:1294220/139‑1 (MQ=255)                cGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGTCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCg                                                                                                            >  1:173148/1‑139 (MQ=255)                                 cttcttcttcttCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                                           <  2:1588811/139‑1 (MQ=255)                                       cttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGt                                                                                     >  2:1109084/1‑139 (MQ=255)                                                gcgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGaa                                                                            >  1:908893/1‑139 (MQ=255)                                                gcgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGaa                                                                            <  2:536739/139‑1 (MQ=255)                                                   cACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCg                                                                         >  1:1090314/1‑139 (MQ=255)                                                     cTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCg                                                                                                                                    >  1:558942/1‑78 (MQ=255)                                                     cTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCg                                                                                                                                    <  2:558942/78‑1 (MQ=255)                                                             tgttgtGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCGGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                                               <  2:1090314/139‑1 (MQ=255)                                                               ttgtGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACAt                                                             >  1:147396/1‑139 (MQ=255)                                                                tgtgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATc                                                            >  1:1039668/1‑139 (MQ=255)                                                                  tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACa                                                              <  1:288154/135‑1 (MQ=255)                                                                  tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACa                                                              >  2:288154/1‑135 (MQ=255)                                                                    cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCg                                                                                                            >  1:1515604/1‑87 (MQ=255)                                                                    cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCg                                                                                                            <  2:1515604/87‑1 (MQ=255)                                                                      tGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCAttt                                                                                                                             >  2:784520/1‑68 (MQ=255)                                                                      tGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCAttt                                                                                                                             <  1:784520/68‑1 (MQ=255)                                                                      tGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑aaa                                                     >  2:1186734/1‑139 (MQ=255)                                                                         aGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAAc                                                  <  2:702953/139‑1 (MQ=255)                                                                                tCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACGTGCAGa                                           >  1:1434960/1‑139 (MQ=255)                                                                                 cGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGaa                                          <  2:894929/139‑1 (MQ=255)                                                                                 cGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGaa                                          <  1:1186734/139‑1 (MQ=255)                                                                                  gATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAACATGCAGaa                                          <  2:147396/139‑1 (MQ=255)                                                                                      atcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCa                                     <  2:1039668/139‑1 (MQ=255)                                                                                      atcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCa                                     <  2:908893/139‑1 (MQ=255)                                                                                               ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                                               <  2:271836/105‑1 (MQ=255)                                                                                               ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                                               >  1:271836/1‑105 (MQ=255)                                                                                                  atcatGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGc                               <  1:1278556/133‑1 (MQ=255)                                                                                                  atcatGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGc                               >  2:1278556/1‑133 (MQ=255)                                                                                                    catGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAg                             >  2:242148/1‑133 (MQ=255)                                                                                                    catGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAg                             <  1:242148/133‑1 (MQ=255)                                                                                                    catGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGacc                       >  2:991893/1‑139 (MQ=255)                                                                                                      tgatgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACGTGCAGAACACCATGGGGCAGGTGaccac                     <  2:1434960/139‑1 (MQ=255)                                                                                                           ttCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAgcg                >  1:85369/1‑139 (MQ=255)                                                                                                            tCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGc                                                                                                                >  1:783296/1‑43 (MQ=255)                                                                                                            tCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGc                                                                                                                <  2:783296/43‑1 (MQ=255)                                                                                                            tCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAgcgc               <  1:991893/139‑1 (MQ=255)                                                                                                              tATTACAGTTAGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGc                               <  2:1290707/121‑1 (MQ=255)                                                                                                              tATTACAGTTAGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGc                               >  1:1290707/1‑121 (MQ=255)                                                                                                                ttACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGaa           <  1:1109084/139‑1 (MQ=255)                                                                                                                ttACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGaa           >  1:501279/1‑139 (MQ=255)                                                                                                                 tACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAg                             <  1:890011/120‑1 (MQ=255)                                                                                                                 tACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAg                             >  2:890011/1‑120 (MQ=255)                                                                                                                     gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGtt                                                                                    >  1:892736/1‑62 (MQ=255)                                                                                                                     gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGtt                                                                                    <  2:892736/62‑1 (MQ=255)                                                                                                                     gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGaa                                          <  2:993366/103‑1 (MQ=255)                                                                                                                     gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTATTTCGATTGGGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGaa                                          >  1:993366/1‑103 (MQ=255)                                                                                                                      ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATg     >  2:854019/1‑139 (MQ=255)                                                                                                                        ggCAGCTGTCGATCATTTTGATTGTGTTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAg                                            >  1:139743/1‑98 (MQ=255)                                                                                                                        ggCAGCTGTCGATCATTTTGATTGTGTTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAg                                            <  2:139743/98‑1 (MQ=255)                                                                                                                        ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGaa                                          >  1:385769/1‑100 (MQ=255)                                                                                                                        ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGaa                                          <  2:385769/100‑1 (MQ=255)                                                                                                                        ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGaccac                     >  2:76748/1‑121 (MQ=255)                                                                                                                        ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGaccac                     <  1:76748/121‑1 (MQ=255)                                                                                                                        ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCt   >  1:335825/1‑139 (MQ=255)                                                                                                                        ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCt   >  1:439929/1‑139 (MQ=255)                                                                                                                         gCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTg  <  1:1556110/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                         GTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGT‑AAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTG  >  NC_000913/966995‑967253

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