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

TTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAG  >  NC_000913/966993‑967237                                                                                                                              |                                                                                                                        ttgtCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCAttt                                                                                                            >  2:143116/1‑139 (MQ=255)           ttACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                  >  1:677852/1‑139 (MQ=255)                                    ttcttcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAtt                                                                               <  1:254833/133‑1 (MQ=255)                                    ttcttcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAtt                                                                               >  2:254833/1‑133 (MQ=255)                                        tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                     <  2:677852/139‑1 (MQ=255)                                               ccGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATc                                                              >  2:57409/1‑139 (MQ=255)                                                cGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                 <  1:1185203/87‑1 (MQ=255)                                                cGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                 >  2:1185203/1‑87 (MQ=255)                                                 ggCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGa                                                            <  2:28553/139‑1 (MQ=255)                                                  gcgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                     >  1:1368754/1‑129 (MQ=255)                                                  gcgcACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                     <  2:1368754/129‑1 (MQ=255)                                                          aTTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTc                                                   >  1:480627/1‑139 (MQ=255)                                                            tACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGc                                                 >  2:262064/1‑139 (MQ=255)                                                                gttgtGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACa                                             <  1:262064/139‑1 (MQ=255)                                                                    tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                 <  2:475229/67‑1 (MQ=255)                                                                    tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGAt                                                                                                                 >  1:475229/1‑67 (MQ=255)                                                                      cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATtgtg                                                                                                    <  2:321230/78‑1 (MQ=255)                                                                      cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATtgtg                                                                                                    >  1:321230/1‑78 (MQ=255)                                                                      cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCg                                                  >  2:610006/1‑128 (MQ=255)                                                                      cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCg                                                  <  1:610006/128‑1 (MQ=255)                                                                      cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACAt                                            <  1:1116599/134‑1 (MQ=255)                                                                      cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACAt                                            >  2:1116599/1‑134 (MQ=255)                                                                           aGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                  >  2:1382348/1‑139 (MQ=255)                                                                              tGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATg                               <  1:143116/139‑1 (MQ=255)                                                                                   cGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                          <  1:57409/139‑1 (MQ=255)                                                                                          cGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCAc                                                                                             >  1:1412641/1‑65 (MQ=255)                                                                                          cGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCAc                                                                                             <  2:1412641/65‑1 (MQ=255)                                                                                          cGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                  >  2:1152618/1‑124 (MQ=255)                                                                                          cGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAAc                                  <  1:1152618/124‑1 (MQ=255)                                                                                                gTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                              >  2:666431/1‑106 (MQ=255)                                                                                                gTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                              <  1:666431/106‑1 (MQ=255)                                                                                                gTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAg             <  1:1382348/139‑1 (MQ=255)                                                                                                 ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCAc                                                                                             <  2:394658/58‑1 (MQ=255)                                                                                                 ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCAc                                                                                             >  1:394658/1‑58 (MQ=255)                                                                                                 ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                     >  2:462840/1‑82 (MQ=255)                                                                                                 ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                     <  1:462840/82‑1 (MQ=255)                                                                                                  tcatcaTGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGt           >  2:209650/1‑139 (MQ=255)                                                                                                      catGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                              >  1:912401/1‑100 (MQ=255)                                                                                                      catGATGTTCTATTACAGTTAGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                              <  2:912401/100‑1 (MQ=255)                                                                                                        tgatgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGaccac     <  2:480627/139‑1 (MQ=255)                                                                                                           tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGAt                                                                                                        >  1:842671/1‑37 (MQ=255)                                                                                                           tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGAt                                                                                                        <  2:842671/37‑1 (MQ=255)                                                                                                           tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                  >  1:901256/1‑43 (MQ=255)                                                                                                           tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCt                                                                                                  <  2:901256/43‑1 (MQ=255)                                                                                                           tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAg  >  2:1190188/1‑139 (MQ=255)                                                                                                                       gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAtt                                                                               <  2:1448652/50‑1 (MQ=255)                                                                                                                       gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGAtt                                                                               >  1:1448652/1‑50 (MQ=255)                                                                                                                       gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                          >  2:220900/1‑103 (MQ=255)                                                                                                                       gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                          <  1:220900/103‑1 (MQ=255)                                                                                                                        ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGtt                                                                   <  1:298158/61‑1 (MQ=255)                                                                                                                        ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGtt                                                                   >  2:298158/1‑61 (MQ=255)                                                                                                                        ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCa                             >  2:1127955/1‑99 (MQ=255)                                                                                                                        ttGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCa                             <  1:1127955/99‑1 (MQ=255)                                                                                                                          ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCg                                                                                           <  2:1487412/35‑1 (MQ=38)                                                                                                                          ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCg                                                                                           >  1:1487412/1‑35 (MQ=38)                                                                                                                          ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                          <  1:763558/52‑1 (MQ=255)                                                                                                                          ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                          >  2:763558/1‑52 (MQ=255)                                                                                                                          ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCa                                                >  2:282387/1‑78 (MQ=255)                                                                                                                          ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCa                                                <  1:282387/78‑1 (MQ=255)                                                                                                                          ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAacac                       >  2:524200/1‑103 (MQ=255)                                                                                                                          ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAacac                       <  1:524200/103‑1 (MQ=255)                                                                                                                              |                                                                                                                        TTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAG  >  NC_000913/966993‑967237

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