Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 2,400,257 IS5 (+) +4 bp 58.3% 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.390) 5/272 NT 53.8% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 2400257 = 5 (0.340)intergenic (‑137/‑783) lrhA/alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2
* ? CP009273 = 2283592NA (NA)8 (0.520) 8/272 NT 60.8% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 2400260 5 (0.340)intergenic (‑140/‑780) lrhA/alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2

CATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACT  >  CP009273/2400118‑2400408
                                                                                                                                              |                                                                                                                                                    
gtataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccaagaagaatttccccgcttattcgcaccttccCTaaa                                                                                                                                                     <  1:109377‑M2/5‑1 (MQ=255)
  tatatTATCACTTACTGGCGGCCCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGa                                                                                                                                              <  2:228241/149‑1 (MQ=255)
     acgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgtctgagtcagctgagaagaatttccccgcttatgcgcaccttccCTAAAAAAGCGACTc                                                                                                                                           <  2:64237‑M2/15‑1 (MQ=255)
     acgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTc                                                                                                                                           <  2:238758‑M2/15‑1 (MQ=255)
                                    ggTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGAt                                                                                                            <  2:277068/149‑1 (MQ=255)
                                          tGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                                                                                      >  2:216092‑M2/1‑101 (MQ=255)
                                                gCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggag                                                                                                     <  1:208839‑M2/144‑50 (MQ=255)
                                                  acacacaAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATcg                                                                                              >  2:91306/1‑148 (MQ=255)
                                                      acaaacaaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatggttggcatctggagccatagaacag                                                                                          >  2:167533‑M2/1‑89 (MQ=255)
                                                        aaacaaGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaataagtccctgatatgagatcatgtttgtcatctggagccatagaacaggg                                                                                        >  2:179487‑M2/1‑87 (MQ=255)
                                                               aaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCat                                                                                 >  2:176487‑M2/77‑149 (MQ=255)
                                                                  aatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgttctttccCTAAAAATGCGTCTCACATTGTTTCGTTATATTGCCTTAAGTAGATCACAGAATATATCTTCAGGGATCGCATATc                                                                              <  2:54534‑M2/76‑1 (MQ=255)
                                                                            aaTGTGCCTCCCGGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTa                                                                    <  2:67034/149‑1 (MQ=255)
                                                                                tgCCTCCCCGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacacgtccctgatatgagatcatgtttgtcatctggagccatagaaccgggttcatcatgagtcatcaac                                                                     >  1:54276‑M2/1‑63 (MQ=255)
                                                                                             gTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTc                                                        >  1:97686/1‑144 (MQ=255)
                                                                                                    ccGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcag                                            >  2:229068‑M2/1‑43 (MQ=255)
                                                                                                       agccgagaagaatttccccgcttattcgcaccttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATAACTATTAAGTTACTCACTCTTTTCTATTTATGAc                                              >  1:35239‑M2/37‑144 (MQ=255)
                                                                                                                cGAAATGAAAAATGCATTTGAGTTGCACTAAggaaggtgcgaagaagtccctgatatgagctcatgtttggcctctggagccatagaacagggttcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaag                                     >  1:208230‑M2/1‑31 (MQ=255)
                                                                                                                  atttccccgcttattcgcaccttccCTAAAAAAACGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGtat                              >  2:168848‑M2/26‑149 (MQ=255)
                                                                                                                                   gcaaagcactaAggaaggcgcgaacaactccctcatatgagatcatgtttgccatctggagccatagaaccttggtcatcatgagtcatcaacttaccttcgccgacagtgaattcagcagtaagcgccgtcagaccagaaaag                  >  1:66096‑M2/6‑12 (MQ=255)
                                                                                                                                      cttccCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTTTCTTTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCtt          <  2:150242‑M2/144‑1 (MQ=255)
                                                                                                                                         cACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACCTGCGCGTGTTTGTTTAAATGTAAATGTGAGTcc            >  1:43936/1‑144 (MQ=255)
                                                                                                                                          aCTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTTTTCATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTc      <  2:185338/149‑1 (MQ=255)
                                                                                                                                             aaaaaaGCGACTCACGTTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGGGTCCTTGTTCCAc   >  2:38446/1‑149 (MQ=255)
                                                                                                                                              aaaaaGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGTGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCAct  <  2:120096/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTCCACT  >  CP009273/2400118‑2400408

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.