Predicted mutation
evidence seq id position mutation freq annotation gene description
JC NC_000913 3,998,555 Δ3 bp 9.5% coding (573‑575/2163 nt) uvrD → DNA‑dependent ATPase I and helicase II

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3998558 =42 (0.960)4 (0.100) 4/240 NT 9.5% coding (576/2163 nt) uvrD DNA‑dependent ATPase I and helicase II
?NC_000913 = 3998554 37 (0.900)coding (572/2163 nt) uvrD DNA‑dependent ATPase I and helicase II

TACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTGGCGGGCGACACCG  >  NC_000913/3998415‑3998703
                                                                                                                                               |                                                                                                                                                 
tggtacatcaacagccagaaagatgaaggcctgcgtccgcatcatattcaaagctacggtaatccggtggagcagacctggcagaaggtgtatcaggcgtatcaggaagcgtgcgaccgcgcgggcctggtggacttcgccgagctgct                                                                                                                                              <  2:77517‑M1/6‑1 (MQ=255)
 gcATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGctgctgctg                                                                                                                                             >  1:23085/2‑149 (MQ=255)
  ggacatcaacagccagaaagatgaaggcctgcgtccgcatcatattcaaagctacggtaatccggtggagcagacctggcagaaggtgtatcaggcgtatcaggaagcgtgtgaccgcgcgggcctggtggacttcgccgaGCtgctgc                                                                                                                                            <  1:359663‑M1/8‑1 (MQ=255)
    tCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCACATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTgcg                                                                                                                                           >  3:104935/1‑148 (MQ=255)
       aCAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCt                                                                                                                                       >  1:45531/1‑149 (MQ=255)
             aGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAg                                                                                                                                 >  1:201038/1‑149 (MQ=255)
              gAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGt                                                                                                                                >  2:287155/1‑149 (MQ=255)
                 aGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTtgt                                                                                                                             <  2:51687/149‑1 (MQ=255)
                  gATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTtgt                                                                                                                             <  4:237914/148‑1 (MQ=255)
                     gAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCt                                                                                                                         >  2:275171/1‑149 (MQ=255)
                      aaGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGATCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACAGCGCGGGCATGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCt                                                                                                                         >  4:241672/1‑148 (MQ=255)
                                    catcatATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCata                                                                                                           >  3:110164/1‑148 (MQ=255)
                                     atcatATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCatat                                                                                                          >  4:323234/1‑148 (MQ=255)
                                     atcatATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCatat                                                                                                          >  3:400212/1‑148 (MQ=255)
                                     atcatATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCatat                                                                                                          <  4:183680/148‑1 (MQ=255)
                                          catattcaaagctacggtaatccggtggagcagacctggcagaaggtgtatcaggcgtatcaggaagcgtgtgaccgcgcgggcctggtggacttcgccgagctg                                                                                                                                                >  3:241132‑M1/102‑105 (MQ=255)
                                          catattcaaagctacggtaatccggtggagcagacctggcagaaggtgtatcaggcgtatcaggaagcgtgtgaccgcgcgggcctggtggacttcgccgagctg                                                                                                                                                <  4:241132‑M1/4‑1 (MQ=255)
                                              aaaGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAacac                                                                                                >  1:256168/1‑149 (MQ=255)
                                               tcaaagctacggtaatccggtggagcagacctggcagaaggtgtatcaggcgtatcaggaagcgtgtgaccgcgcgggcctggtggacttcgccgaGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAacac                                                                                                <  4:135874‑M1/52‑1 (MQ=255)
                                                 gCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTa                                                                                              >  4:126966/1‑148 (MQ=255)
                                                            ccGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGttt                                                                                  <  1:198444/149‑1 (MQ=255)
                                                            ccGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGtt                                                                                   <  3:126966/148‑1 (MQ=255)
                                                              tccggtggagcagacctggcagaaggtgtatcaggcgtatcaggaagcgtgtgaccgcgcgggcctggtggacttcgccgaGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACcgc                                                                                          <  2:79034‑M1/58‑1 (MQ=255)
                                                              tccggtggagcagacctggcagaaggtgtatcaggcgtatcaggaagcgtgtgaccgcgcgggcctggtggacttcgccgaGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACcgc                                                                                          >  1:79034‑M1/82‑139 (MQ=255)
                                                                      aGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATc                                                                        >  2:101317/1‑149 (MQ=255)
                                                                      aGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGCGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATc                                                                        <  1:77517/149‑1 (MQ=255)
                                                                         ccTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCtg                                                                     <  2:201038/149‑1 (MQ=255)
                                                                               aGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGAc                                                               >  1:68944/1‑149 (MQ=255)
                                                                                   gaaggtgtatcaggcgtatcaggaagcgtgtgaccgcgcgggcctggtggacttcgccgaGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGaa                                                            >  3:135874‑M1/61‑148 (MQ=255)
                                                                                       tATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTcc                                                        <  3:5543/148‑1 (MQ=255)
                                                                                       tATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCa                                                       >  1:360003/1‑149 (MQ=255)
                                                                                               gTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATAcc                                                <  4:110164/148‑1 (MQ=255)
                                                                                               gTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACca                                               <  1:101317/149‑1 (MQ=255)
                                                                                                tATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACcaa                                              <  1:320718/149‑1 (MQ=255)
                                                                                                tATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACca                                               <  4:249060/148‑1 (MQ=255)
                                                                                                   cAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACcaaca                                            <  3:341334/148‑1 (MQ=255)
                                                                                                      gAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACAt                                        >  2:5333/1‑149 (MQ=255)
                                                                                                        aGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTc                                      <  2:256168/149‑1 (MQ=255)
                                                                                                            tgtgACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGt                                   >  4:60760/1‑148 (MQ=255)
                                                                                                               gACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTAcg                                <  3:60760/148‑1 (MQ=255)
                                                                                                                 ccGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGt                             <  1:269105/149‑1 (MQ=255)
                                                                                                                         gCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGc                     <  2:359663/149‑1 (MQ=255)
                                                                                                                            tggtggACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCctg                  >  2:254819/1‑149 (MQ=255)
                                                                                                                                 gACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTgg              >  4:396174/1‑148 (MQ=255)
                                                                                                                                    ttCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTggcggg          >  1:43103/1‑149 (MQ=255)
                                                                                                                                    ttCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTggcggg          >  2:147554/1‑149 (MQ=255)
                                                                                                                                     tCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTggcgggc         <  1:275171/149‑1 (MQ=255)
                                                                                                                                       gCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGaa                                                            <  4:111094/96‑1 (MQ=255)
                                                                                                                                       gCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGaa                                                            >  3:111094/1‑96 (MQ=255)
                                                                                                                                       gCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTGGCGGGCGa       <  1:402773/149‑1 (MQ=255)
                                                                                                                                            gctgctgctGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTGGCGGGCGACACCg  <  2:23085/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TACATCAACAGCCAGAAAGATGAAGGCCTGCGTCCGCATCATATTCAAAGCTACGGTAATCCGGTGGAGCAGACCTGGCAGAAGGTGTATCAGGCGTATCAGGAAGCGTGTGACCGCGCGGGCCTGGTGGACTTCGCCGAGCTGCTGCTGCGCGCTCACGAGTTGTGGCTTAACAAGCCGCATATCCTGCAACACTACCGCGAACGTTTTACCAATATCCTGGTGGACGAATTCCAGGATACCAACAACATTCAGTACGCGTGGATCCGCCTGCTGGCGGGCGACACCG  >  NC_000913/3998415‑3998703

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.