Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 3,311,380 IS4 (+) +12 bp noncoding (156‑167/174 nt) psrO → novel sRNA, function unknown

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3311380 =0 (0.000)27 (1.350) 17/132 0.1 100% noncoding (156/174 nt) psrO novel sRNA, function unknown
?NC_000913 4502090 = NA (NA)noncoding (1426/1426 nt) IS4 repeat region
* ? NC_000913 = 33113910 (0.000)17 (0.850) 13/132 0.3 100% noncoding (167/174 nt) psrO novel sRNA, function unknown
?NC_000913 = 4503515 NA (NA)noncoding (1/1426 nt) IS4 repeat region

GTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCC  >  NC_000913/4502024‑4502157
                                                                  |                                                                   
gctcatcgagcgcctgggtctgcgtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTa                                                                    >  2:588155‑M2/67‑68 (MQ=255)
  tcatcgagcgcctgggtctgcgtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAAt                                                                  >  2:28401‑M2/65‑68 (MQ=255)
  tcatcgagcgcctgggtctgcgtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAAt                                                                  <  2:647486‑M2/4‑1 (MQ=255)
  aaactcGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATATGCATCTTGGTTCTGTAAt                                                                  >  2:143289/4‑68 (MQ=255)
  aTTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAAt                                                                  <  1:445362/68‑1 (MQ=255)
     ctcGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGcc                                                               >  1:564610/1‑68 (MQ=255)
       gagcgcctgggtctgcgtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAAtgccga                                                             >  2:587985‑M2/60‑66 (MQ=255)
       gagcgcctgggtctgcgtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAATGCcag                                                             >  1:416196‑M2/60‑68 (MQ=255)
       cGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCAGATTGGTTCTGTAATGCcag                                                             >  1:381717/1‑68 (MQ=255)
             ctgggtctgcgtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAATGCCAGTCAGtt                                                       <  1:690766‑M2/15‑1 (MQ=255)
             ctgggtctgcgtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAATGCCAGTCAGtt                                                       <  1:517279‑M2/15‑1 (MQ=255)
             ctgggtctgcgtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAATGCCAGTCAGtt                                                       <  1:398421‑M2/15‑1 (MQ=255)
             ctgggtctgcgtcgctaattcttgcgagtctcagaaaagggggcctgagtggcTAATGCCAGTCAGtt                                                       <  1:458801‑M2/15‑1 (MQ=255)
             cAACCAGATAATTTATCAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGtt                                                       <  1:641300/68‑1 (MQ=255)
             cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGtt                                                       <  1:39712/68‑1 (MQ=255)
             cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGtt                                                       <  2:547805/68‑1 (MQ=255)
             cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGtt                                                       <  1:561020/68‑1 (MQ=255)
             cAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCAGCTTGGTTCTGTAATGCCAGTCAGtt                                                       <  1:525188/68‑1 (MQ=255)
              aaCCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTa                                                      >  2:448030/1‑68 (MQ=255)
                      aaTTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAactg                                              >  1:716739/1‑68 (MQ=255)
                       gtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAATGCCAGTCAGtt                                                       >  1:751073‑M2/44‑58 (MQ=255)
                       gtcgctaattcttgcgagtttcagaaaagggggcctgagtggcTAATGCCAGTCAGtt                                                       <  2:751074‑M2/15‑1 (MQ=255)
                       aTTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAactga                                             <  2:114790/68‑1 (MQ=255)
                       aTTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAactga                                             >  2:630942/1‑68 (MQ=255)
                         ttAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGtt                                                       <  1:149793/56‑1 (MQ=255)
                         ttAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGtt                                                       >  2:149793/1‑56 (MQ=255)
                          taccAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTAAAGCAacagactg                                          >  1:67167/4‑68 (MQ=255)
                          tAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAactgaatg                                          >  2:330764/1‑68 (MQ=255)
                            aCAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGc                                        >  1:656691/1‑68 (MQ=255)
                            aCAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGc                                        >  1:561698/1‑68 (MQ=255)
                            aCAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGc                                        >  2:519473/1‑68 (MQ=255)
                              aGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAAt                                                                  >  1:76180/1‑40 (MQ=255)
                              aGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAAt                                                                  <  2:76180/40‑1 (MQ=255)
                                 cttgcgagtttcagaaaagggggcctgagtggcTAATGCCAGTCAGTTAAGCAACTGACTGGCTCttt                                   >  1:465513‑M2/34‑68 (MQ=255)
                                  ttgcgagtttcagaaaagggggcctgagtggcTAATGCCAGTCAGTTAAGCAACTGACTGGCTCtttt                                  <  1:154352‑M2/36‑1 (MQ=255)
                                      cTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCgg                              >  1:95475/1‑68 (MQ=255)
                                        gtttcagaaaagggggcctgagtggcTACTGCCAGACAGATAAGCAACTGCCTGGCTCTATTTAgggg                            >  1:169194‑M2/27‑68 (MQ=255)
                                        gtttcagaaaagggggcctgagtggcTAATGCCAGTCAGTTAGGCAACTGACTGGCTCTTTTTAgggg                            >  1:170040‑M2/27‑68 (MQ=255)
                                         tttcagaaaagggggcctgagtggcTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGc                           <  2:626033‑M2/43‑1 (MQ=255)
                                         tttcagaaaagggggcctgagtggcTAATGCCAGTCAGTTAAGCAACTGACTGGCGCTTTTTCGGGGc                           <  2:578468‑M2/43‑1 (MQ=255)
                                          ttcagaaaagggggcctgagtggcTAATGCCTGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCt                          <  2:117424‑M2/44‑1 (MQ=255)
                                                aaagggggcctgagtggcTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTgggg                    <  2:321482‑M2/50‑1 (MQ=255)
                                                aaagggggcctgagtggcTAATACCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTgggg                    <  2:623322‑M2/50‑1 (MQ=255)
                                                aTAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTgggg                    <  1:35162/68‑1 (MQ=255)
                                                  aGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTa                  <  1:476022/68‑1 (MQ=255)
                                                   gggggcctgagtggcTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTAt                 >  2:272199‑M2/16‑68 (MQ=255)
                                                    ggggcctgagtggcTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTAtt                <  2:446520‑M2/54‑1 (MQ=255)
                                                     cATCTTGGTGCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTAttt               <  2:135682/68‑1 (MQ=255)
                                                      aTCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTc              <  2:745658/68‑1 (MQ=255)
                                                         ttGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAg           >  2:81852/1‑68 (MQ=255)
                                                           gagtggcTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAggg         <  2:84368‑M2/61‑1 (MQ=255)
                                                           gagtggcTAATGCCAGTCAGGTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAggg         <  2:626749‑M2/61‑1 (MQ=255)
                                                             ttCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGTGTATTTCCAGGGcc       <  1:513165/68‑1 (MQ=255)
                                                              tCTGTAATGCCAGTCAGTAAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGcat      >  2:408976/1‑66 (MQ=255)
                                                                tGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCctct    <  1:448218/68‑1 (MQ=255)
                                                                gcTTATGCCAGGCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGTGTATTTCCAGGGCctct    <  1:354130‑M2/66‑1 (MQ=255)
                                                                gcTAATGCCTGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGTGCTGTGGGGTATTTCCAGGGCctct    <  2:409762‑M2/66‑1 (MQ=255)
                                                                cgTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGCATTTCCAGGGCctct    <  1:597733/67‑1 (MQ=255)
                                                                 cTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGGGGGGTATTTCCAGGGCCtctc   <  2:218312‑M2/67‑1 (MQ=255)
                                                                  tAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTcc  <  1:289543/68‑1 (MQ=255)
                                                                  |                                                                   
GTATTCTCGTTGGCAACCAGATAATTTAACAGTTGATGCTGAAACGAAATAGGCATCTTGGTTCTGTAATGCCAGTCAGTTAAGCAACTGACTGGCTCTTTTTCGGGGCTGTGGGGTATTTCCAGGGCCTCTCC  >  NC_000913/4502024‑4502157

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