Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 2,400,257 IS5 (+) +4 bp 59.1% 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)17 (0.860) 16/270 NT 64.7% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 2400257 = 9 (0.470)intergenic (‑137/‑783) lrhA/alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2
* ? CP009273 = 2283592NA (NA)9 (0.450) 8/270 NT 49.2% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 2400260 9 (0.470)intergenic (‑140/‑780) lrhA/alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2

GCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCC  >  CP009273/2400110‑2400405
                                                                                                                                                    |                                                                                                                                                    
aaaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                      <  1:350927‑M2/1‑1 (MQ=255)
aaaaacgggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                      <  1:156013‑M2/1‑1 (MQ=255)
 aaaactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCt                                                                                                                                                     <  1:204974‑M2/2‑1 (MQ=255)
   aactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTaa                                                                                                                                                   <  1:119986‑M2/4‑1 (MQ=255)
   aactggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTaa                                                                                                                                                   >  1:79782‑M2/146‑149 (MQ=255)
   aactggggagaacgccttaaatggcgaacaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTaa                                                                                                                                                   <  1:273239‑M2/4‑1 (MQ=255)
      tggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccc                                                                                                                                                      >  2:263008‑M2/143‑143 (MQ=255)
        cATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaagg                                                                                                                                             >  1:340179‑M2/1‑143 (MQ=255)
         ggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGc                                                                                                                                             <  1:96316‑M2/10‑1 (MQ=255)
            taacgccttaaatggcgaagaaacggtctaaataggctgatgcaaggcatttacgggcgaaaaaatcggcgcaaacatgaagaaatgaaatgactgagtcaaccgagaagaattgacccgcttatacgcaccttccctacatt                                                                                                                                                >  2:305019‑M2/137‑139 (MQ=255)
              cgatcaTTTACTGGCGGCTCATACTGGGCGGGTTAATGCTGGGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAGCCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaa                                                                                                                                       >  1:95073‑M2/3‑137 (MQ=255)
                  cttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACAt                                                                                                                                    >  1:202375‑M2/131‑149 (MQ=255)
                        tggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAACAAAGCGACTCACAt                                                                                                                                    >  2:223142‑M2/125‑143 (MQ=255)
                           gcggcTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTAGCCGTTTCATCCGGAACCCGGCCGAGAGTAAACATGCATCTGCGTTGCACTCAggaagtagcgcgcaagtccctgcca                                                                                                                          <  1:77129‑M2/149‑26 (MQ=255)
                                    tattgtgtGGTTAAATGCTTTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGGtata                                                                                                                        >  2:250362/8‑143 (MQ=255)
                                     gtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagcagagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATaa                                                                                                                       >  2:36126‑M2/112‑143 (MQ=255)
                                                 tgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATcaca                                                                                                     >  1:351062‑M2/100‑149 (MQ=255)
                                                   attcaaggcatgtacgggagaaaaaagcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGa                                                                                                   <  1:2192‑M2/52‑1 (MQ=255)
                                                     tcaaggcatttatgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATca                                                                                                       <  2:293761‑M2/48‑1 (MQ=255)
                                                          gcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATc                                                                                            <  1:223099‑M2/59‑1 (MQ=255)
                                                          cacacacaAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGaa                                                                                                  <  2:124560/143‑1 (MQ=255)
                                                               acaaacaaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAg                                                                                       >  1:242720/1‑149 (MQ=255)
                                                                    caaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATc                                                                                  <  1:76061/149‑1 (MQ=255)
                                                                     agaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGACGTACCTCACAAAATATATCTTCAgggaact                                                                                 >  1:174856‑M2/80‑146 (MQ=255)
                                                                     agaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCg                                                                                 >  1:358416‑M2/80‑149 (MQ=255)
                                                                      gaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATACCTTCAgg                                                                                      >  2:315033‑M2/79‑143 (MQ=255)
                                                                           ttCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgagaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatga                                                                           >  1:246210‑M2/1‑77 (MQ=255)
                                                                             ccATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                                                                               >  2:305955‑M2/1‑75 (MQ=255)
                                                                                 ttGTAACGTGCCTCCCGGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatga                                                                           >  2:283866‑M2/1‑71 (MQ=255)
                                                                                    aaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGtt                                                                  >  1:98876‑M2/65‑149 (MQ=255)
                                                                                       tgtgCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCAAATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTact                                                               >  1:253442/1‑149 (MQ=255)
                                                                                            cTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaactt                                                                >  2:228474‑M2/1‑60 (MQ=255)
                                                                                                          cATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGa                                            <  1:79974/149‑1 (MQ=255)
                                                                                                           aTCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACCAACATTGAGCCGGTATAATGCCTGCAGTACATAGATCTCAACAACTTCATATAGCGCAGATCTATTAAGTTACTCACTCTTTTCTAtt                                                 >  2:232585/1‑143 (MQ=255)
                                                                                                             gtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTa                                               <  2:316977‑M2/104‑1 (MQ=255)
                                                                                                                 gccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGAc                                           >  2:21424‑M2/36‑143 (MQ=255)
                                                                                                                        aCGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGCGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGcgcg                                    <  2:79630/143‑1 (MQ=255)
                                                                                                                          gcaagtAAAAATGAATATGAGTTGCACTACggaaggtgcgaaccagtccctgatatgagctcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgt                            >  1:351497‑M2/3‑30 (MQ=255)
                                                                                                                               ccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGt                             >  2:307289‑M2/22‑143 (MQ=255)
                                                                                                                                aaaaaTGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGta                            <  2:20176/143‑1 (MQ=255)
                                                                                                                                    cttattcgcaccttccCTAAAAAAGCGGCTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTAtaaa                        >  2:277459‑M2/17‑143 (MQ=255)
                                                                                                                                          cgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAg            >  1:281688‑M2/11‑149 (MQ=255)
                                                                                                                                               ctccactAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttc       >  1:126281‑M2/4‑9 (MQ=255)
                                                                                                                                                ttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAg            >  2:182285‑M2/5‑143 (MQ=255)
                                                                                                                                                    cTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggggccagagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagattttcttgtc  <  1:346766‑M2/149‑146 (MQ=255)
                                                                                                                                                    cTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagatttt        >  2:274302‑M2/1‑4 (MQ=255)
                                                                                                                                                    cTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaagagatttt        >  2:320690‑M2/1‑4 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCC  >  CP009273/2400110‑2400405

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