Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 2,400,257 IS5 (+) +4 bp 43.5% intergenic (‑137/‑780) lrhA ← / → alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 2282398 =NA (NA)6 (0.290) 6/276 NT 30.9% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 2400257 = 13 (0.660)intergenic (‑137/‑783) lrhA/alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2
* ? CP009273 = 2283592NA (NA)14 (0.690) 13/276 NT 51.1% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 2400260 13 (0.660)intergenic (‑140/‑780) lrhA/alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2

TTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCC  >  CP009273/2400124‑2400405
                                                                                                                                     |                                                                                                                                                    
ttatCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTCATGTACCTCCCTGGCAATTTACTTTATTTTCTTGACGAAAGTAAAATTGCATATGAGTTGCACTAAGAAAGCGACTca                                                                                                                                       >  2:244757/1‑149 (MQ=255)
    ttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACAtt                                                                                                                                   >  1:203284‑M2/130‑149 (MQ=255)
    cccTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACCAGCAATTACATTTGTAATATGCCCCCCTGGCAGTTTCACCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACAtt                                                                                                                                   <  2:160634/147‑1 (MQ=255)
               ggcTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggcaggtgcgaacaagtccctgatatg                                                                                                                        <  2:168396‑M2/149‑28 (MQ=255)
                 cTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAAt                                                                                                                      >  2:158506/1‑149 (MQ=255)
                 aaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAAt                                                                                                                      <  1:333068‑M2/33‑1 (MQ=255)
                          aactaggctgatgcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGt                                                                                                             <  1:24561‑M2/42‑1 (MQ=255)
                               gTTAATGCTGTGCACACACAAACAAGCAATTAGCATTGTAATGTGCCTCGCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccccgataggagatcaggtttgtagt                                                                                                        <  2:175010‑M2/149‑44 (MQ=255)
                                gctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATca                                                                                                       >  2:175819‑M2/102‑149 (MQ=255)
                                    tGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                                                                                   <  2:290327‑M2/149‑49 (MQ=255)
                                    tGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                                                                                   >  1:67640‑M2/1‑101 (MQ=255)
                                             atttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCtt                                                                                          <  1:187589‑M2/61‑1 (MQ=255)
                                              acacaAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                                                                         <  2:209116‑M2/149‑59 (MQ=255)
                                                acaaaccaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAg                                                                                       <  2:153722/149‑1 (MQ=255)
                                                 caaacagGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAgg                                                                                      <  1:54257/149‑1 (MQ=255)
                                                   gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATCTATCTTCcgggc                                                                                    >  2:321671‑M2/83‑148 (MQ=255)
                                                      aaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCg                                                                                 >  2:141175/1‑149 (MQ=255)
                                                               cATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtc                                                                        >  2:191938‑M2/1‑74 (MQ=255)
                                                                 tttGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgaggcat                                                                      <  2:330178‑M2/149‑78 (MQ=255)
                                                                          tgCCGCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAgggaggtgcgagcaagtccctgatatgcgatcctgtttgtcatctggagccatataacggggttcatcatgagtcatcaacttacc                                                             >  1:191697‑M2/1‑63 (MQ=255)
                                                                          tgCCGCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGCCTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTtctcc                                                             <  1:118423/149‑2 (MQ=255)
                                                                           gCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTACAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcccctggaaccagaaaacagggctcatcatgagtcataacgctccct                                                            <  1:10365‑M2/149‑88 (MQ=255)
                                                                                  tGGCAGTTTCATCCTGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTACggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaccagggttcatcatgagtcatcaacttaccttcgccga                                                     >  1:315817‑M2/1‑55 (MQ=255)
                                                                                    gCAGTTTCATCCGGAACCCGGACGGAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgaca                                                   >  1:144927‑M2/1‑53 (MQ=255)
                                                                                           caagcggagCCCGGACGAAAGTAAAAATGCATCTGAGTTGCACTAAggaaggcgcgaacaagtccctgatatgagatcatgtttgtcatctggaggcatagaacagggttcatcatgagacatcaacttaccttcgccgacagtgaatt                                            >  1:197348‑M2/5‑46 (MQ=255)
                                                                                                              aGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcag                         >  1:153733‑M2/1‑27 (MQ=255)
                                                                                                              aGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGtaaca                         >  1:103838/1‑146 (MQ=255)
                                                                                                              aGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGGGGTTGTAtaa                         >  2:113261/1‑149 (MQ=255)
                                                                                                               gTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTAtaaa                        <  1:213240/149‑1 (MQ=255)
                                                                                                                    attGCATTTGGTTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTAtaaatgtaa                   <  2:280305/147‑1 (MQ=255)
                                                                                                                    aaTGCATATGAGTTGCACTAACAAAGCGCCTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTAtaaatgtaa                   <  2:5730/149‑1 (MQ=255)
                                                                                                                             gAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagata          <  2:173767‑M2/149‑138 (MQ=255)
                                                                                                                             gAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTc          <  1:88915/149‑1 (MQ=255)
                                                                                                                                     cTAAggaaggtgcgaacaagtccctggtatgagatcatgtttgtcatctggagccatagaacggggctcatgatgagacatcaacttaccttcgacgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtc  <  1:139169‑M2/149‑146 (MQ=255)
                                                                                                                                     cTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTcc  >  1:280513/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
TTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCC  >  CP009273/2400124‑2400405

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