Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 2,400,257 IS5 (+) +4 bp 87.2% 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)18 (1.250) 15/274 NT 85.3% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 2400257 = 3 (0.210)intergenic (‑137/‑783) lrhA/alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2
* ? CP009273 = 2283592NA (NA)23 (1.590) 23/274 NT 88.2% noncoding (1/1195 nt) IS5 repeat region
?CP009273 = 2400260 3 (0.210)intergenic (‑140/‑780) lrhA/alaA transcriptional repressor of flagellar, motility and chemotaxis genes/valine‑pyruvate aminotransferase 2

CCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCC  >  CP009273/2282249‑2282536
                                                                                                                                                     |                                                                                                                                          
ccAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCCACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATaag                                                                                                                                            >  1:180262/1‑150 (MQ=255)
              atGCTGTTCTCCGCCAGGATTGATTATCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATaa                                                                                                                              >  1:195805/1‑150 (MQ=255)
                   gctctcCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAATGTGCGAATAAGCggg                                                                                                                         <  1:239090/148‑1 (MQ=255)
                                                  ggCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat                                                                                              <  2:200144/146‑1 (MQ=255)
                                                     ggTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGGCTCGGtcatttc                                                                                           >  2:62047/1‑146 (MQ=255)
                                                     ggTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAAAACCATCAACTGTCTCACCATCCTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGACTAAGCGGGGAAATTTTTCTCGGATAGCTCAGtcacttt                                                                                           >  2:96459/1‑145 (MQ=255)
                                                            tAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCTAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCGGACTCAGTCATTTCATTTCTTCATg                                                                                >  1:162019/1‑150 (MQ=255)
                                                                           tCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAttttt                                                                 >  1:121221/1‑150 (MQ=255)
                                                                                                            ccATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTAACTCCGGCAGTTCAGTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAg                                <  1:127825/150‑1 (MQ=18)
                                                                                                                             tGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTGTCTTCGc                   >  2:163177/1‑146 (MQ=255)
                                                                                                                                          ggTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcn  <  1:242010/150‑2 (MQ=25)
                                                                                                                                                     |                                                                                                                                          
CCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCC  >  CP009273/2282249‑2282536

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