Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,676,572 G→A P85L (CCG→CTG)  yphD ← putative sugar ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,676,5720GA100.0% 186.9 / NA 56P85L (CCG→CTG) yphDputative sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/19);  total (37/19)

AGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCA  >  NC_000913/2676476‑2676665
                                                                                                |                                                                                             
aGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCTGCCAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGcc                                                                                            >  1:1279285/1‑100 (MQ=255)
   gcgcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCaa                                                                                         >  2:379756/1‑100 (MQ=255)
     gcCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACg                                                                                       >  1:69060/1‑100 (MQ=255)
      ccTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAACATGCCAGGCTCACCGAGACAAAAGCCACCATCAGcacac                                                                                         >  1:1428043/1‑96 (MQ=255)
        tAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTg                                                                                    <  1:424758/100‑1 (MQ=255)
        tAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGACACACAGAGACAAAAGCAACCATCAGCCCAACGCTg                                                                                    >  2:572312/1‑100 (MQ=255)
        tAACAGCAACACCAGCAGGCACGCCACCGCCAGCAGAACTTCAAATTGCAGCAAACATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCa                                                                                          >  1:1257547/1‑94 (MQ=255)
          aCAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGAc                                                                                  <  1:1318881/100‑1 (MQ=255)
            aGCAACACCCGCAGGCACGCAACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACa                                                                                 >  1:1779376/1‑99 (MQ=255)
             ncAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGACACAAAAGCCACCATCAGCCCAACGCTGACATc                                                                               >  2:379101/2‑100 (MQ=255)
              cAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCa                                                                              >  2:759772/1‑100 (MQ=255)
                 cacCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAAtt                                                                           <  1:379754/100‑1 (MQ=255)
                 cacCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAAtt                                                                           <  1:323817/100‑1 (MQ=255)
                     agcagGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCAc                                                                       >  1:651955/1‑100 (MQ=255)
                     agcagGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGAACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCAc                                                                       <  2:12499/100‑1 (MQ=255)
                        agGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTacac                                                                       >  1:1657610/1‑97 (MQ=255)
                        agGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCgg                                                                    <  2:856748/100‑1 (MQ=255)
                         gccACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGTCCAGCCCGGACAGCAATTTCACCGga                                                                   >  1:458661/3‑100 (MQ=255)
                            aCGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGat                                                                >  2:559632/1‑100 (MQ=255)
                             cGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGCCATCAATTTCACCGgagaca                                                               >  1:1017348/1‑98 (MQ=255)
                             cGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGata                                                               <  1:1482305/100‑1 (MQ=255)
                              gCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGAAGataa                                                              >  1:1497498/1‑100 (MQ=255)
                               ccACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGataat                                                             >  2:1641574/1‑100 (MQ=255)
                               ccACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATAAGACCAACGCTGACATCAATTTCACCGGAGataat                                                             >  2:953182/1‑100 (MQ=255)
                                cACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGataata                                                            <  1:709421/100‑1 (MQ=255)
                                 aCCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGAtaataa                                                           >  2:115954/1‑100 (MQ=255)
                                      cAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGc                                                      >  2:1550794/1‑100 (MQ=255)
                                       aGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCg                                                     >  2:855549/1‑100 (MQ=255)
                                          ggAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCa                                                  >  2:1545302/1‑100 (MQ=255)
                                             aCTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCg                                               <  1:839478/100‑1 (MQ=255)
                                                tCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGccc                                            <  2:730908/100‑1 (MQ=255)
                                                tCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGccc                                            <  2:700325/100‑1 (MQ=255)
                                                       gcagcaAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAggcgg                                      <  2:1279296/99‑1 (MQ=255)
                                                        cagcaAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCa                                    >  1:1045354/1‑100 (MQ=255)
                                                        cagcaAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAAGTTCACCAGAGAAAATAATCAGCGTCAGCGCCCAgccggc                                     >  1:618521/1‑99 (MQ=255)
                                                         agcaAAAATGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCaa                                   >  2:913163/1‑100 (MQ=255)
                                                                 tGCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAAtg                           <  1:559628/100‑1 (MQ=255)
                                                                  gCCAGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAAtgg                          >  1:1579966/1‑100 (MQ=255)
                                                                   cgaGGCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAAtggt                         >  2:1228191/3‑100 (MQ=255)
                                                                      ggCACACCGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTggc                      <  2:490352/100‑1 (MQ=255)
                                                                      ggCACACCGAGACAAAAGCCACCAACAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAACGTCAACGCACCGGCGACAATCCCaaacg                          >  1:1445859/1‑93 (MQ=255)
                                                                        cacaCCGAGACAAAAGCCACCATCAGCCCAGCGCTGACAGCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTggcgn                    <  1:379099/100‑2 (MQ=255)
                                                                             ctataCAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATcg               <  2:1710785/96‑1 (MQ=255)
                                                                             cGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATcg               <  2:1779389/100‑1 (MQ=255)
                                                                             cGAGACAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATcg               >  1:1680400/1‑100 (MQ=255)
                                                                              gagaCAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATcgc              >  1:1744362/1‑100 (MQ=255)
                                                                              gagaCAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATcgc              >  1:1739136/1‑100 (MQ=255)
                                                                              gagaCAAAAGACACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATcgc              >  1:1739122/1‑100 (MQ=255)
                                                                                gaCAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCgcgc            >  2:1017445/1‑100 (MQ=255)
                                                                                gaCAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCgcgc            >  2:1412142/1‑100 (MQ=255)
                                                                                 aCAAAAGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgca           <  1:452235/100‑1 (MQ=255)
                                                                                    aaaGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAACGGTGGCGGCATCGCgcagca        <  1:84232/100‑1 (MQ=255)
                                                                                     aaGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCgcagc         >  1:837854/1‑98 (MQ=255)
                                                                                      aGCCACCATCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGAAGCGGCATAGCGCAGCACg      >  2:453754/1‑100 (MQ=255)
                                                                                        ccaccaTCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGtt    >  2:1100233/1‑100 (MQ=255)
                                                                                          accaTCAGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCa  >  2:1579134/1‑100 (MQ=255)
                                                                                                |                                                                                             
AGCGCGCCTAACAGCAACACCAGCAGGCACGCCACCGCCAGCGGAACTTCAAATTGCAGCAAAAATGCCAGGCACACCGAGACAAAAGCCACCATCGGCCCAACGCTGACATCAATTTCACCGGAGATAATAATCAGCGTCATCGCCCAGGCGGCAATCCCAATGGTGGCGGCATCGCGCAGCACGTTCA  >  NC_000913/2676476‑2676665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: