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

TTAAGGAGTTGTTATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGG  >  CP009273/542123‑542419
                                                                                                                                                    |                                                                                                                                                    
ttAAGGAGTGGTTATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGttat                                                                                                                                                      <  1:302823/149‑1 (MQ=255)
 tAAGGAGTTGTTATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGttatt                                                                                                                                                     <  1:292792/149‑1 (MQ=255)
    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)
                                                                                                                                                    |                                                                                                                                                    
TTAAGGAGTTGTTATGATCCACGCCTTTATTAAAAAAGGGTGTTTTCAGGATTCGGTCAGTTTAATGATTATTTCACGAAAACTCAGCGAATCAGAAAATGTTGATGATGTTTCCGTAATGATGGGTACGCCCGCCAATAAAGCGTTATTAGATACCACAGGTTTCTGGCATGACGATTTTAATAACGCCACGCCGAACGATATTTGCGTGGCAATTCGTAGCGAAGCGGCGGATGCGGGGATCGCGCAGGCGATTATGCAGCAGCTTGAAGAGGCGCTAAAACAACTGGCGCAGGG  >  CP009273/542123‑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.