Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 3,966,718 +GAT coding (302/1260 nt) rho → transcription termination factor

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 39667152 (0.080)29 (1.180)
+GAT
28/286 0.1 93.6% coding (299/1260 nt) rho transcription termination factor
?NC_000913 3966716 = 2 (0.080)coding (300/1260 nt) rho transcription termination factor

GGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAA‑‑‑GATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAC  >  NC_000913/3966573‑3966861
                                                                                                                                              |                                                                                                                                                     
ggCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattc                                                                                                                                               <  2:356445‑M1/151‑9 (MQ=255)
ggCGTAATGGAGATATTGCAGGATGGATTTGGCTTCCTCCGTTCAGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattc                                                                                                                                               >  1:58751‑M1/1‑143 (MQ=255)
 gCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcg                                                                                                                                              >  2:394234‑M1/1‑142 (MQ=255)
  cGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcg                                                                                                                                              >  2:87199‑M1/1‑141 (MQ=255)
   gTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcc                                                                                                                                            >  2:220762‑M1/1‑140 (MQ=255)
      cTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccg                                                                                                                                          <  1:337085‑M1/150‑14 (MQ=255)
       tGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgcc                                                                                                                                        >  1:290448‑M1/1‑136 (MQ=255)
                   aGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtg                                                                                                                             >  2:163165‑M1/1‑124 (MQ=255)
                     gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaac                                                                                                                          >  2:14683‑M1/1‑122 (MQ=255)
                           tttGGTTTCCTCCGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGATGACATCTACTTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgcta                                                                                                                      <  2:88345‑M1/149‑34 (MQ=255)
                                                     ctcctcCCCCGCCGGTCCTGATACCATTTACTTTCCCCCTAGCAAAACCCGCCGTTTCAACCCCCGCACTGGTGATACCATCTCTGGTAAGATGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                                                                 <  1:76613/144‑1 (MQ=255)
                                                       cctccctcGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaa                                                                                        <  2:210951‑M1/151‑64 (MQ=255)
                                                          ccctCGCCGGTCCTGATGACATTTGCTTTTCCCCTAGCAAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaa                                                                                        <  1:159651‑M1/147‑64 (MQ=255)
                                                          aCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaactt                                                                                     >  1:123659‑M1/1‑85 (MQ=255)
                                                              cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGCTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgac                                                                                 >  2:79436‑M1/1‑81 (MQ=255)
                                                                    tCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacct                                                                           >  2:213640‑M1/1‑75 (MQ=255)
                                                                              aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgccc                                                                 >  2:282407‑M1/1‑65 (MQ=255)
                                                                               tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccg                                                                >  1:132042‑M1/1‑64 (MQ=255)
                                                                                     tttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaa                                                             <  2:237295‑M1/148‑91 (MQ=255)
                                                                                     tttCCCCTAGCAAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgccc                                                                 <  2:269088‑M1/144‑87 (MQ=255)
                                                                                     ttccccctacCAAAACCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgca                                                              <  2:96828‑M1/144‑90 (MQ=255)
                                                                                                 aaaTCAGCAGTTTCAACCTCCCCACGGGTGCTACCCTACCTGGTAAGACGATTCTCCCGCCGCAAGAAGGTGAAAGCTATTTGGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACATGAAATCGACCGCAACAAAATCCTCTTTg                                               >  1:292268/1‑150 (MQ=255)
                                                                                                 aaaCCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctcttt                                                <  2:244372‑M1/149‑104 (MQ=255)
                                                                                                    tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgaga                                            >  2:311580‑M1/1‑43 (MQ=255)
                                                                                                          gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgagaacttaac                                     >  2:278436‑M1/1‑37 (MQ=255)
                                                                                                              cAACCTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgagaacttaaccccg                                 >  1:144639‑M1/1‑33 (MQ=255)
                                                                                                                 ccTCCGCACTGGTGATACCATCTCTGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgagaacttaaccccgctg                              >  1:276272‑M1/1‑30 (MQ=255)
                                                                                                                                       tctGGTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgagaacttaaccccgctgcacgcaaactctcgtctgcgt         >  1:122542‑M1/1‑8 (MQ=255)
                                                                                                                                         tGGTAAgatgtttcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgagaacttaaccccgctgcacgcaaactctcgtctgcgtat       <  2:246019‑M1/150‑145 (MQ=255)
                                                                                                                                          ggTAAgatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgagaacttaaccccgctgcacgcaaactctcgtctgcgtat       <  1:179437‑M1/149‑145 (MQ=255)
                                                                                                                                          ggTAAgatgattcccccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgagaacttaaccccgctgcacgcaaactctcgtctgcg          <  1:225833‑M1/146‑142 (MQ=255)
                                                                                                                                             aagatgattcgcccgccgaaagaaggtgaacgctattttgcgctgctgaaagttaacgaagttaacttcgacaaacctgaaaacgcccgcaacaaaatcctctttgagaacttaaccccgctgcacgcaaactctcgtctgcgtatggaac  <  2:77245‑M1/151‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAA‑‑‑GATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAC  >  NC_000913/3966573‑3966861

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