Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 2,406,778 IS5 (+) +4 bp intergenic (‑137/‑780) 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 2288919 =NA (NA)21 (1.090) 17/278 0.1 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 2406778 = 0 (0.000)intergenic (‑137/‑783) lrhA/alaA DNA‑binding transcriptional dual regulator LrhA/glutamate‑‑pyruvate aminotransferase AlaA
* ? NC_000913 = 2290113NA (NA)18 (0.930) 14/278 0.2 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 = 2406781 0 (0.000)intergenic (‑140/‑780) lrhA/alaA DNA‑binding transcriptional dual regulator LrhA/glutamate‑‑pyruvate aminotransferase AlaA

ACACACAGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATT  >  NC_000913/2406534‑2407009
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
ggttgctgaaaaacgataaccaactggcgatgttattcacgctggccaacctgtttcgggcggcccaaatgatacgtcagtgggagagatctccctaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTaaaa                                                                                                                                                                                                                                     >  1:23056‑M2/245‑250 (MQ=255)
                                                       tcgggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAatgac                                                                                                                                                                                                       <  1:477207‑M2/36‑3 (MQ=255)
                                                         gggcggaccaaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggactaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGt                                                                                                                                                                                                                <  1:6781‑M2/27‑1 (MQ=255)
                                                                 caaatgatacgtcagtgggagagatctcactaaaaactggggataacgccttacatggcgaagaaccggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGAc                                                                                                                                                                                                                              >  1:208725‑M2/180‑192 (MQ=255)
                                                                      aggtcgaggtTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCTGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaagg                                                                                                                                                                                                                                 >  1:157584‑M2/1‑178 (MQ=255)
                                                                              agtgggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGtt                                                                                                                                                      <  1:233508‑M2/85‑1 (MQ=255)
                                                                                 aaTTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACCCAAACAAGAAATTCCATTTGTAATGTGCATCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaaaaagtccctgatatgcgatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaactt                                                                                                                                                    <  1:440462‑M2/250‑84 (MQ=255)
                                                                                  ggagagatctcactaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTactc                                                                                                                                                  >  1:218582‑M2/163‑251 (MQ=255)
                                                                                              ctaaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAg                                                                                                                                                                                                                                  >  1:421806‑M2/151‑159 (MQ=255)
                                                                                                          gataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATgcg                                                                                                                          <  1:253893‑M2/113‑1 (MQ=255)
                                                                                                                              ggcTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccaga                                                                                                       <  1:346250‑M2/250‑129 (MQ=255)
                                                                                                                                         aaataggctgagtcaaggcatttacgggagaaaaaatcggcgcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGAGCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGAc                                                                                                                               <  1:190739‑M2/108‑1 (MQ=255)
                                                                                                                                           cTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttc                                                                                           >  1:159393‑M2/1‑109 (MQ=255)
                                                                                                                                             ggTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTCCACTAAggaaggttcgaacaagtccctgatatgag                                                                                                                                                                                                          >  1:153441‑M2/1‑107 (MQ=255)
                                                                                                                                                tAATGATGTGCACACACAAACAAGCAATTCCATTCGTAAAGTGCCTCCCTAGCAGTTTCATCGGGAACCCGAACGAAAGTAAAAATGCATATGAGTTGCACTAAggacggtgcgaacaagtccctgatatgagatcatgattgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgcattcagcagtaagcgccgtcagaccagaaaagagatttaattgtccc                                                                                    <  1:193545‑M2/251‑148 (MQ=255)
                                                                                                                                                            caaacaaacaagcaATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggagcgaacaagtccc                                                                                                                                                                                                                   <  1:65204‑M2/109‑21 (MQ=255)
                                                                                                                                                                      gaaaaaatcggctcaaacatgaagaaatgaaatgactgaatcagccgagaagaatttacccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATACATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTc                                                              >  1:335632‑M2/79‑250 (MQ=255)
                                                                                                                                                                                gctcaaacatgaagaaaagaaatgactgaggcagccgagaagaatgtctccgcttattcgcaccttccCTAAAAAAGCGACGCACATTGGTCCGTTAGAAGGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCGCTCTTTTc                                                                                                                                         <  1:497058‑M2/98‑1 (MQ=255)
                                                                                                                                                                                    aaacatgaagaaatgaaatgactgagtcagccgagaagactttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATAACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACACGAACACGt                                                >  1:79638‑M2/65‑250 (MQ=255)
                                                                                                                                                                                          naagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAAc                                          >  1:311195‑M2/59‑249 (MQ=255)
                                                                                                                                                                                           ccTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgantcannnnnttnn                                                                                                                                                  <  2:306801‑M2/146‑86 (MQ=255)
                                                                                                                                                                                           ccTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagn                                                                                                                     >  1:306801‑M2/1‑61 (MQ=255)
                                                                                                                                                                                           ccTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcg                                                                                                                                             <  1:136592‑M2/151‑91 (MQ=255)
                                                                                                                                                                                                   gaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGa                                                      >  1:259597‑M2/50‑229 (MQ=255)
                                                                                                                                                                                                   aaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGANNTAGatcnn                                                                                                                                                                                          >  2:33727‑M2/50‑96 (MQ=255)
                                                                                                                                                                                                   aaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTcc                                                                                     <  1:33727‑M2/149‑1 (MQ=255)
                                                                                                                                                                                                                 aaCCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggct                   >  1:673‑M2/1‑39 (MQ=255)
                                                                                                                                                                                                                         ngaAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc                                                                                                                                                            <  1:311003‑M2/105‑76 (MQ=255)
                                                                                                                                                                                                                         cGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc                                                                                                                                                            >  2:311003‑M2/1‑31 (MQ=255)
                                                                                                                                                                                                                         cGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgataggaggtcatgtttgtcatctggagccatagaacagggttcatcatgagacatcaacttaccttcgccgacagtgaattcagcagtaagcgccttcagaccagaaaagagattttctgggcccgcatggggcagattctgccagggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatgg           >  1:353384‑M2/1‑31 (MQ=255)
                                                                                                                                                                                                                          gAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggc          >  1:408864‑M2/1‑30 (MQ=255)
                                                                                                                                                                                                                             aGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggagcagattctgccatggcaaaacatggtggaagtcatcgagccgttttaccccaaggctggtaatggccgg       >  1:31284‑M2/1‑27 (MQ=255)
                                                                                                                                                                                                                               cccctcttattcgcatcttcccTAAAAAAGCGACTCACATTGTTCCGTTAGAATGCCTGAAGTAGATCACAGAAGATATCTTCAGGGATCGGATAGCTAGTAAGTTACTCACTCTTTTCTATTTATGACATGCGTTTGTTTGTATAAATGGAAAGTGGAGTCCGGGTTCCACTCACGTGCAGCAGCGCTGGTCATACGCGAACACGTACCAAc                                          <  1:236309/192‑1 (MQ=255)
                                                                                                                                                                                                                                 ccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGGAGATCACAGAATATATCTTCAGGGATCGAATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCAtt  >  1:152996‑M2/20‑251 (MQ=255)
                                                                                                                                                                                                                                        atatGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatganncatcaannt                                                                                                                                                    <  2:25440‑M2/99‑84 (MQ=255)
                                                                                                                                                                                                                                        atatGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtcccgcatggag                                                                            >  1:25440‑M2/1‑16 (MQ=255)
                                                                                                                                                                                                                                         tcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTtactcan                                                                                                                                                <  1:317566‑M2/91‑2 (MQ=255)
                                                                                                                                                                                                                                         tcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTactcac                                                                                                                                                >  2:317566‑M2/12‑102 (MQ=255)
                                                                                                                                                                                                                                          cgcaccttccCTAAAAAAGCGACTCACATTNTTCCGNTNTNNTGNNTGAAGTAg                                                                                                                                                                                               <  2:323404‑M2/44‑1 (MQ=255)
                                                                                                                                                                                                                                          cgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATATTGCCTGAAGTAg                                                                                                                                                                                               >  1:323404‑M2/11‑54 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
ACACACAGCGGCAGCTGCGGCAGCAAAAGTGTTCAGATCGGCAACAGCAACAAATGTTCTCAGCAGATCGAGGTCGAGGTTAATTATCGGACGATTTGCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACTCTCGTGCAGCATCGCTGGTCATACGCGAACACGTACCAACAGCAATGGTGAGGCGGCATCAAGAGCAGGGATCCGTCATT  >  NC_000913/2406534‑2407009

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.