Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 3,015,771 IS150 (–) +4 bp coding (1192‑1195/1677 nt) ECB_02816 → KpsD protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 2775877NA (NA)49 (0.940) 33/66 0.2 98.0% noncoding (1/1443 nt) IS150 repeat region
?REL606 3015771 = 1 (0.200)coding (1192/1677 nt) ECB_02816 KpsD protein
* ? REL606 = 30157741 (0.150)55 (1.000) 32/70 0.3 98.2% coding (1195/1677 nt) ECB_02816 KpsD protein
?REL606 = 3894996 NA (NA)noncoding (1443/1443 nt) IS150 repeat region

Only 100 of 117 total aligned reads displayed.

TGAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGC  >  REL606/3015736‑3015806
                                   |                                   
 gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACt                                    <  1:478905/36‑1 (MQ=255)
 gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACt                                    >  1:1723748/1‑36 (MQ=255)
 gAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACt                                    <  2:52450/36‑1 (MQ=255)
  catttcatcgtccaacaaaatgggcgcagtacaCtt                                   >  1:3707352‑M2/34‑36 (MQ=255)
  aaTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACtt                                   <  2:2140748/36‑1 (MQ=255)
  aaTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACtt                                   <  1:872366/36‑1 (MQ=255)
  aaTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACtt                                   >  1:1318575/1‑36 (MQ=255)
  aaTCCCTCACTGCAAAAGATGGAGAAAGCCTCACtt                                   <  1:2254253/36‑1 (MQ=19)
   atttcatcgtccaacaaaatgggtgcagtacaccgg                                  >  2:3225739‑M2/33‑33 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacaccgg                                  <  2:2431258‑M2/4‑4 (MQ=255)
   atttcatcgtccaacaaaatgggtgcagtacacaac                                  >  2:1686687‑M2/33‑33 (MQ=255)
   aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACttt                                  >  2:2378829/1‑36 (MQ=255)
   aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACttt                                  <  1:188746/36‑1 (MQ=255)
   aTCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACttt                                  >  1:3297775/1‑36 (MQ=255)
     ttcatcgtccagcaaaatgggtgcagtacaCTTTCt                                <  1:3007280‑M2/6‑1 (MQ=255)
     ctctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtg                                >  2:328517‑M1/1‑34 (MQ=255)
      tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                               >  2:1469917‑M1/1‑33 (MQ=255)
      tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                               >  2:372091‑M1/1‑33 (MQ=255)
      tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                               <  1:1140254‑M1/36‑4 (MQ=255)
      tctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgt                               <  2:1300651‑M1/36‑4 (MQ=255)
      tcatcgtccgacaaaatgggtgcagtacaCTTTCTg                               >  1:1294855‑M2/30‑36 (MQ=255)
       ctcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgta                              <  2:2493685‑M1/36‑5 (MQ=255)
       catcgtccaacaaaatgggtgcagtacaCTTTCTGc                              >  1:791289‑M2/29‑36 (MQ=255)
        tcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtac                             >  1:1176597‑M1/1‑31 (MQ=255)
        tcGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtac                             >  1:1278072‑M1/1‑31 (MQ=255)
        atcgtccaacaaaatgggtgcagtacaCTTTCTGcc                             >  2:3663076‑M2/28‑36 (MQ=255)
         tcgtccaacaaaatgggtgcagtacaCTTTCTGccc                            >  1:2109187‑M2/27‑36 (MQ=255)
         cGCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtact                            <  1:3464218‑M1/36‑7 (MQ=255)
          gCTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactg                           >  2:127834‑M1/1‑29 (MQ=255)
          gCTGCAAAAGCTGGAGGAAGCCTCACTCTtgtactg                           >  1:3234402‑M1/1‑29 (MQ=255)
           gtccaacaaaatgggtgccgtacaCTTTCTGCCCAg                          >  2:2007225‑M2/25‑36 (MQ=255)
           gtccaacaaaatgggtgcagtacaCTTTCTGCCCAg                          <  2:1626024‑M2/12‑1 (MQ=255)
           cTGCAAAAGCTGGAGGAAGCCTCACTTTtgtactga                          >  1:2794001‑M1/1‑28 (MQ=255)
           cTGCAAAAGATGGAGGAAGCCTCACTTTtgtactga                          >  1:1716309‑M1/1‑28 (MQ=255)
            tccaacaaaatgggtgcagtacaCTTTCTGCCCAGt                         <  1:312984‑M2/13‑1 (MQ=255)
            tccaacaaaatgggtgcagtacaCTTTCTGCCCAGt                         >  1:3464077‑M2/24‑36 (MQ=255)
             gCAAAAGCTGGAGGAAGCCTCACTTTtgtactgacc                        <  1:676952‑M1/36‑11 (MQ=255)
             ccaacaaaatgggtgcagtacaCTTTCTGCCCAGTc                        <  2:1109873‑M2/14‑1 (MQ=255)
              caacaaaatgggtgcagtacaCTTTCTGCCCAGTCa                       >  2:3248431‑M2/22‑36 (MQ=255)
              caacaaaatgggtgcagtacaCTTTCTGCCCAGTCa                       >  2:2309774‑M2/22‑36 (MQ=255)
               aagagctGGAGGAAGCCTCACTTTtgtactgacccc                      <  1:652336‑M1/33‑13 (MQ=255)
               aaaaGCTGGATGAAGCCTCACTTTtgtactgacccc                      <  1:681680‑M1/36‑13 (MQ=255)
               aaaaGCTGGAGGAAGCCTCACTTTtgtactgacccc                      <  1:2299269‑M1/36‑13 (MQ=255)
                acaaaattggtgcagtacaCTTTCTGCCCAGTCATc                     <  2:173679‑M2/17‑1 (MQ=255)
                aaaGCTGGAGGAAGCCTCACTTTtgtactgacccca                     <  1:2612123‑M1/36‑14 (MQ=255)
                aaaGCTGGAGGAAGCCTCACTTTtgtactgacccca                     >  2:1959424‑M1/1‑23 (MQ=255)
                 aaGCTGGAGGAAGCCTCACTTTtgtactgaccccaa                    >  2:1567932‑M1/1‑22 (MQ=255)
                  aaaatgggtgcagtacaCTTTCTGCCCAGTTATcca                   >  1:3118564‑M2/18‑36 (MQ=255)
                  aaaatgggtgcagtacaCTTTCTGCCCAGTCATcca                   <  1:2841118‑M2/19‑1 (MQ=255)
                  aaaatgggtgcagtacaCTTTCTGCCCAGTCATcca                   >  1:2688366‑M2/18‑36 (MQ=255)
                  aaaatgggtgcagtacaCTTTCTGCCCAGTCATcca                   >  1:26606‑M2/18‑36 (MQ=255)
                  aGCTGGAGGAAGCCTCACTTTagtactgaccccaaa                   >  1:3570011‑M1/1‑21 (MQ=255)
                   gCTGGAGGAAGCCTCACTTTtgtactgaccccaaaa                  >  2:2477106‑M1/1‑20 (MQ=255)
                   gCTGGAGGAAGCCTCACTTTtgtactgaccccaaaa                  >  1:2365250‑M1/1‑20 (MQ=255)
                   aaatgggtgcagtacaCTTTCTGCCCAGTCATccac                  <  1:1787012‑M2/20‑1 (MQ=255)
                   aaatgggtgcagtacaCTTTCTGCCCAGTCATccac                  >  1:2235403‑M2/17‑36 (MQ=255)
                   aaaggggtgcagtacaCTTTCTGCCCAGTCATccac                  <  2:2072391‑M2/20‑1 (MQ=255)
                    aatgggcgcagtacaCTTTCTGCCCAGTCATccacc                 >  2:2628651‑M2/16‑36 (MQ=255)
                     tGGAGGAAGCCTCACTTTtgtactgaccccaaaaag                >  1:3094531‑M1/1‑18 (MQ=255)
                      tgggtgcagtaccCTTTCTGCCCAGTCATCCACCaa               >  2:1784357‑M2/14‑36 (MQ=255)
                      tgggtgcagtacaCTTTCTGCCCAGTCATCCACCaa               <  2:430658‑M2/23‑1 (MQ=255)
                      ggaggaAGCCTCACTTTtgtactgaccccaaaaagt               >  1:3520078‑M1/1‑17 (MQ=255)
                      ggaggaAGCCTCACTTTtgtactgaccccaaaaagt               >  2:3431683‑M1/1‑17 (MQ=255)
                       gggtgcagtacaCTTTCTGCCCAGTCATCCACCaaa              <  1:892269‑M2/24‑1 (MQ=255)
                        ggtgcagtacaCTTTCTGCCCAGTCATCCAccaacg             >  1:455081‑M2/12‑34 (MQ=255)
                         gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga            >  2:3676929‑M2/11‑36 (MQ=255)
                         gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga            >  1:1122133‑M2/11‑36 (MQ=255)
                         gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga            >  1:289823‑M2/11‑36 (MQ=255)
                         gtgcagtacaCTTTCTGCCCAGTCATCCACCAaaga            >  1:764288‑M2/11‑36 (MQ=255)
                         ggaAGCCTCACTTTtgtactgaccccaaaaagttgg            >  1:2976673‑M1/1‑14 (MQ=255)
                         ggaAGCCTCACTTTtgtactgaccccaaaaagttgg            <  1:3223186‑M1/36‑23 (MQ=255)
                          gaAGCCTCACTTTtgtactgaccccaaaacgttgga           >  2:756664‑M1/1‑13 (MQ=255)
                          gaAGCCTCACTTTtgtactgaccccaaaaagttgga           >  2:3699850‑M1/1‑13 (MQ=255)
                          gaAGCCTCACTTTtgtactgaccccaaaaagttgga           <  2:2550722‑M1/36‑24 (MQ=255)
                           gcagtacaCTTTCTGCCCAGTCATCCACCAaagaag          <  2:2675430‑M2/28‑1 (MQ=255)
                           aaGCCTCACTTTtgtactgaccccaaaaggttggac          <  2:1616879‑M1/36‑25 (MQ=255)
                            cagtacaCTTTCTGCCCAGTCATCCACCAaagaaga         <  1:91716‑M2/29‑1 (MQ=255)
                            aGCCTCACTTTtgtactgaccccaacaagttggaca         >  1:848799‑M1/1‑11 (MQ=255)
                            aGCCTCACTTTtgtactgaccccaaaaagttggaca         >  2:1372317‑M1/1‑11 (MQ=255)
                             agtacaCTTTCTGCCCAGTCATCCACCAaagaagaa        <  1:3279116‑M2/30‑1 (MQ=255)
                               tacaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGc      >  1:3111387‑M2/5‑36 (MQ=255)
                               tacaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGc      >  1:1779211‑M2/5‑36 (MQ=255)
                               tacaCTTGCTGCCCAGTCATCCACCAAAGAAGAAGc      <  2:3561696‑M2/32‑1 (MQ=255)
                               cTCACTTTtgtactgaccccaaaaagttggacagtt      <  2:1397468‑M1/36‑29 (MQ=255)
                               cTCACTTTtgtactgaccccaaaaagttggacagtt      >  1:2861377‑M1/1‑8 (MQ=255)
                                tCACTTTtgtactgaccccaaaaagttggacagtta     >  1:2435280‑M1/1‑7 (MQ=255)
                                acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc     >  1:3060729‑M2/4‑36 (MQ=255)
                                acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc     <  1:1645864‑M2/33‑1 (MQ=255)
                                acaCTTTCTGCCCAGTCATCCACCAAAGAAGAAGcc     <  1:2628850‑M2/33‑1 (MQ=255)
                                acACTATCAGCCCAGACATCCACCAAAGAAGAagcc     <  1:3209652/35‑1 (MQ=25)
                                 caCTTTCTGCCCAGTCATCCACCAAAGACGAAGCCa    >  2:796419‑M2/3‑36 (MQ=255)
                                 cACTTTtgtactgaccccaaaaagttggacagttaa    >  1:577402‑M1/1‑6 (MQ=255)
                                  aCTTTtgtactgaccccaaaaagttggacagttaaa   >  2:1268307‑M1/1‑5 (MQ=255)
                                  aCTTTtgtactgaccccaaaaagttggacagttaaa   <  2:289993‑M1/36‑32 (MQ=255)
                                  aCTTTtgtactgaccccaaaaagttggacagttaaa   >  1:396217‑M1/1‑5 (MQ=255)
                                  aCTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAg   >  2:2734471‑M2/2‑36 (MQ=255)
                                   cTTTtgtactgaccccaaaaagttggacagttaaac  <  1:2692004‑M1/36‑33 (MQ=255)
                                   cTTTCTGCCCAGTCATCCACCAAAGAAGAagccagc  >  2:922292/1‑36 (MQ=255)
                                   cTTTCTGCCCAGTCATCCACCAAAGAAGAagccagc  <  1:2673107/36‑1 (MQ=255)
                                   cTTTCTGCCCAGTCATCCACCAAAGAAGAagccagc  <  2:2262820/36‑1 (MQ=255)
                                   |                                   
TGAATCTCTCGCTGCAAAAGCTGGAGGAAGCCTCACTTTCTGCCCAGTCATCCACCAAAGAAGAAGCCAGC  >  REL606/3015736‑3015806

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG