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

TCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTT  >  CP009273/2282270‑2282532
                                                                                                                                |                                                                                                                                      
tctcCGCCAGGATTGCTTTTCGTATTCTTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACACAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATaa                                                                                                                          <  1:363812/143‑1 (MQ=255)
      ccAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGTATATAAGGAAGGTGCGAATAAGCGTGGAAAttc                                                                                                              <  2:465631/149‑1 (MQ=255)
        aGGATTGATTTTCGTATTATTGGCGGTCCCATAACCCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATTCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                            <  2:587018/149‑1 (MQ=255)
           attgattTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAAtt                                                                                                               <  1:107607/143‑1 (MQ=255)
                ttttCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTg                                                                                                    >  2:16852/1‑149 (MQ=255)
                   tCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACt                                                                                                 >  2:315650/1‑149 (MQ=255)
                     gTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa                                                                                               >  2:579929/1‑149 (MQ=255)
                         tattGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGCTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtca                                                                                           <  2:567761/149‑1 (MQ=255)
                              gCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtc                                                                                            <  1:224461/143‑1 (MQ=255)
                                     cATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTCTCTCACCATACTCAACAGTTGATGCCAATCCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                                                                               <  2:448505/149‑1 (MQ=255)
                                                                  acaacaCCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaaa                                                  <  2:23607/149‑1 (MQ=255)
                                                                     acacCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGt                                                     >  1:285808/1‑143 (MQ=255)
                                                                        ccATCAACTTGCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTGTTCTCCCGTAAATGCCtt                                            <  2:467940/149‑1 (MQ=255)
                                                                           tCAACTGTCTCACCATCCTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGaa                                         <  2:363376/149‑1 (MQ=255)
                                                                             aaCTGTCTCACCATACTCAACAGTTGATGCCAATGCCATGGCTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTGCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                             >  1:235832/1‑143 (MQ=255)
                                                                              aCTGGCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                      <  2:226687/149‑1 (MQ=255)
                                                                                    tcACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                      >  1:33919/1‑143 (MQ=255)
                                                                                       ccATACTCAACAGTAGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCCTGAATCAGCCTATTTa                             <  2:399308/149‑1 (MQ=255)
                                                                                           aCTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                               <  1:523443/143‑1 (MQ=255)
                                                                                               aaccacttttaccAATTCTATGGGAGAATATAATGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGttt                     <  2:460772/138‑1 (MQ=14)
                                                                                               aaCAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGttt                     >  2:527076/1‑149 (MQ=255)
                                                                                                            aaTGCTATGGGTGAATATAAGGAAGGGGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTaa        <  2:324519/149‑1 (MQ=255)
                                                                                                                cTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt          >  1:171951/1‑143 (MQ=40)
                                                                                                                cTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTATTCATGTTTGATACGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt          >  1:202885/1‑143 (MQ=255)
                                                                                                                        gAATATAAGGAAGGTGCGCAGATGGGGTGAGATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAGTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGtt  >  1:433447/1‑143 (MQ=11)
                                                                                                                                |                                                                                                                                      
TCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTT  >  CP009273/2282270‑2282532

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