Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,018,760 G→A *476* (TAG→TAA)  rmuC → DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,018,7600GA100.0% 223.8 / NA 69*476* (TAG→TAA) rmuCDNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/37);  total (32/37)

ACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTGGTTTT  >  NC_000913/4018628‑4018893
                                                                                                                                     |                                                                                                                                     
aCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                                 <  1:289614/139‑1 (MQ=255)
   cgAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                              <  2:700563/139‑1 (MQ=255)
           aTCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTa                                                                                                                      <  1:310083/139‑1 (MQ=255)
                               gCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAggg                                                                                                  >  2:1201663/1‑139 (MQ=255)
                                 ggTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                                    <  1:257410/119‑1 (MQ=255)
                                 ggTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                                    >  2:257410/1‑119 (MQ=255)
                                    gAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCg                                                                                             >  2:694328/1‑139 (MQ=255)
                                       ccAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACg                                                                                          >  1:1473329/1‑139 (MQ=255)
                                        cAGGATGAAGAGTATCGACTTCGGTCGGGTCCGGAGAAGCCGAATGATGAAGCTTATCAACGCGATG‑AGGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGc                                                                                         <  2:119188/139‑1 (MQ=255)
                                         aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                                    <  2:1504516/111‑1 (MQ=255)
                                         aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                                    >  1:1504516/1‑111 (MQ=255)
                                         aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        >  1:129348/1‑139 (MQ=255)
                                         aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCCAATGATGACGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        <  2:1426723/139‑1 (MQ=255)
                                            aTGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                     <  2:1473329/139‑1 (MQ=255)
                                                   gTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              <  1:1415956/139‑1 (MQ=255)
                                                   gTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              <  1:694328/139‑1 (MQ=255)
                                                    tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGAATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATC‑AGGGCATCGACGCCCAATCTGTTac                                                                             >  2:93986/1‑139 (MQ=255)
                                                    tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                                 >  2:481553/1‑87 (MQ=255)
                                                    tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                                 <  1:481553/87‑1 (MQ=255)
                                                    tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTac                                                                             <  2:906921/139‑1 (MQ=255)
                                                              ggtcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              <  2:754505/128‑1 (MQ=255)
                                                              ggtcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              >  1:754505/1‑128 (MQ=255)
                                                                             agcCGAATGATGAAGCTTATCAACGCGATGAATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATC‑AGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                    <  1:93986/139‑1 (MQ=255)
                                                                              gcCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   >  2:1429297/1‑139 (MQ=255)
                                                                              gcCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   >  1:1277175/1‑139 (MQ=255)
                                                                               ccGAATTAAGAAGCTTATCAAAGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCCAGGGCAACGGCCCCCAATCTGTTACCCTTCTGGAACAATTTTTTGATGAGCa                                                  <  1:1141720/139‑1 (MQ=255)
                                                                                       tgaAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTAAACTTCTGGAACAATTTTTTGATGAGAAGGCATTga                                          >  1:819748/1‑139 (MQ=255)
                                                                                        gaAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTCTGATGAGCAGGCATTgag                                         <  2:835246/139‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGa                                                                                                                             <  1:36425/54‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGa                                                                                                                             >  2:36425/1‑54 (MQ=255)
                                                                                         aaGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                       >  2:397471/1‑76 (MQ=255)
                                                                                         aaGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                       <  1:397471/76‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        >  2:450162/1‑91 (MQ=255)
                                                                                         aaGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        <  1:450162/91‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              <  1:843713/101‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              >  2:843713/1‑101 (MQ=255)
                                                                                               aTCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCCTTGGGAGTAGTTAAGcc                                                                                                                  <  2:862313/59‑1 (MQ=255)
                                                                                               aTCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                                  <  1:1139530/59‑1 (MQ=255)
                                                                                               aTCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                                  >  2:1139530/1‑59 (MQ=255)
                                                                                               aTCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                                  >  1:862313/1‑59 (MQ=255)
                                                                                                tCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                    >  1:555826/1‑120 (MQ=255)
                                                                                                tCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                    <  2:555826/120‑1 (MQ=255)
                                                                                                 cAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                                <  1:1429297/139‑1 (MQ=255)
                                                                                                    cgcgATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATc                                                                                   >  1:34399/1‑85 (MQ=255)
                                                                                                    cgcgATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATc                                                                                   <  2:34399/85‑1 (MQ=255)
                                                                                                    cgcgATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   >  2:662853/1‑117 (MQ=255)
                                                                                                    cgcgATG‑ATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   <  1:662853/117‑1 (MQ=255)
                                                                                                         tg‑atgaATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   <  2:1266317/112‑1 (MQ=255)
                                                                                                         tg‑atgaATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   >  1:1266317/1‑112 (MQ=255)
                                                                                                                 tataATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        <  2:1446960/68‑1 (MQ=255)
                                                                                                                 tataATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        >  1:1446960/1‑68 (MQ=255)
                                                                                                                    aaTCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                                >  2:1173769/1‑121 (MQ=255)
                                                                                                                    aaTCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                                <  1:1173769/121‑1 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAacgac               >  2:1256359/1‑136 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAacgac               <  1:1256359/136‑1 (MQ=255)
                                                                                                                      tCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCa            >  2:1297098/1‑139 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                     <  2:1495857/63‑1 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                     >  1:1495857/1‑63 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCAt                                             <  2:1014240/103‑1 (MQ=255)
                                                                                                                         gcagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCAt                                             >  1:1014240/1‑103 (MQ=255)
                                                                                                                          cagTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACttt        <  2:129348/139‑1 (MQ=255)
                                                                                                                           agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        <  1:965117/58‑1 (MQ=255)
                                                                                                                           agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        >  2:965117/1‑58 (MQ=255)
                                                                                                                           agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTac                                                                             >  2:667846/1‑69 (MQ=255)
                                                                                                                           agTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTac                                                                             <  1:667846/69‑1 (MQ=255)
                                                                                                                               gcgcTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                     >  2:1416678/1‑57 (MQ=255)
                                                                                                                               gcgcTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                     <  1:1416678/57‑1 (MQ=255)
                                                                                                                                cgcTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTGGtttt  >  1:813908/1‑139 (MQ=255)
                                                                                                                                  cTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCa                                                  <  2:129952/89‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
ACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATG‑ATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTGGTTTT  >  NC_000913/4018628‑4018893

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