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‑‑‑GATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGT  >  NC_000913/3966573‑3966866
                                                                                                                                                  |                                                                                                                                                      
ggcgtactggagatattgcaggatggatttggtttcctccgttccgcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTc                                                                                                                                                    <  2:356445‑M2/5‑1 (MQ=255)
ggcgtaatggagatattgcaggatggatttggcttcctccgttcagcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTc                                                                                                                                                    >  1:58751‑M2/147‑151 (MQ=255)
 gcgtactggagatattgcaggatggatttggtttcctccgttccgcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTcg                                                                                                                                                   >  2:394234‑M2/146‑151 (MQ=255)
  cgtactggagatattgcaggatggatttggtttcctccgttccgcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTcg                                                                                                                                                   >  2:87199‑M2/145‑150 (MQ=255)
   gtactggagatattgcaggatggatttggtttcctccgttccgcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTcgcc                                                                                                                                                 >  2:220762‑M2/144‑151 (MQ=255)
      ctggagatattgcaggatggatttggtttcctccgttccgcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTcgcccg                                                                                                                                               <  1:337085‑M2/10‑1 (MQ=255)
       tggagatattgcaggatggatttggtttcctccgttccgcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTcgcccgcc                                                                                                                                             >  1:290448‑M2/140‑151 (MQ=255)
                   aggatggatttggtttcctccgttccgcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTg                                                                                                                                  >  2:163165‑M2/128‑150 (MQ=255)
                     gatggatttggtttcctccgttccgcagacagctcctacctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAAc                                                                                                                               >  2:14683‑M2/126‑151 (MQ=255)
                           tttggtttcctccgttccgcagacagctcctccctcgccggtcctgatgacatctacttttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTa                                                                                                                           <  2:88345‑M2/30‑1 (MQ=255)
                                                     ctcctcCCCCGCCGGTCCTGATACCATTTACTTTCCCCCTAGCAAAACCCGCCGTTTCAACCCCCGCACTGGTGATACCATCTCTGGTAAGATGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                                                                      <  1:76613/144‑1 (MQ=255)
                                                       cctccctcgccggtcctgatgacatctacgtttcccctagccaaacccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                                                             <  2:210951‑M2/60‑1 (MQ=255)
                                                          ccctcgccggtcctgatgacatttgcttttcccctagcaaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                                                                             <  1:159651‑M2/60‑1 (MQ=255)
                                                          acctcgccggtcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                                                                          >  1:123659‑M2/89‑151 (MQ=255)
                                                              cgccggtcctgatgacatctacgtttcccctagccaaatccgccgcttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                                                                      >  2:79436‑M2/85‑151 (MQ=255)
                                                                    tcctgatgacatctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                                                                >  2:213640‑M2/79‑151 (MQ=255)
                                                                              atctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                                                                      >  2:282407‑M2/69‑151 (MQ=255)
                                                                               tctacgtttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                                                                     >  1:132042‑M2/68‑151 (MQ=255)
                                                                                     tttcccctagccaaatccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                                                                  <  2:237295‑M2/87‑1 (MQ=255)
                                                                                     tttcccctagcaaaacccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                                                                      <  2:269088‑M2/83‑1 (MQ=255)
                                                                                     ttccccctaccaaaacccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                                                                   <  2:96828‑M2/86‑1 (MQ=255)
                                                                                                 aaacccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt                                                     <  2:244372‑M2/100‑1 (MQ=255)
                                                                                                 aaaTCAGCAGTTTCAACCTCCCCACGGGTGCTACCCTACCTGGTAAGACGATTCTCCCGCCGCAAGAAGGTGAAAGCTATTTGGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACATGAAATCGACCGCAACAAAATCCTCTTTg                                                    >  1:292268/1‑150 (MQ=255)
                                                                                                    tccgccgtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga                                                 >  2:311580‑M2/47‑150 (MQ=255)
                                                                                                          gtttcaacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAAc                                          >  2:278436‑M2/41‑151 (MQ=255)
                                                                                                              caacctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                      >  1:144639‑M2/37‑151 (MQ=255)
                                                                                                                 cctccgcactggtgataccatctctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTg                                   >  1:276272‑M2/34‑151 (MQ=255)
                                                                                                                                       tctggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGt              >  1:122542‑M2/12‑150 (MQ=255)
                                                                                                                                         tggtaagatGTTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTAt            <  2:246019‑M2/141‑1 (MQ=255)
                                                                                                                                          ggtaagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTAt            <  1:179437‑M2/141‑1 (MQ=255)
                                                                                                                                          ggtaagatGATTCCCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCg               <  1:225833‑M2/138‑1 (MQ=255)
                                                                                                                                             aagatGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAc       <  2:77245‑M2/146‑1 (MQ=255)
                                                                                                                                               gatGATACCCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACg      <  2:279580‑M2/147‑1 (MQ=255)
                                                                                                                                                  gATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt  >  2:389679/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAA‑‑‑GATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGT  >  NC_000913/3966573‑3966866

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.