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→TAArmuCDNA recombination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/37);  total (32/37)

ACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTGGTTTT  >  NC_000913/4018628‑4018893
                                                                                                                                     |                                                                                                                                     
aCGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                                 <  1:289614/139‑1 (MQ=255)
   cgAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTggg                                                                                                                              <  2:700563/139‑1 (MQ=255)
           aTCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTa                                                                                                                      <  1:310083/139‑1 (MQ=255)
                               gCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAggg                                                                                                  >  2:1201663/1‑139 (MQ=255)
                                 ggTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                                    <  1:257410/119‑1 (MQ=255)
                                 ggTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                                    >  2:257410/1‑119 (MQ=255)
                                    gAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCg                                                                                             >  2:694328/1‑139 (MQ=255)
                                       ccAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACg                                                                                          >  1:1473329/1‑139 (MQ=255)
                                        cAGGATGAAGAGTATCGACTTCGGTCGGGTCCGGAGAAGCCGAATGATGAAGCTTATCAACGCGATGAGGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGc                                                                                         <  2:119188/139‑1 (MQ=255)
                                         aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                                    <  2:1504516/111‑1 (MQ=255)
                                         aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAg                                                                                                                    >  1:1504516/1‑111 (MQ=255)
                                         aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        >  1:129348/1‑139 (MQ=255)
                                         aGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCCAATGATGACGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        <  2:1426723/139‑1 (MQ=255)
                                            aTGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCaa                                                                                     <  2:1473329/139‑1 (MQ=255)
                                                   gTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              <  1:1415956/139‑1 (MQ=255)
                                                   gTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              <  1:694328/139‑1 (MQ=255)
                                                    tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGAATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCAGGGCATCGACGCCCAATCTGTTac                                                                             >  2:93986/1‑139 (MQ=255)
                                                    tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                                 >  2:481553/1‑87 (MQ=255)
                                                    tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCAtt                                                                                                                                 <  1:481553/87‑1 (MQ=255)
                                                    tATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTac                                                                             <  2:906921/139‑1 (MQ=255)
                                                              ggtcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              <  2:754505/128‑1 (MQ=255)
                                                              ggtcggtTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              >  1:754505/1‑128 (MQ=255)
                                                                             agcCGAATGATGAAGCTTATCAACGCGATGAATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                    <  1:93986/139‑1 (MQ=255)
                                                                              gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   >  2:1429297/1‑139 (MQ=255)
                                                                              gcCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   >  1:1277175/1‑139 (MQ=255)
                                                                               ccGAATTAAGAAGCTTATCAAAGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCCAGGGCAACGGCCCCCAATCTGTTACCCTTCTGGAACAATTTTTTGATGAGCa                                                  <  1:1141720/139‑1 (MQ=255)
                                                                                       tgaAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTAAACTTCTGGAACAATTTTTTGATGAGAAGGCATTga                                          >  1:819748/1‑139 (MQ=255)
                                                                                        gaAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTCTGATGAGCAGGCATTgag                                         <  2:835246/139‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGa                                                                                                                             <  1:36425/54‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGa                                                                                                                             >  2:36425/1‑54 (MQ=255)
                                                                                         aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                       >  2:397471/1‑76 (MQ=255)
                                                                                         aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATc                                                                                                       <  1:397471/76‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        >  2:450162/1‑91 (MQ=255)
                                                                                         aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGcc                                                                                        <  1:450162/91‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              <  1:843713/101‑1 (MQ=255)
                                                                                         aaGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTa                                                                              >  2:843713/1‑101 (MQ=255)
                                                                                               aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCCTTGGGAGTAGTTAAGcc                                                                                                                  <  2:862313/59‑1 (MQ=255)
                                                                                               aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                                  <  1:1139530/59‑1 (MQ=255)
                                                                                               aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                                  >  2:1139530/1‑59 (MQ=255)
                                                                                               aTCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGcc                                                                                                                  >  1:862313/1‑59 (MQ=255)
                                                                                                tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                    >  1:555826/1‑120 (MQ=255)
                                                                                                tCAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAg                                                    <  2:555826/120‑1 (MQ=255)
                                                                                                 cAACGCGATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGAt                                <  1:1429297/139‑1 (MQ=255)
                                                                                                    cgcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATc                                                                                   >  1:34399/1‑85 (MQ=255)
                                                                                                    cgcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATc                                                                                   <  2:34399/85‑1 (MQ=255)
                                                                                                    cgcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   >  2:662853/1‑117 (MQ=255)
                                                                                                    cgcgATGATGAATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   <  1:662853/117‑1 (MQ=255)
                                                                                                         tgatgaATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   <  2:1266317/112‑1 (MQ=255)
                                                                                                         tgatgaATATAATCAGCAGTCGCGCTAACCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGc                                                   >  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)
                                                                                                                                     |                                                                                                                                     
ACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAAGAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGAACAATTTTTTGATGAGCAGGCATTGAGATGGTGGATAAGTCACAAGAAACGACGCACTTTGGTTTT  >  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: