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

CGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATG  >  NC_000913/966988‑967250                                                                                                                                   |                                                                                                                                     cGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                              >  1:1679834/1‑139 (MQ=255)        gtCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGATGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTg                                                                                                                       >  2:1386850/1‑139 (MQ=255)          cACGTATTACCTACGATTCCGCACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGTCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGAt                                                                                                                     <  2:325493/139‑1 (MQ=255)                ttACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                               <  1:655748/139‑1 (MQ=255)                  aCCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTgg                                                                                                             >  2:445364/1‑139 (MQ=255)                             cGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGtt                                                                                                  >  1:1303114/1‑139 (MQ=255)                                    gTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTg                                                                                           >  1:1044768/1‑139 (MQ=255)                                    gTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTg                                                                                           >  1:29295/1‑139 (MQ=255)                                             tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                  <  2:1679834/139‑1 (MQ=255)                                             tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTATCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                  <  2:1303114/139‑1 (MQ=255)                                              cttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGt                                                                                 >  1:1501149/1‑139 (MQ=255)                                                 cttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGt                                                                              >  2:34180/1‑139 (MQ=255)                                                       gcgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                              <  2:370459/85‑1 (MQ=255)                                                       gcgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                              >  1:1301988/1‑85 (MQ=255)                                                       gcgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                              <  2:1301988/85‑1 (MQ=255)                                                       gcgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                              >  1:370459/1‑85 (MQ=255)                                                              gATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                  <  2:51576/122‑1 (MQ=255)                                                              gATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                  >  1:51576/1‑122 (MQ=255)                                                                         tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAg                                                      >  1:331377/1‑139 (MQ=255)                                                                           cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAtt                                                                                            >  2:902785/1‑99 (MQ=255)                                                                           cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAtt                                                                                            <  1:902785/99‑1 (MQ=255)                                                                                aGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAtt                                                                                            >  2:807152/1‑94 (MQ=255)                                                                                aGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAtt                                                                                            <  1:807152/94‑1 (MQ=255)                                                                                   tGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATg                                            <  1:1235908/139‑1 (MQ=255)                                                                                     cGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGa                                                                         <  1:1347100/108‑1 (MQ=255)                                                                                     cGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCa                                          >  1:620311/1‑139 (MQ=255)                                                                                     cGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCa                                          <  2:1044768/139‑1 (MQ=255)                                                                                     cGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATGGTTTCGATTGCGATTCGCGTTGTATCGa                                                                         >  2:1347100/1‑108 (MQ=255)                                                                                        cGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                       <  2:1501149/139‑1 (MQ=255)                                                                                        cGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                       <  1:34180/139‑1 (MQ=255)                                                                                          atcatcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAaca                                     >  2:393258/1‑139 (MQ=255)                                                                                             atcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCa                                  <  2:386580/139‑1 (MQ=255)                                                                                             atcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCa                                  <  2:780467/139‑1 (MQ=255)                                                                                                ggCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATgg                               <  1:1386850/139‑1 (MQ=255)                                                                                                    tGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCa                           >  1:574569/1‑139 (MQ=255)                                                                                                        catcatGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTg                       <  2:331377/139‑1 (MQ=255)                                                                                                            atgatgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATc                                                                           <  2:411756/83‑1 (MQ=255)                                                                                                            atgatgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATc                                                                           >  1:411756/1‑83 (MQ=255)                                                                                                               atgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAt                                                                                             <  1:1644772/62‑1 (MQ=255)                                                                                                               atgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAt                                                                                             >  2:1644772/1‑62 (MQ=255)                                                                                                                tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACa                                                          <  1:209385/96‑1 (MQ=255)                                                                                                                tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACa                                                          >  2:209385/1‑96 (MQ=255)                                                                                                                tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAg               <  1:393258/139‑1 (MQ=255)                                                                                                                tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAg               >  1:429999/1‑139 (MQ=255)                                                                                                                  ttCTATTACAGTTGGCAGCTGTCGATCATTTTGATtgtg                                                                                                                 >  2:163388/1‑39 (MQ=255)                                                                                                                  ttCTATTACAGTTGGCAGCTGTCGATCATTTTGATtgtg                                                                                                                 <  1:163388/39‑1 (MQ=255)                                                                                                                            gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGa                                                                         >  2:1589967/1‑69 (MQ=255)                                                                                                                            gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGa                                                                         <  1:1589967/69‑1 (MQ=255)                                                                                                                            gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                       <  2:691354/103‑1 (MQ=255)                                                                                                                            gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                       >  1:691354/1‑103 (MQ=255)                                                                                                                            gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAg               >  1:1424054/1‑127 (MQ=255)                                                                                                                            gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAg               <  2:1424054/127‑1 (MQ=255)                                                                                                                             ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                                       <  2:614571/54‑1 (MQ=255)                                                                                                                             ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                                       >  1:614571/1‑54 (MQ=255)                                                                                                                             ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATg  >  2:1350017/1‑139 (MQ=255)                                                                                                                             ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATg  >  2:1001870/1‑139 (MQ=255)                                                                                                                               ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGtt                                                                   >  2:1000214/1‑72 (MQ=255)                                                                                                                               ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGtt                                                                   <  1:1000214/72‑1 (MQ=255)                                                                                                                                   |                                                                                                                                     CGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATG  >  NC_000913/966988‑967250

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