Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 2,406,780 IS1 (–) +8 bp intergenic (‑139/‑774) lrhA ← / → alaA DNA‑binding transcriptional dual regulator LrhA/glutamate‑‑pyruvate aminotransferase AlaA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1978503 =NA (NA)10 (0.770) 9/288 0.2 100% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 2406787 0 (0.000)intergenic (‑146/‑774) lrhA/alaA DNA‑binding transcriptional dual regulator LrhA/glutamate‑‑pyruvate aminotransferase AlaA
* ? NC_000913 = 1979270NA (NA)15 (1.150) 11/290 0.1 100% noncoding (1/768 nt) IS1 repeat region
?NC_000913 2406780 = 0 (0.000)intergenic (‑139/‑781) lrhA/alaA DNA‑binding transcriptional dual regulator LrhA/glutamate‑‑pyruvate aminotransferase AlaA

GCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAG  >  NC_000913/2406580‑2406992
                                                                                                                                                                                                               |                                                                                                                                                                                                             
gagaagcggtgtaagtgaactgcagttgccatgttttacggcagtgagagccgagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtccgtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAAAAAAGCGACTCACATTGTTCCGtt                                                                                                                                                                                             >  1:210713‑M2/201‑226 (MQ=255)
                         nccncttntnngnnctcngtgagagcagagatagcactgnngtccggcggnccttttgccgttacgccccaccccgtcagtagctgaacaggagggacagctganagaaacagannccnctggancnnntnnnaanctccatcatacactaaatcngtnngttggcngcataaccAAAAAAGCGACTCACAtt                                                                                                                                                                                                     <  2:74636‑M2/18‑1 (MQ=255)
                            ccatgttttacggcagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccacgggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAg                                                                                                                                                                                <  1:190088‑M2/39‑1 (MQ=255)
                                         cagtgagagcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAAAAAAGCGACTCACAtt                                                                                                                                                                                                     >  1:74636‑M2/160‑177 (MQ=255)
                                                tttGCACTTATCATATATTATCACTTACTGGCGGCTCATACCGAGCTGGTTAAAGATGTGCCCACACCAACAAGCAATTCCATTTGTAAAGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGAggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctacaccgattttgagaacgacagc                                                                                                                    <  1:180731‑M2/251‑92 (MQ=255)
                                                ngcagagatagcgctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTAt                                                                                                                    >  1:84666‑M2/153‑251 (MQ=255)
                                                        tATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCAGTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttc                                                                                                              >  1:231136‑M2/1‑152 (MQ=255)
                                                         aTCATATATTATCACTTACTGGCGGCACATACTGATCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgtttagttcagataatgcccgatgagtttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcc                                                                                                           <  1:90552‑M2/251‑101 (MQ=255)
                                                               gatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATAGCTTCAGGGATCGCATAGCTATTAAGTTACTCACTCTGTTCTATTTATGACATGcgcg                                                                                                      <  1:215645‑M2/113‑1 (MQ=255)
                                                                       cTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCGTCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttacgttcagataatgccctatgactttgtcatgcagctccaccgattttgagaacgacagcgacttcc                                                                                                             >  1:321886‑M2/1‑137 (MQ=255)
                                                                            ncttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAggg                                                                                                                                                       >  1:205356‑M2/125‑188 (MQ=255)
                                                                                                                  caaacaaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccacc                                                                                                                                      >  1:279563‑M2/1‑94 (MQ=255)
                                                                                                                  caaacaaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgata                                                                                                                                                                                       <  1:100892‑M2/118‑25 (MQ=255)
                                                                                                                      nngannnnngatngaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatccccAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAggg                                                                                                                                                       <  2:205356‑M2/64‑1 (MQ=255)
                                                                                                                      caaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGGTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccan                                                                                                                                       >  1:67595‑M2/1‑90 (MQ=255)
                                                                                                                                   ttGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgtctccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgagttggagaacgacagcgacttccgtcccagccgtgc                                                                                                >  1:249687‑M2/1‑77 (MQ=255)
                                                                                                                                     aacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGt                               >  1:10559‑M2/68‑251 (MQ=255)
                                                                                                                                                cttgagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAAAAAAGCGACTCACATTGTTc                                                                                                                                                                                                 >  1:234063‑M2/57‑78 (MQ=255)
                                                                                                                                                         ctcaaacacaccctcatacactaaatcagtaagttggcagcatcaccAAACAAGCGACTCCCATTGTTCCGTTATAATGc                                                                                                                                                                                      >  1:54320‑M2/48‑80 (MQ=255)
                                                                                                                                                          tcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAACAAAGCGCATCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTc                                                                                                                       >  1:318760‑M2/47‑142 (MQ=255)
                                                                                                                                                                    aCCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctcaattcgctgcgtatatcgcttgctgattacgtgcagctttcccttcaggcgggattcatac  <  1:175313‑M2/249‑206 (MQ=255)
                                                                                                                                                                                           tggcagcatcaccAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAatat                                                                                                                                                                 >  1:6569‑M2/14‑67 (MQ=255)
                                                                                                                                                                                           tggcagcatcaccAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAatat                                                                                                                                                                 <  2:6569‑M2/54‑1 (MQ=255)
                                                                                                                                                                                                    caccAAAAAAGCGACTCACATTGATCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATAACTATTAAGTTACTCACTCTTTTCTATTTATGACATGCTCGTTTTTGGATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTgg                                               <  1:319481‑M2/168‑1 (MQ=255)
                                                                                                                                                                                                    cACTAAAAAAGCggtaatgactccaacttattgatagtgttttatgtgcagataatgccccatggctttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagc                                                                                                     >  1:313360‑M2/1‑12 (MQ=255)
                                                                                                                                                                                                               |                                                                                                                                                                                                             
GCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAG  >  NC_000913/2406580‑2406992

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