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

GCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTC  >  CP009273/2400110‑2400404
                                                                                                                                                   |                                                                                                                                                   
aaaactggggataacgccttaaatggcgaagaaacggtctaaatagactgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCt                                                                                                                                                    >  2:132393‑M2/148‑149 (MQ=255)
     tggggataacgccttaaatggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgaaccttccCt                                                                                                                                                    >  1:154376‑M2/143‑144 (MQ=255)
        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)
                                                                                                                                                   |                                                                                                                                                   
GCACTTATCATATATTATCACTTACTGGCGGCTCATACTGAGCTGGTTAATGCTGTGCACACACAAACAAGCAATTCCATTTGTAATGTGCCTCCCTGGCAGTTTCATCCGGAACCCGGACGAAAGTAAAAATGCATATGAGTTGCACTAAAAAAGCGACTCACATTGTTCCGTTATAATGCCTGAAGTAGATCACAGAATATATCTTCAGGGATCGCATATCTATTAAGTTACTCACTCTTTTCTATTTATGACATGCGCGTGTTTGTATAAATGTAAATGTGAGTCCTTGTTC  >  CP009273/2400110‑2400404

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.