Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 2,666,074 IS1 (+) +9 bp coding (640‑648/1281 nt) csiE → stationary phase‑inducible protein CsiE

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 279163 =NA (NA)7 (0.540) 6/290 0.5 100% noncoding (768/768 nt) IS1 repeat region
?NC_000913 2666074 = 0 (0.000)coding (640/1281 nt) csiE stationary phase‑inducible protein CsiE
* ? NC_000913 = 291401NA (NA)21 (1.620) 18/288 0.0 100% noncoding (1/768 nt) IS1 repeat region
?NC_000913 = 2666082 0 (0.000)coding (648/1281 nt) csiE stationary phase‑inducible protein CsiE

ACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGG  >  NC_000913/279002‑279355
                                                                                                                                                                 |                                                                                                                                                                                                
acacAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAACAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGCTTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGCTTTTGAGAACGAc                                                                                                           >  1:180319/1‑249 (MQ=255)
          gcTGACGGTGGTAATGGCGCGCACAAAGCTGAAAACCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGt                                                                                               >  1:115044/1‑251 (MQ=255)
             gACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTccc                                                                                            >  1:186791/1‑251 (MQ=255)
                                                           aTCGCCGCGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATGTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTc                                                                                                                                                              <  1:179968/139‑1 (MQ=255)
                                                                 gggCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCGCCACCGATTGTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGc                                                        <  1:219022/235‑1 (MQ=255)
                                                                             ggCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCAGCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGc                            <  1:110203/251‑1 (MQ=255)
                                                                              gcagcaGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATAGATTTGGCCCGGTAATGACTCCACCTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCt                           >  1:276362/1‑251 (MQ=255)
                                                                                        ttgCGTAATGGTACCCAGGGAACCGGCGCGCAGGCTAAGTGAGGCTTCATCCAGCGATTTATCGATTTGGCCCGTTAATGACTCCAACGTATTGATAGTGTTTTATGTTCAGAt                                                                                                                                                          <  1:287026/113‑1 (MQ=37)
                                                                                                       cAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGTTAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAgg  >  1:91198/1‑251 (MQ=255)
                                                                                                             aCCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGttn                                                                                                                                                                     <  1:79682/82‑2 (MQ=255)
                                                                                                                            gagtgagGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCgcct                                                                                                    >  1:135143/1‑129 (MQ=255)
                                                                                                                              gtgagGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAgg                                                                                >  1:17072/1‑150 (MQ=255)
                                                                                                                                                      cGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTNNGAGAACNNc                                                                                                           >  2:1508/1‑99 (MQ=255)
                                                                                                                                                              cccGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGCCTTTGTCATGCTGCTCCCCCGATTTTGAGCACGACAGCGCCTTc                                                                                                  >  1:321569/1‑100 (MQ=11)
                                                                                                                                                                 |                                                                                                                                                                                                
ACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGG  >  NC_000913/279002‑279355

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