Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 3,888,720 Δ6 bp intergenic (+211/‑5) tnaC → / → tnaA tnaAB operon leader peptide/tryptophanase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 38887193 (0.090)20 (0.590) 18/484 0.7 87.0% intergenic (+210/‑11) tnaC/tnaA tnaAB operon leader peptide/tryptophanase
?NC_000913 3888726 = 3 (0.090)intergenic (+217/‑4) tnaC/tnaA tnaAB operon leader peptide/tryptophanase

TTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGC  >  NC_000913/3888551‑3888974
                                                                                                                                                                               |                                                                                                                                                                                                                                                        
attagattcaatgtgatctatttgtttgctatatcttaattttgccttttgcaaaggtcctctctcgtttatttacttgttttagtaaatgatggtgcttgcatatatatctggcgaattaatcggtatagcagctgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTcgc                                                                                                                                                                               >  1:151611‑M2/176‑251 (MQ=255)
     attcaatgtgatctatttgtttgctatatcttaattttgccttttgcaaaggtcatctctcgtttatttacttgttttagtaaatgatggtgcttgcatatatatctggcgaattaatcggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCtt                                                                                                                                                                           >  1:43751‑M2/171‑250 (MQ=255)
         actgtggtctatttgtttgctatgtctgaattttgccttttgcaaaggtcatctctcgtttatttacttgttttagtaaatgatggtgcttgcatatatatctggcgaattaatcggtatagcagatgaaatattcatagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATATCACTGAACCATTCAGCATTCGTGTTATTGAGCCAGTAACACGTACCACTCGCGATTATCGt                                                                                                                                                                      >  2:271172‑M2/167‑251 (MQ=255)
                                        tttgccttttgcaaaggtcatctctcgtttatttacttgttttagtaaatgatggtgcttgcatatatctctcgcgaattaatcggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTACTGGAAAACTTTAAACATCACCCTGAACCGTTCCGCATTCGTGTTATTGAGCAAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTCTGAAcc                                                                                                                                       >  1:39257‑M2/136‑251 (MQ=255)
                                                ttgcataggccatctctcgtttatttacttgttgtagtaaatgatggtgcttgcatatatatctggcgaattaatcggtagagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCGCTGAACCTTTCCGCATCCGTGTGATTTAGCCAGTGAAACGGACCTCTCGCGCTTATCGTGAAGACGCATTGATTGAATCCGGTGTGAACCCGTTCct                                                                                                                                <  2:340718‑M2/123‑1 (MQ=255)
                                                    cATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGCAATTAAAATAAATGAAGGATTATGTAATGGAAAACCTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAAGCCGGTATGAACCCGTTCctgctg                                                                                                                            <  2:125694/250‑1 (MQ=255)
                                                      aggtcatctctcgtttatttacttgttttaataaatgatggtgcttgcatatatatctggcgaattaatcggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGAt                                                                                                                         >  2:103467‑M2/122‑251 (MQ=255)
                                                                  gtttatttacttgttttagtaaatgatggtgcttgcatatatatctggcgaattaatcggtatagcagatgtaatattcacagggatcactgtaattaaaagaaatgaaTTTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATTGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTGCCTGCTGGATAGCGAAGatgct                                                                                                             <  2:230477‑M2/142‑3 (MQ=255)
                                                                     tatttacttgttttagtaaatgatggtgcttgcatatatatctggcgaattaatcggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGttttt                                                                                                          >  2:182478‑M2/107‑251 (MQ=255)
                                                                      ttGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCa                                                                                                                                                                                    >  2:97587/1‑176 (MQ=255)
                                                                      ttGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCa                                                                                                                                                                                    <  1:97587/176‑1 (MQ=255)
                                                                       tttacttgttttagtaaatgatggtgcttgaatatatctctggcgaagtcatcggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAACCTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAACACGTACAACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGCAGATGTTTTTAt                                                                                                        >  1:325093‑M2/105‑251 (MQ=255)
                                                                                    ataaatgatggtgcttgcatatatatctggcgaattaatcggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTGCTGAcc                                                                                           <  1:103467‑M2/160‑1 (MQ=255)
                                                                                           atggtgcttgcatatatatctggcgaattaatcggtatagcagctgtaatattcacagtgatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTAACCGACAGCg                                                                                    >  1:122815‑M2/85‑251 (MQ=255)
                                                                                                gcttgcatatatatctggcgaattaatcggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCAcc                                                                               >  2:52320‑M2/80‑251 (MQ=255)
                                                                                                              ctggcgacttaatcgggatagcagatgtaatattcacagggatcactgtacttaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTCTTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGTTTTACTGACCGACAGCGGCACCGGGGCGGTGACgcc                                                                 >  2:141647‑M2/66‑250 (MQ=255)
                                                                                                                ggcgaattaatcggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTAGCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCaga                                                               <  2:43751‑M2/188‑1 (MQ=255)
                                                                                                                            ggtatagcagatgtaatattcacagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTACTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTg                                                   <  1:246520‑M2/200‑1 (MQ=255)
                                                                                                                                  gcagaggtaatcttcacagggatccctgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCACTTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCgatg                                              >  2:19779‑M2/46‑250 (MQ=255)
                                                                                                                                                  cagggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGc                             >  1:123094‑M2/30‑251 (MQ=255)
                                                                                                                                                    gggatcactgtaattaaaataaatgaaTGTAATGGAAAACTTTAAACATCTCCCTGAACAGATCCGCATTAGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATAGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGAAAAGCATGCAGGCTGCGATGATGCGCGGAGACGAGGCCt                           >  2:307186‑M2/28‑251 (MQ=255)
                                                                                                                                                     ggatcaatgtaattaaaataaatgaaTGTACTGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGGATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCCGCGAGGATGCGCGGCGACGAAGCCt                           >  2:113208‑M2/27‑250 (MQ=255)
                                                                                                                                                                             aaTGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGc  >  1:155113‑M2/3‑251 (MQ=255)
                                                                                                                                                                               |                                                                                                                                                                                                                                                        
TTCAATGTGATCTATTTGTTTGCTATATCTTAATTTTGCCTTTTGCAAAGGTCATCTCTCGTTTATTTACTTGTTTTAGTAAATGATGGTGCTTGCATATATATCTGGCGAATTAATCGGTATAGCAGATGTAATATTCACAGGGATCACTGTAATTAAAATAAATGAAGGATTATGTAATGGAAAACTTTAAACATCTCCCTGAACCGTTCCGCATTCGTGTTATTGAGCCAGTAAAACGTACCACTCGCGCTTATCGTGAAGAGGCAATTATTAAATCCGGTATGAACCCGTTCCTGCTGGATAGCGAAGATGTTTTTATCGATTTACTGACCGACAGCGGCACCGGGGCGGTGACGCAGAGCATGCAGGCTGCGATGATGCGCGGCGACGAAGCCTACAGCGGCAGTCGTAGCTACTATGC  >  NC_000913/3888551‑3888974

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