Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC CP009273 542,271 IS5 (–) +4 bp 35.0% coding (136‑139/1668 nt) fdrA → putative NAD(P)‑binding acyl‑CoA synthetase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 269666 =NA (NA)9 (0.510) 8/276 NT 40.2% noncoding (1195/1195 nt) IS5 repeat region
?CP009273 = 542274 13 (0.750)coding (139/1668 nt) fdrA putative NAD(P)‑binding acyl‑CoA synthetase
* ? CP009273 = 270860NA (NA)5 (0.280) 5/276 NT 27.2% noncoding (320/1221 nt) IS30 repeat region
?CP009273 542271 = 13 (0.750)coding (136/1668 nt) fdrA putative NAD(P)‑binding acyl‑CoA synthetase

GGAGTTGTTATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGG  >  CP009273/542127‑542419
                                                                                                                                                   |                                                                                                                                                 
ggAGTTGTTATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGa                                                                                                                                                  >  1:216297/1‑149 (MQ=255)
    ttgttATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACTAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATAGGTAAGCCCGCCAATAAAGAGTTATTAGATAca                                                                                                                                              >  1:99254/1‑148 (MQ=255)
       ttATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGGGTTATTAGATACcaca                                                                                                                                           >  1:38096/1‑149 (MQ=255)
             tCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCCGCGTATCAGAAAATGTTCATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACGACAGGTTTc                                                                                                                                     >  2:43440/1‑149 (MQ=255)
                tctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcagatcagggacttgttcgcaccttccTTAGATACAACAGGATTCTgg                                                                                                                                  >  1:137129‑M2/129‑149 (MQ=255)
                        attaAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCCGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGggacggtgggaataaacggggagat                                                                                                                          <  1:711‑M2/149‑26 (MQ=255)
                            aaaaaGGGTGTTTTCAGGATTCGGTCAGTTCAATTATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGAtttt                                                                                                                      <  2:135021/149‑1 (MQ=255)
                                        ttCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGAGGATGTTTCCGTAATGATGGGTACGCCCCCCAATAAAGCGTTATGAGATACCACAGGTTTCTGGCATGACGACTTTAATAacgcaacg                                                                                                          >  1:311166/1‑149 (MQ=255)
                                           aGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCTGAAGATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCg                                                                                                       >  1:205929/1‑149 (MQ=255)
                                                    ggcgacggtaaggttatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGATACCACAGGTTTATGGCATGACGATTTTAATAACGCCACGCCGAACGATAtt                                                                                              <  1:67458‑M2/57‑1 (MQ=255)
                                                     cAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcattt                                                                                             <  2:140356‑M2/149‑55 (MQ=255)
                                                             tGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGCTTTTAATACCCCCCCGCCGAACGATATTTGCGTGGCa                                                                                     >  2:66538/1‑149 (MQ=255)
                                                                            gatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAgcg                                                                      <  2:282020‑M2/81‑1 (MQ=255)
                                                                                 aGCGAATCAGAAAATGTGGATGATGTTTCCGTAATGCTGGGTACGCCCGCCAATAAAGCGTTATTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctccc                                                                 >  1:275881‑M2/1‑67 (MQ=255)
                                                                                      accAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCgg                                                            <  2:216297/147‑1 (MQ=255)
                                                                                      aTCAGAAAATGTTGATGATGTTTCCGTACTGATGGGTCCGCCTGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCgg                                                            <  1:249370/149‑1 (MQ=255)
                                                                                                  tgatgaTGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatca                                                >  2:35662‑M2/1‑50 (MQ=255)
                                                                                                  tgatgaTGTTCCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcaggtttgagccgcttttttctcccgtaaatgccttgaatca                                                >  2:166855‑M2/1‑50 (MQ=255)
                                                                                                            tCCGTAATGATGGGTACGCCCGCCTATAAAGCGTTATTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttag                                      <  2:232198‑M2/149‑110 (MQ=255)
                                                                                                               gTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagacc                                   <  2:5761‑M2/149‑113 (MQ=255)
                                                                                                                    gatgGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAGGCGGCGGATGCGGGGGTCGCGCAGGCGATTAAGcagcaggcag                              >  1:186591/1‑145 (MQ=255)
                                                                                                                       gggTACGCCAGCCAGTAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGTCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAag                           >  1:270572/1‑149 (MQ=255)
                                                                                                                          tACGCCCGCCAATAAAGCGTTATTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgcc                        <  2:107277‑M2/149‑124 (MQ=255)
                                                                                                                            cgccaaacaattgggCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGACGCGGCGGACGCGGGGATCGGGCAGGCGATTATGCAGCAGCTTGAAGAGgcg                      >  1:131071/15‑149 (MQ=255)
                                                                                                                              gacttgttagtaccttctTTAGATACCAGCGGTTTCTGGCCTGACAATTGTAGTAACGCCGCGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGAGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCt                    <  2:24320‑M2/131‑1 (MQ=255)
                                                                                                                               ccgccAATAAAGCGTTATTAGACACCACAGGTTTCTGGCATCACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTa                   <  2:99254/149‑1 (MQ=255)
                                                                                                                                  tgttcgcaccttccTTAGACACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTaaaa                <  1:108382‑M2/135‑1 (MQ=255)
                                                                                                                                   cAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAaac               <  1:262012/149‑1 (MQ=255)
                                                                                                                                      tAAAGCGTTATTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgtt            >  2:146218‑M2/1‑14 (MQ=255)
                                                                                                                                                ttaGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAggg  >  2:186956/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGAGTTGTTATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGG  >  CP009273/542127‑542419

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.