Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,555,491 Δ41,523 bp [yjiC]lgoD 37 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4555491 4597013 41523 34 [0] [0] 33 [yjiC]–lgoD [yjiC],ytiC,ytiD,idlP,iraD,hypT,iadA,yjiG,yjiH,kptA,yjiJ,yjiK,ytiA,yjiL,yjiM,yjiN,mdtM,rpnD,yjiR,yjiS,yjiT,yjiV,mcrC,mcrB,symE,symR,hsdS,hsdM,hsdR,mrr,yjiA,yjiX,btsT,tsr,lgoT,lgoR,lgoD

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45554900 (0.000)31 (1.060) 26/96 0.0 100% coding (831/831 nt) yjiC uncharacterized protein YjiC
?NC_000913 4597014 = 0 (0.000)intergenic (+2/+136) lgoD/opgB L‑galactonate oxidoreductase/phosphoglycerol transferase I

TTTAGCTATCTGTATAACCTGAATGTTAGTACTCATTCTTCCTGGTAGTTATTTACCAATATAATTCCATTCACCATTTTTTAATTCAAACAGTTCTG  >  NC_000913/4555442‑4555539
                                                |                                                 
tttAGCTATCTGTATAACCTGAATGTTAGTACTCATTCTTCCTGGTAGTa                                                  <  1:1285326‑M1/50‑2 (MQ=255)
tttAGCTATCTGTATAACCTGAATGTTAGTACTCATTCTTCCTGGTAGT                                                   <  1:2229774‑M1/49‑1 (MQ=255)
     cTATCTGTATAACCTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggc                                              >  1:1261282‑M1/1‑44 (MQ=255)
          tGTATAACCTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagat                                         >  1:2035560‑M1/1‑39 (MQ=255)
            tataACCTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattat                                      >  1:750454‑M1/1‑37 (MQ=255)
            tataACCTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattat                                      >  1:291748‑M1/1‑37 (MQ=255)
            tataACCTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagatta                                       <  1:984962‑M1/49‑13 (MQ=255)
            tataACCTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagat                                         <  1:1028474‑M1/47‑11 (MQ=255)
              taACCTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattatt                                     <  1:226440‑M1/49‑15 (MQ=255)
                 ccTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaag                                 <  1:657650‑M1/50‑19 (MQ=255)
                  cTGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagt                                <  1:398327‑M1/50‑20 (MQ=255)
                   tGAATGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagtt                               >  1:1921462‑M1/1‑30 (MQ=255)
                     aaTGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttgc                             >  1:723326‑M1/1‑28 (MQ=255)
                     aaTGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttgc                             >  1:1278019‑M1/1‑28 (MQ=255)
                     aaTGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttg                              <  1:841406‑M1/49‑22 (MQ=255)
                      aTGTTAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttgc                             >  1:1972760‑M1/1‑27 (MQ=255)
                         ttAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttgccatg                         >  1:2033095‑M1/1‑24 (MQ=255)
                         ttAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttgccat                          <  1:1630561‑M1/49‑26 (MQ=255)
                         ttAGTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttgccat                          <  1:940057‑M1/49‑26 (MQ=255)
                            gTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttgccatgcag                      <  1:321554‑M1/50‑30 (MQ=255)
                            gTACTCATTCTTCCTGGTAGTaaggcgagattattaaagttgccatgca                       <  1:1468067‑M1/49‑29 (MQ=255)
                               cTCATTCTTCCTGGTAGTaaggcgagattattaaagttgccatgcagcgt                   >  1:1281035‑M1/1‑18 (MQ=255)
                                tCATTCTTCCTGGTAGTaaggcgagattattaaagttgccatgcagcgtc                  >  1:1630957‑M1/1‑17 (MQ=255)
                                  aTTCTTCCTGGTAGTaaggcgagattattaaagttgccatgcagcgtcc                 <  1:1850667‑M1/49‑35 (MQ=255)
                                   ttcttcCTGGTAGTaaggcgagattattaaagttgccatgcagcgtccgg               >  1:1082802‑M1/1‑14 (MQ=255)
                                      ttcCTGGTAGTaaggcgagattattaaagttgccatgcagcgtccgggg             >  1:902518‑M1/1‑11 (MQ=255)
                                        ccTGGTAGTaaggcgagattattaaagttgccatgcagcgtccggggaa           <  1:1972730‑M1/49‑41 (MQ=255)
                                            gtagtaaggcgagattattaaagttgccatgcagcgtccggggaagtgtt      <  1:1943192‑M1/50‑46 (MQ=255)
                                            gtagtaaggcgagattattaaagttgccatgcagcgtccggggaagtgtt      >  1:735533‑M1/1‑5 (MQ=255)
                                            gtagtaaggcgagattattaaagttgccatgcagcgtccggggaagtgtt      <  1:1611936‑M1/50‑46 (MQ=255)
                                            gtagtaaggcgagattattaaagttgccatgcagcgtccggggaagtgt       <  1:428435‑M1/49‑45 (MQ=255)
                                               gtaaggcgagattattaaagttgccatgcagggtccggggaagtgttgcg   >  1:1210396‑M1/1‑2 (MQ=255)
                                                taaggcgagattattaaagttgccatgcagcgtccggggaagtgttgggc  >  1:1155278‑M1/1‑1 (MQ=255)
                                                taaggcgagattattaaagttgccatgcagcgtccggggaagtgttgggc  <  1:2071333‑M1/50‑50 (MQ=255)
                                                |                                                 
TTTAGCTATCTGTATAACCTGAATGTTAGTACTCATTCTTCCTGGTAGTTATTTACCAATATAATTCCATTCACCATTTTTTAATTCAAACAGTTCTG  >  NC_000913/4555442‑4555539

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ 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.