Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 104,664 A→G H228R (CAC→CGC)  ftsA → ATP‑binding cell division FtsK recruitment protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913104,6640AG100.0% 150.7 / NA 42H228R (CAC→CGC) ftsAATP‑binding cell division FtsK recruitment protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/29);  total (13/29)

GAAAGTTGACCAACTGATATTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCACACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGA  >  NC_000913/104521‑104799
                                                                                                                                               |                                                                                                                                       
gAAAGTTGACCAACTGATATTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTaa                                                                                                                                   <  2:350058/150‑1 (MQ=255)
  aaGTTGACCAACTGATATTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGt                                                                                                                                >  1:276667/1‑151 (MQ=255)
  aaGTTGACCAACTGATATTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGt                                                                                                                                <  2:508409/151‑1 (MQ=255)
       gACCAACTGATATTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAAtt                                                                                                                            <  1:655389/150‑1 (MQ=255)
           aaCTGATATTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTa                                                                                                                       <  2:216975/151‑1 (MQ=255)
           aaCTGATATTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTa                                                                                                                       <  1:39776/151‑1 (MQ=255)
                atatTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTg                                                                                                                  <  1:225668/151‑1 (MQ=255)
                atatTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTg                                                                                                                  <  2:399522/151‑1 (MQ=255)
                 tatTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTg                                                                                                                  <  2:339708/150‑1 (MQ=255)
                     tGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCaa                                                                                                              <  1:813721/150‑1 (MQ=255)
                       ccGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGt                                                                                                           >  2:181657/1‑151 (MQ=255)
                        cGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGt                                                                                                           >  1:184690/1‑150 (MQ=255)
                        cGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTc                                                                                                          >  2:580321/1‑151 (MQ=255)
                           aCTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTg                                                                                                       <  2:412524/151‑1 (MQ=255)
                               gCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCa                                                                                                   <  1:386189/151‑1 (MQ=255)
                                    aaGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGat                                                                                              <  2:349915/151‑1 (MQ=255)
                                    aaGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGat                                                                                              <  2:839562/151‑1 (MQ=255)
                                      gTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGatat                                                                                            <  1:489501/151‑1 (MQ=255)
                                        tATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATc                                                                                           <  2:415988/150‑1 (MQ=255)
                                            cGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTa                                                                                      <  2:699727/151‑1 (MQ=255)
                                             ggTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTAc                                                                                     <  1:802256/151‑1 (MQ=255)
                                                aTTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGcc                                                                                  <  2:289283/151‑1 (MQ=255)
                                                     cGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTgg                                                                             <  2:248624/151‑1 (MQ=255)
                                                        aaGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCAc                                                                          >  1:50207/1‑151 (MQ=255)
                                                             gAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCAcgccg                                                                      <  1:816353/150‑1 (MQ=255)
                                                                      cTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACg                                                            >  2:256604/1‑151 (MQ=255)
                                                                         ggTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCg                                                         <  1:154041/151‑1 (MQ=255)
                                                                            gtCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAg                                                      <  1:201478/151‑1 (MQ=255)
                                                                               tGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCg                                                    >  1:627863/1‑150 (MQ=255)
                                                                               tGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCg                                                    >  2:796359/1‑150 (MQ=255)
                                                                                     gtcgATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTaa                                               >  2:355640/1‑149 (MQ=255)
                                                                                                 ggtggtACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGt                                  <  2:333663/150‑1 (MQ=255)
                                                                                                    ggtACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTtgt                               <  1:405691/150‑1 (MQ=255)
                                                                                                      tACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTgcg                            >  2:168078/1‑151 (MQ=255)
                                                                                                          aTGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTgg                        >  2:91536/1‑151 (MQ=255)
                                                                                                             gATATCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGtt                     <  2:145129/151‑1 (MQ=255)
                                                                                                              atatCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTc                    <  1:444750/151‑1 (MQ=255)
                                                                                                                atCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGcc                                                          <  2:119371/111‑1 (MQ=255)
                                                                                                                atCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGcc                                                          >  1:119348/1‑111 (MQ=255)
                                                                                                                 tCGCCGTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCAt                 >  1:564903/1‑151 (MQ=255)
                                                                                                                      gTTTATACCGGTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTg            <  2:349474/151‑1 (MQ=255)
                                                                                                                                gTGGGGCATTGCGCCGCACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATga  <  2:658591/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                       
GAAAGTTGACCAACTGATATTTGCCGGACTGGCATCAAGTTATTCGGTATTGACGGAAGATGAACGTGAACTGGGTGTCTGCGTCGTCGATATCGGTGGTGGTACAATGGATATCGCCGTTTATACCGGTGGGGCATTGCGCCACACTAAGGTAATTCCTTATGCTGGCAATGTCGTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGA  >  NC_000913/104521‑104799

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