Predicted mutation
evidence seq id position mutation annotation gene description
JC NC_000913 3,754,317 +GCG coding (1138/1149 nt) yiaY ← putative oxidoreductase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 37543173 (0.020)129 (0.850)
+GCG
81/144 0.4 97.8% coding (1138/1149 nt) yiaY putative oxidoreductase
?NC_000913 3754318 = 3 (0.020)coding (1137/1149 nt) yiaY putative oxidoreductase

Only 100 of 142 total aligned reads displayed.

GGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTGC‑‑‑ATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTC  >  NC_000913/3754246‑3754393
                                                                           |                                                                           
gggttgcagagatagttgaaaaagcggctagcaatttgccagccgttgtggaaatgatgattacatcgctgctcga                                                                             <  2:687698‑M2/1‑1 (MQ=255)
gggttacagagatagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcga                                                                             >  2:2489785‑M2/76‑76 (MQ=255)
gggttacagagataggtgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcga                                                                             >  2:3862883‑M2/76‑76 (MQ=255)
 ggttacagagatagttgaaaaagcggctaacaatttgccggccgttgtggaaatgatgattacatcgctgcgcgAt                                                                            <  2:2468299‑M2/2‑1 (MQ=255)
 ggttacagagatagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgAt                                                                            >  1:180323‑M2/75‑76 (MQ=255)
   ttacagagatagtggaaaaagcggctaacaatttgccagccgttgtggaaaggatgattacatcgctgcgcgATaa                                                                          <  2:2040934‑M2/4‑1 (MQ=255)
    tacagagatagttgaaaaagaggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATaaa                                                                         <  2:4044724‑M2/5‑1 (MQ=255)
     gcagagatagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAAt                                                                        >  2:2640004‑M2/71‑76 (MQ=255)
      aagagatagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATc                                                                       <  2:5041877‑M2/7‑1 (MQ=255)
        gagatagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGc                                                                     <  1:4954303‑M2/9‑1 (MQ=255)
          gatagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCca                                                                   >  2:3297603‑M2/66‑76 (MQ=255)
           atagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAGATCGCcac                                                                  >  1:793636‑M2/65‑76 (MQ=255)
           atagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCcac                                                                  <  1:2735607‑M2/12‑1 (MQ=255)
            tagttgaaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCcaca                                                                 >  2:4331605‑M2/64‑76 (MQ=255)
                  aaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATttct                                                           <  1:1787796‑M2/19‑1 (MQ=255)
                  aaaaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATttct                                                           >  1:3841905‑M2/58‑76 (MQ=255)
                   aaatgcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATttctt                                                          <  1:1377026‑M2/20‑1 (MQ=255)
                    aaagcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATttcttc                                                         >  1:5357004‑M2/56‑76 (MQ=255)
                    aaagcggctaacaatttgccagccgttgtggaaatgatgattacatccctgcgcgACAAATCGCCACACTttcttc                                                         >  2:5267044‑M2/56‑76 (MQ=255)
                      agcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGt                                                       >  1:182019‑M2/54‑76 (MQ=255)
                      agcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGt                                                       >  2:5435669‑M2/54‑76 (MQ=255)
                       gcggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTg                                                      <  1:5261779‑M2/24‑1 (MQ=255)
                        cggctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGa                                                     >  2:3846156‑M2/52‑76 (MQ=255)
                          gctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGt                                                   >  1:4124125‑M2/50‑76 (MQ=255)
                          gCTAACACTTTGCCAGCCGCTGAGGAAATGATGATTACATCCCTGCGCGATAACACGCCACAATCTCTTCGCGAGt                                                   >  2:5045408/1‑76 (MQ=255)
                           ttaataattggccagctgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGtt                                                  <  1:1441150‑M2/28‑1 (MQ=255)
                           ctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGtt                                                  >  2:3981118‑M2/49‑76 (MQ=255)
                           ctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGtt                                                  <  2:2319395‑M2/28‑1 (MQ=255)
                           ctaacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTAGTGAGtt                                                  >  2:1167178‑M2/49‑76 (MQ=255)
                            taacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTg                                                 >  1:2471867‑M2/48‑76 (MQ=255)
                            taacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTg                                                 <  1:5674117‑M2/29‑1 (MQ=255)
                            taacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCATCACAATTTATTCGTGAGTTg                                                 >  1:3352925‑M2/48‑76 (MQ=255)
                            taacaatttgccagccgttgcggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTg                                                 >  1:3410038‑M2/48‑76 (MQ=255)
                             aacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGc                                                <  2:501846‑M2/30‑1 (MQ=255)
                             aacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGc                                                >  2:2471451‑M2/47‑76 (MQ=255)
                             aacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGc                                                >  2:4702347‑M2/47‑76 (MQ=255)
                             aacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGc                                                >  1:1027566‑M2/47‑76 (MQ=255)
                             aacaatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGc                                                <  2:3161335‑M2/30‑1 (MQ=255)
                             aacaattggccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGc                                                <  2:3430509‑M2/30‑1 (MQ=255)
                               caatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCt                                              <  1:1634522‑M2/32‑1 (MQ=255)
                               caatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCt                                              >  2:4538238‑M2/45‑76 (MQ=255)
                               caatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCt                                              >  2:4389107‑M2/45‑76 (MQ=255)
                                aatttgcgagccgttgtggaaatgatgattccatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTg                                             <  2:300795‑M2/33‑1 (MQ=255)
                                aatttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTg                                             >  1:1028341‑M2/44‑76 (MQ=255)
                                 atttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCCTGAGTTGCCTgg                                            >  2:21029‑M2/43‑76 (MQ=255)
                                  tttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGa                                           >  1:3288634‑M2/42‑76 (MQ=255)
                                   ttgccggccgttgtggaaatgatgtttacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGATTTGCCTGGAt                                          <  2:3127107‑M2/36‑1 (MQ=255)
                                   ttgccagccgttgtgggaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGGTGCCTGGAt                                          >  1:3229996‑M2/41‑76 (MQ=255)
                                   ttgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGAt                                          >  2:423216‑M2/41‑76 (MQ=255)
                                    tgccagccgttgtggaaatgatgattacatcgctgcgctATAAATCGCCCCAATTTCTTCGTGAGTTGcctggctc                                         <  1:3493101‑M2/37‑4 (MQ=255)
                                    tgccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATc                                         <  1:1547176‑M2/37‑1 (MQ=255)
                                     gccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCg                                        >  1:4467334‑M2/39‑76 (MQ=255)
                                     gccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCg                                        >  2:4167349‑M2/39‑76 (MQ=255)
                                     gccagcaggtgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCg                                        >  1:4308636‑M2/39‑76 (MQ=255)
                                      ccagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCgg                                       <  1:2182110‑M2/39‑1 (MQ=255)
                                       cagccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCggg                                      >  1:2072872‑M2/37‑76 (MQ=255)
                                        agccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGt                                     <  1:3663461‑M2/41‑1 (MQ=255)
                                         gccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGtt                                    <  2:5732661‑M2/42‑1 (MQ=255)
                                          ccgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTa                                   <  2:2909241‑M2/43‑1 (MQ=255)
                                           cgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAg                                  >  1:4054276‑M2/33‑76 (MQ=255)
                                           cgttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAg                                  >  2:1113508‑M2/33‑76 (MQ=255)
                                            gttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGt                                 <  2:4556994‑M2/45‑1 (MQ=255)
                                             ttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTa                                <  1:3862884‑M2/46‑1 (MQ=255)
                                             ttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTa                                <  1:2289397‑M2/46‑1 (MQ=255)
                                             ttgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTa                                >  1:1100374‑M2/31‑76 (MQ=255)
                                              tgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTaa                               <  1:4040257‑M2/47‑1 (MQ=255)
                                              tgtggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTaa                               >  2:2321625‑M2/30‑76 (MQ=255)
                                               ggggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTaaa                              <  2:946276‑M2/48‑1 (MQ=255)
                                                tggaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAg                             >  1:5836505‑M2/28‑76 (MQ=255)
                                                  gaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGcc                           >  2:5788310‑M2/26‑76 (MQ=255)
                                                  gaaatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGcc                           >  1:4719138‑M2/26‑76 (MQ=255)
                                                    aatgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCcac                         <  1:4826691‑M2/53‑1 (MQ=255)
                                                     atgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCcaca                        <  1:5714162‑M2/54‑1 (MQ=255)
                                                      tgatgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAg                       <  1:3529227‑M2/55‑1 (MQ=255)
                                                        atgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGc                     <  2:3835667‑M2/57‑1 (MQ=255)
                                                         tgattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCa                    >  2:263087‑M2/19‑76 (MQ=255)
                                                          gattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCAt                   <  2:149556‑M2/59‑1 (MQ=255)
                                                          gattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCAt                   >  2:1369390‑M2/18‑76 (MQ=255)
                                                          gattacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCAt                   >  2:5629008‑M2/18‑76 (MQ=255)
                                                           attacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATc                  >  2:5238335‑M2/17‑76 (MQ=255)
                                                           attacatcgctgcgcgATAAATCGCCACAATTTCCTCGTGAGGTGCCTGGATCGGGTTAGTAAAGCCACAGGCATc                  >  2:2014749‑M2/17‑76 (MQ=255)
                                                            ttacatcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCt                 <  1:3037783‑M2/61‑1 (MQ=255)
                                                               catcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTc              <  2:2263279‑M2/64‑1 (MQ=255)
                                                               catcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTc              >  1:3072098‑M2/13‑76 (MQ=255)
                                                                atcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCgg                                       <  1:5150402‑M2/39‑1 (MQ=255)
                                                                 tcgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAg            <  2:182019‑M2/66‑1 (MQ=255)
                                                                  cgctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAgg           >  2:2165750‑M2/10‑76 (MQ=255)
                                                                   gctgcgcgATAAATCGCCAGAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAggg          <  1:4280493‑M2/68‑1 (MQ=255)
                                                                   gctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAggg          >  1:798010‑M2/9‑76 (MQ=255)
                                                                   gctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAggg          >  2:4509867‑M2/9‑76 (MQ=255)
                                                                   gctgcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAggg          <  2:5192411‑M2/68‑1 (MQ=255)
                                                                      gcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCAt       <  1:3736159‑M2/71‑1 (MQ=255)
                                                                      gcgcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCAt       >  2:463841‑M2/6‑76 (MQ=255)
                                                                       c‑‑‑ATAATTCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAACGACACAGTCATCTTTCAGCGCACTAg  >  1:2785827/1‑75 (MQ=255)
                                                                        gcgATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTa     <  2:2777726‑M2/73‑1 (MQ=255)
                                                                          gATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGt   >  2:3283768‑M2/2‑76 (MQ=255)
                                                                          gATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGt   >  1:2921275‑M2/2‑76 (MQ=255)
                                                                           aTAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAATCCACAGTCATCTTTCAGGGCATTAGTc  <  1:2673996/76‑1 (MQ=255)
                                                                           aTAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTc  <  2:3405732/76‑1 (MQ=255)
                                                                           aTAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTc  <  1:2511696/76‑1 (MQ=255)
                                                                           |                                                                           
GGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTGC‑‑‑ATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTC  >  NC_000913/3754246‑3754393

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