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_000913967,1180AG100.0% 147.1 / NA 41Q166Q (CAA→CAGmsbAlipid ABC transporter permease/ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/18);  total (23/18)

GGTACGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTG  >  NC_000913/966984‑967253
                                                                                                                                      |                                                                                                                                       
ggTACGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGt                                                                                                                                     <  1:697761/139‑1 (MQ=255)
 gTACGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTc                                                                                                                                    >  1:577448/1‑139 (MQ=255)
            tCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTTAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGa                                                                                                                         <  2:448652/139‑1 (MQ=255)
            tCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGACCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGa                                                                                                                         >  1:222574/1‑139 (MQ=255)
                 tATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATtgtg                                                                                                                    >  2:706884/1‑139 (MQ=255)
                        tACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCAc                                                                                                             >  2:11435/1‑139 (MQ=255)
                                           cttcttcttcttCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGAt                                                                                          <  2:577448/139‑1 (MQ=255)
                                                tcttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcgcg                                                                                     <  1:706884/139‑1 (MQ=255)
                                                 cttcttcCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGt                                                                                    >  2:443659/1‑139 (MQ=255)
                                                                          tgtgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATc                                                           <  1:11435/139‑1 (MQ=255)
                                                                          tgtgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATc                                                           >  1:262263/1‑139 (MQ=255)
                                                                            tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcg                                                                                       >  2:596912/1‑109 (MQ=255)
                                                                            tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcg                                                                                       <  1:596912/109‑1 (MQ=255)
                                                                            tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAg                                                         >  1:2555/1‑139 (MQ=255)
                                                                            tgCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAg                                                         >  2:484994/1‑139 (MQ=255)
                                                                              cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTa                                                       >  2:614472/1‑139 (MQ=255)
                                                                              cGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTa                                                       >  2:381890/1‑139 (MQ=255)
                                                                                 gAAGGTGCGTCGATCATCGACCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTaaaa                                                    <  2:222574/139‑1 (MQ=255)
                                                                                      tGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATg                                               <  2:470474/139‑1 (MQ=255)
                                                                                        cGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCATCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCa                                             <  2:521518/139‑1 (MQ=255)
                                                                                           cGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGaa                                          <  1:381890/139‑1 (MQ=255)
                                                                                               catcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGAt                                                                                                                        >  1:178330/1‑57 (MQ=255)
                                                                                               catcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGAt                                                                                                                        <  2:178330/57‑1 (MQ=255)
                                                                                               catcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACAcc                                      >  2:114071/1‑139 (MQ=255)
                                                                                                atcGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCa                                     <  2:37804/139‑1 (MQ=255)
                                                                                                       tGTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGCAAAAACATGCAGAACACCATGGGGCa                              >  2:169341/1‑139 (MQ=255)
                                                                                                        gTTCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAg                             >  1:665463/1‑139 (MQ=255)
                                                                                                         ttCATCATGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAgg                            >  2:262633/1‑139 (MQ=255)
                                                                                                           catcatGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTg                          <  1:262633/139‑1 (MQ=255)
                                                                                                              catGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcg                                                                                       <  1:153748/75‑1 (MQ=255)
                                                                                                              catGATGTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcg                                                                                       >  2:153748/1‑75 (MQ=255)
                                                                                                                   tgTTCTATTACAGTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAg                  <  1:639511/139‑1 (MQ=255)
                                                                                                                               gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcg                                                                                       >  2:497821/1‑58 (MQ=255)
                                                                                                                               gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTcg                                                                                       <  1:497821/58‑1 (MQ=255)
                                                                                                                               gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                                              >  2:414128/1‑83 (MQ=255)
                                                                                                                               gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAAc                                                              <  1:414128/83‑1 (MQ=255)
                                                                                                                               gTTGGCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAAt      >  1:676767/1‑139 (MQ=255)
                                                                                                                                  ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGtt                                                                                   <  1:528713/59‑1 (MQ=255)
                                                                                                                                  ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGtt                                                                                   >  2:528713/1‑59 (MQ=255)
                                                                                                                                  ggCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCt   >  1:11147/1‑139 (MQ=255)
                                                                                                                                   gCAGCTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTg  >  1:732646/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GGTACGCTGTTGTCACGTATTACCTACGATTCCGAACAGGTTGCTTCTTCTTCTTCCGGCGCACTGATTACTGTTGTGCGTGAAGGTGCGTCGATCATCGGCCTGTTCATCATGATGTTCTATTACAGTTGGCAACTGTCGATCATTTTGATTGTGCTGGCACCGATTGTTTCGATTGCGATTCGCGTTGTATCGAAGCGTTTTCGCAACATCAGTAAAAACATGCAGAACACCATGGGGCAGGTGACCACCAGCGCAGAACAAATGCTG  >  NC_000913/966984‑967253

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