Predicted mutation
evidence seq id position mutation freq annotation gene description
JC REL606 980,293 3 bp x 2 7.6% duplication rpsA → 30S ribosomal protein S1

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 980293 =62 (1.220)5 (0.100) 5/62 NT 7.6% coding (1207/1674 nt) rpsA 30S ribosomal protein S1
?REL606 = 980295 60 (1.040)coding (1209/1674 nt) rpsA 30S ribosomal protein S1

TTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAA  >  REL606/980258‑980328
                                   |                                   
ttCACCTGTCTGACATCTCCTGGAACGTTGCAGGCg                                     >  1:6331964/1‑36 (MQ=255)
ttCACCTGTCTGACATCTCCTGGAACGTTGCAGGCg                                     <  1:5476711/36‑1 (MQ=255)
acctgtctgacatctcctggaacgttgcaggcgaag                                     >  1:889233‑M1/36‑36 (MQ=255)
 tCACCTGTCTGACATCTCCTGGAACGTTGCAGGCga                                    >  1:874383/1‑36 (MQ=255)
  ctgtctgacatctcctggaacgttgcaggcgaagaa                                   >  1:6580088‑M1/34‑36 (MQ=255)
  ctgtctgacatctcctggaacgttgcaggcgaagaa                                   >  1:5313766‑M1/34‑36 (MQ=255)
  cACCTGTCTGACATCTCCTGGAACGTTGCAGGCgaa                                   >  1:5699381/1‑36 (MQ=255)
   tgtctgacatctcatggaacgttgcaggcgaagaag                                  <  1:7348174‑M1/4‑1 (MQ=255)
   aCCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAg                                  >  1:889233‑M2/1‑35 (MQ=255)
    gtctgacatctcctggaacgttgcaggcgaagaaga                                 <  1:86322‑M1/5‑1 (MQ=255)
     cTGTCTGACATCTCCTGGAACGTTGCAGGCGAAgaa                                >  1:6580088‑M2/1‑33 (MQ=255)
     cTGTCTGACATCTCCTGGAACGTTGCAGGCGAAgaa                                >  1:5313766‑M2/1‑33 (MQ=255)
      tGTCTGACATCTCATGGAACGTTGCAGGCGAAgaag                               <  1:7348174‑M2/36‑5 (MQ=255)
       gTCTGACATCTCCTGGAACGTTGCAGGCGAAgaaga                              <  1:86322‑M2/36‑6 (MQ=255)
       gTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGc                              >  1:5675926/1‑36 (MQ=255)
       gTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGc                              <  1:4794591/36‑1 (MQ=255)
        tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                             <  1:389441/36‑1 (MQ=255)
        tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                             >  1:875609/1‑36 (MQ=255)
        tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                             >  1:467556/1‑36 (MQ=255)
        tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                             >  1:4985247/1‑36 (MQ=255)
        tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                             >  1:1442373/1‑36 (MQ=255)
        tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                             <  1:1875374/36‑1 (MQ=255)
         cTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAg                            <  1:2178128/36‑1 (MQ=255)
         cTGACATCTCCGGGAACGGGGCAGGCGAAGAAGCAg                            <  1:1317026/36‑1 (MQ=25)
          tGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGt                           <  1:3730626/36‑1 (MQ=255)
          tGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGt                           <  1:5096602/36‑1 (MQ=255)
           gACATCTCCTGGATCGTTGCAGGCGAAGAAGCAGtt                          >  1:4308600/1‑36 (MQ=255)
           gACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGtt                          <  1:7309796/36‑1 (MQ=255)
           gACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGtt                          >  1:5119138/1‑36 (MQ=255)
           gACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGtt                          <  1:159519/36‑1 (MQ=255)
            aCATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTc                         <  1:573239/36‑1 (MQ=255)
            aCATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTc                         <  1:7461603/36‑1 (MQ=255)
              aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCgc                       >  1:722175/1‑35 (MQ=255)
              aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                       >  1:4617156/1‑36 (MQ=255)
              aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                       >  1:4926456/1‑36 (MQ=255)
              aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                       >  1:6550311/1‑36 (MQ=255)
              aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                       <  1:6934056/36‑1 (MQ=255)
              aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                       <  1:1260336/36‑1 (MQ=255)
               tctcCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTg                      <  1:3854859/36‑1 (MQ=255)
               tctcCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTg                      >  1:1508024/1‑36 (MQ=255)
               tctcCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTg                      <  1:3111086/36‑1 (MQ=255)
                 tcCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGaa                    >  1:1221509/1‑36 (MQ=255)
                 tcCTGGAACGTTGCAGGCGAAGAAGCAGCTCGTGaa                    >  1:816265/1‑36 (MQ=255)
                  ccTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAAt                   >  1:1029167/1‑36 (MQ=255)
                  ccTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAAt                   >  1:3529215/1‑36 (MQ=255)
                  ccTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAAt                   <  1:3506854/36‑1 (MQ=255)
                   cTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATa                  <  1:4984285/36‑1 (MQ=255)
                   cTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATa                  >  1:2244559/1‑36 (MQ=255)
                    tgtaaCGTTGCAGGCGAAGAAGCAGTTCGTGAATAc                 <  1:5276293/33‑1 (MQ=255)
                    tGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATAc                 >  1:115078/1‑36 (MQ=255)
                    aacgttgcaggcgaaGAAGAAGCAGTTCGTGAATAc                 <  1:3199481‑M1/21‑1 (MQ=255)
                     ggAACGTTGCAGCCGAAGCAGCCGTTCGTGAATACa                >  1:3220265/1‑36 (MQ=25)
                      gAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACaa               >  1:2975846/1‑36 (MQ=255)
                      gAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACaa               <  1:6179724/36‑1 (MQ=255)
                      cgttgcaggcgaaGAAGAAGCAGTTCGTGAATACaa               <  1:1958077‑M1/23‑1 (MQ=255)
                       gttgcaggcgaaGAAGAAGCAGTTCGTGAATACaaa              >  1:2409065‑M1/13‑36 (MQ=255)
                       aaCGTTGCAGGCGAAgaagaagcagttcgtgaatac              <  1:3199481‑M2/36‑22 (MQ=255)
                       aaCGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaa              <  1:4333071/36‑1 (MQ=255)
                       aaCGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaa              <  1:6816537/36‑1 (MQ=255)
                       aaCGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaa              >  1:4599821/1‑36 (MQ=255)
                         cGTTGCAGGCGAAgaagaagcagttcgtgaatacaa            <  1:1958077‑M2/36‑24 (MQ=255)
                         cGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaa            >  1:6897058/1‑36 (MQ=255)
                         cGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaa            <  1:46245/36‑1 (MQ=255)
                         cGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaa            >  1:5071238/1‑36 (MQ=255)
                          ggtGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa           >  1:171340/3‑36 (MQ=255)
                          gTTGCAGGCGAAgaagaagcagttcgtgaatacaaa           >  1:2409065‑M2/1‑12 (MQ=255)
                          gTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa           >  1:2438305/1‑36 (MQ=255)
                          gTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa           >  1:3629766/1‑36 (MQ=255)
                          gTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa           >  1:5192115/1‑36 (MQ=255)
                          gTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa           >  1:5191002/1‑36 (MQ=255)
                           ttGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAg          >  1:4635835/1‑36 (MQ=255)
                             gCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGc        >  1:639870/1‑36 (MQ=255)
                              cAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGcg       >  1:3034711/1‑36 (MQ=255)
                               aGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGcga      >  1:2208047/1‑36 (MQ=255)
                               aGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGcga      >  1:4928278/1‑36 (MQ=255)
                               aGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGcga      >  1:490166/1‑36 (MQ=255)
                                ggCGAAGAAGCAGTTCGTGAATACAAAAAAGGcgac     <  1:722881/36‑1 (MQ=255)
                                gaaGAAGAAGCAGTTCGTGAATACAAAAAAGGcgac     <  1:5908452‑M1/33‑1 (MQ=255)
                                 gCGAAGAAGCAGTTCGTGAATACAAAAAAGGcgacg    >  1:19486/1‑36 (MQ=255)
                                  cGAAGAAGCAGTTCGTGAATACAAAAAAGGcgacga   <  1:1662365/36‑1 (MQ=255)
                                  cGAAGAAGCAGTTCGTGAATACAAAAAAGGcgacga   >  1:2823115/1‑36 (MQ=255)
                                   gaagaagaagcagttcgtgaatacaaaaaaggcgac  <  1:5908452‑M2/36‑34 (MQ=255)
                                   GAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGaa  <  1:7313437‑M1/36‑1 (MQ=255)
                                   GAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGaa  <  1:3322760‑M1/36‑1 (MQ=255)
                                   GAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGaa  <  1:1867108‑M1/36‑1 (MQ=255)
                                   |                                   
TTCACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAA  >  REL606/980258‑980328

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG