Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 2,882,823 IS5 (+) +4 bp coding (1313‑1316/1509 nt) casA ← CRISP RNA (crRNA) containing Cascade antiviral complex protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2290113NA (NA)166 (1.790) 82/198 0.0 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 = 2882826 0 (0.000)coding (1313/1509 nt) casA CRISP RNA (crRNA) containing Cascade antiviral complex protein
* ? NC_000913 2882823 =0 (0.000)122 (1.320) 78/198 0.0 100% coding (1316/1509 nt) casA CRISP RNA (crRNA) containing Cascade antiviral complex protein
?NC_000913 = 3131340 NA (NA)noncoding (1195/1195 nt) IS5 repeat region

AAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGAT  >  NC_000913/3131248‑3131439
                                                                                            |                                                                                                   
aaGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATc                                                                                              <  1:1790975/100‑1 (MQ=11)
aaGGCATTTACGGGAGAAAAAATCGGCTCAAACATAAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATc                                                                                              >  1:1249401/1‑100 (MQ=11)
 aGGTATTTACGGTAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATc                                                                                              <  2:1033458/99‑1 (MQ=2)
 aGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCa                                                                                             <  1:692327/100‑1 (MQ=31)
   gCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAgg                                                                                           <  1:1420367/100‑1 (MQ=33)
   gCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAgg                                                                                           <  2:1221426/100‑1 (MQ=33)
    cATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGt                                                                                          >  2:1120563/1‑100 (MQ=35)
     aTTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTc                                                                                         <  2:216044/100‑1 (MQ=35)
        tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATa                                                                                      <  2:1193122/100‑1 (MQ=40)
          cGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACg                                                                                    >  1:318120/1‑100 (MQ=255)
             gagaACAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCtt                                                                                 <  2:1760431/100‑1 (MQ=255)
               gaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCg                                                                               <  2:1201000/100‑1 (MQ=255)
                aaaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCga                                                                              <  2:2112405/100‑1 (MQ=255)
                aaaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCga                                                                              <  1:282412/100‑1 (MQ=255)
                 aaaaaTCGGCTCAAACCTGAAGAAATGAAATGACGGAGTCAGCCGAGAAGAATTTCCGGGCTGATTCGGACCGTCCCTGAATCAGGGCAGACGCTTCgag                                                                             >  2:1951397/1‑100 (MQ=14)
                 aaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCgag                                                                             <  2:216535/100‑1 (MQ=255)
                   aaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGAt                                                                           <  2:350054/100‑1 (MQ=255)
                   aaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGAt                                                                           >  1:479489/1‑100 (MQ=255)
                    aaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATa                                                                          <  1:2162451/100‑1 (MQ=255)
                    aaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATa                                                                          <  1:1338293/100‑1 (MQ=255)
                     aTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATAc                                                                         <  1:1895730/100‑1 (MQ=255)
                     aTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATAc                                                                         >  1:247975/1‑100 (MQ=255)
                      tCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACt                                                                        <  1:1735782/100‑1 (MQ=255)
                         gCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTaa                                                                     >  1:1432967/1‑100 (MQ=255)
                          cTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAAc                                                                    <  2:2001994/100‑1 (MQ=255)
                           tCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACg                                                                   <  2:1202645/100‑1 (MQ=255)
                           tCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACg                                                                   <  1:1606746/100‑1 (MQ=255)
                            cAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGc                                                                  <  1:454600/100‑1 (MQ=255)
                             aaaCATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGcc                                                                 <  1:1081422/100‑1 (MQ=255)
                                cATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCaaa                                                              <  1:2309383/100‑1 (MQ=255)
                                 aTGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCTAAac                                                             <  1:1062203/100‑1 (MQ=255)
                                  tGAAGAAAGGAAATGACTGAGTCAGGCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGGCaaa                                                              <  2:2003234/98‑1 (MQ=255)
                                   gaagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAacac                                                           >  1:381718/1‑100 (MQ=255)
                                    aagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACAcc                                                          >  1:2067715/1‑100 (MQ=255)
                                     agaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCa                                                         >  2:1516596/1‑100 (MQ=255)
                                      gaaatgaaatGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAg                                                        >  1:699907/1‑100 (MQ=255)
                                           gaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGaaa                                                   <  1:57604/100‑1 (MQ=255)
                                           gaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGaaa                                                   <  1:1494400/100‑1 (MQ=255)
                                           gaaatgaCTGAGACAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGaaa                                                   >  1:1536295/1‑100 (MQ=255)
                                                gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGc                                              <  2:63482/100‑1 (MQ=255)
                                                 aCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCg                                             <  1:1620807/100‑1 (MQ=255)
                                                  cTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCgg                                            >  2:415386/1‑100 (MQ=255)
                                                  cTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCgg                                            <  1:1003100/100‑1 (MQ=255)
                                                        cAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGcc                                      <  2:1210946/100‑1 (MQ=255)
                                                        cAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGcc                                      <  1:1120565/100‑1 (MQ=255)
                                                         aGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACACCAAACACCAGCGAAATGAGCGGCTTGCcg                                     <  1:519598/100‑1 (MQ=255)
                                                                aagaaTTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCa                              <  1:879683/100‑1 (MQ=255)
                                                                        ccccGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCATGCCGACGATCAACACcgcg                      <  2:1432963/100‑1 (MQ=255)
                                                                         cccGCTTATTCGCACCTTCCCTAAATCAGGTCCTACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGa                     >  1:1749296/1‑100 (MQ=255)
                                                                         cccGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGa                     >  1:1880048/1‑100 (MQ=255)
                                                                                 ttCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGccacc             >  1:2100306/1‑100 (MQ=255)
                                                                                   cGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAg           >  2:1721650/1‑100 (MQ=255)
                                                                                   cGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAg           >  1:1250143/1‑100 (MQ=255)
                                                                                     cACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGaa         >  1:1834201/1‑100 (MQ=255)
                                                                                          tCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCg    <  1:1283602/100‑1 (MQ=255)
                                                                                          tCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCg    >  1:1295966/1‑100 (MQ=255)
                                                                                           cccTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGACTACCGa   >  2:401456/1‑100 (MQ=255)
                                                                                            ccTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGAt  >  2:203092/1‑100 (MQ=255)
                                                                                            |                                                                                                   
AAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGAT  >  NC_000913/3131248‑3131439

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