Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,001,645 Δ5 bp coding (220‑224/951 nt) corA → magnesium/nickel/cobalt transporter

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4001645 4001649 5 43 [0] [0] 45 corA magnesium/nickel/cobalt transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 40016440 (0.000)43 (0.950) 37/282 0.2 100% coding (219/951 nt) corA magnesium/nickel/cobalt transporter
?NC_000913 4001650 = 0 (0.000)coding (225/951 nt) corA magnesium/nickel/cobalt transporter

TCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGT  >  NC_000913/4001515‑4001792
                                                                                                                                       |                                                                                                                                              
attgatcttgtcgaaccggacgacgacgagcgactgcgcgaacaatctgaacttggccagggcatggcaaaccgcccggaactggaagacatcgaagcatcggcacgttgctttgaggacgacgacggcctgcatCTCcttct                                                                                                                                         >  1:283893‑M2/136‑143 (MQ=255)
 ttgctcttgtcgaaccggacgccgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggagctggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttctt                                                                                                                                        <  2:481881‑M2/9‑1 (MQ=255)
     tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttc                                                                                                                                    >  2:50823‑M2/131‑143 (MQ=255)
       ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                                  >  1:712465‑M2/129‑143 (MQ=255)
       ttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                                  >  1:592542‑M2/129‑143 (MQ=255)
       gtgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCcttcttcttctt                                                                                                                                  <  1:592584‑M2/15‑1 (MQ=255)
        tgtcgaaccggacgacgacgagcgcctgcgcgtacaatctgaagttggccagagcctggcaacccgcccggaacgggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCttt                                                                                                                                 <  2:164836‑M2/16‑1 (MQ=255)
          tcgcatcggacgccggcgggcgactgcgcgtacaatctgaccttggcccgagcctgtcaacccgcccggagctggaagacatcgaagcatcggcacgtgtctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGa                                                                                                                               <  2:663893‑M2/18‑1 (MQ=255)
           cgacccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                              <  2:188121‑M2/19‑1 (MQ=255)
                        gacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAagagc                                                                                                                 >  1:727452‑M2/112‑141 (MQ=255)
                                   gcgcgtacaatctgaacttggctagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgcagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAc                                                                                                      <  2:697599‑M2/43‑1 (MQ=255)
                                   gcccgtactatcgtaacttggccagagccgggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAc                                                                                                      <  2:226500‑M2/43‑1 (MQ=255)
                                    cgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACt                                                                                                     <  2:510490‑M2/44‑1 (MQ=255)
                                     gcgtacaatctgaccttggcccgagcctggcaacccgctccgcactggaagtcatcgaagcatctgcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGCAGATCACGCCGGTAACTc                                                                                                    <  2:494203‑M2/45‑1 (MQ=255)
                                       gtacgatctgagcttggccggagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTATGAAGACGCGGAAGATCACGCCGGTAACTCCa                                                                                                  <  2:101732‑M2/47‑1 (MQ=255)
                                       gtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTGCTTCTTTGAAGATGCGGAAGATCAAGACGGTAACTCCa                                                                                                  >  1:366863‑M2/97‑143 (MQ=255)
                                         aaaatctgaacttgggcagagcctggcgacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACt                                                                                                <  2:41277‑M2/49‑1 (MQ=255)
                                            atctgaactttgccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                             <  1:618314‑M2/52‑1 (MQ=255)
                                            atctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                             <  2:278434‑M2/52‑1 (MQ=255)
                                            acctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACtgtg                                                                                             <  1:688033‑M2/52‑1 (MQ=255)
                                             tctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttggttggagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGCTGCGGAAGATCACGCCGGGAACTGCGGTGTgg                                                                                            >  1:452920‑M2/91‑143 (MQ=255)
                                                    ttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAATATCACGCCGGTAACTCCACTGTGGCATTTAc                                                                                     >  1:499593‑M2/84‑143 (MQ=255)
                                                      ggccagagcctggcaacccgcccggaactggaagtcagcgaagcatgggcacgtggctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTGACTCCACTGTGGCATTTACCa                                                                                   >  1:830085‑M2/82‑143 (MQ=255)
                                                            agcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCAGTGTGGCATTTACCATCCGTg                                                                             >  1:74450‑M2/76‑143 (MQ=255)
                                                              cctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAAGTCCACTGTGGCATTTACCATCCGTGAt                                                                           >  1:388106‑M2/74‑143 (MQ=255)
                                                                         gcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                                <  1:496529‑M2/81‑1 (MQ=255)
                                                                         gccccgaactggaagacatcgcagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                                <  2:408265‑M2/81‑1 (MQ=255)
                                                                                   ggaagactccggagcatcggcatgtttctctgacgacgacgacggcctgcatCTCCTTCTTCTCCTTTTAAGATGCGGAAGATCACGCCGGTAACTCCACTGGGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgt                                                      <  2:815607‑M2/91‑1 (MQ=255)
                                                                                      agacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgag                                                   >  2:119862‑M2/50‑143 (MQ=255)
                                                                                       gacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtgcgc                                                  >  1:24300‑M2/49‑143 (MQ=255)
                                                                                            cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTCACTCCACTGTGGCATTTACCATCCGTGATGGTCGCCTGTTTGCTctgcgcgagcgtgca                                             >  2:462620‑M2/44‑141 (MQ=255)
                                                                                            cgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTCGCATTTCCAATCCGGGATGTTTTTCTGTTTACTCTGCGTGAGCGTGaa                                             <  1:350669‑M2/100‑1 (MQ=255)
                                                                                              aagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGtgaaat                                           >  1:691097‑M2/42‑141 (MQ=255)
                                                                                                 catcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGGCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                        <  2:157765‑M2/105‑1 (MQ=255)
                                                                                                 catcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                        <  1:644599‑M2/105‑1 (MQ=255)
                                                                                                         cgttgctgtgaagacgacgacggcctgaatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCTTGAGCGTGAACTGCCCGCTTTTc                                <  1:428415‑M2/113‑1 (MQ=255)
                                                                                                          gtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGACCTGCCCGCTTTTCg                               >  1:130819‑M2/30‑143 (MQ=255)
                                                                                                              ctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGGAACTCCACTGTGGCATTTACCATCCGGGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCGCCTTTTCGTCTg                           >  1:326199‑M2/26‑143 (MQ=255)
                                                                                                                ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGCACTGCCCGCTTTTCGGCCGTa                         >  1:570359‑M2/24‑143 (MQ=255)
                                                                                                                ttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTa                         <  1:50823‑M2/120‑1 (MQ=255)
                                                                                                                  gaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATc                       <  1:235195‑M2/122‑1 (MQ=255)
                                                                                                                                cctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTCCTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTGTGTATCGAATGCGTGCCCGt         >  1:340054‑M2/8‑143 (MQ=255)
                                                                                                                                   tcaccTCCTTCTGCGTCTTTGAAGCTGCGCAAGATCACGCCGGTAACTCCACTGTGGCATGTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGc      <  2:657178/139‑1 (MQ=255)
                                                                                                                                   gcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGCCTGTGTCGTATGCGTGCCCGCAGc      >  1:139367‑M2/5‑143 (MQ=255)
                                                                                                                                       cTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGt  >  1:481114/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                              
TCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGT  >  NC_000913/4001515‑4001792

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ 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.