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

CACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAAT  >  REL606/980260‑980330
                                   |                                   
ctgtctgacatctcctggaacgttgcaggcgaagaa                                     >  1:5313766‑M1/34‑36 (MQ=255)
ctgtctgacatctcctggaacgttgcaggcgaagaa                                     >  1:6580088‑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:5313766‑M2/1‑33 (MQ=255)
   cTGTCTGACATCTCCTGGAACGTTGCAGGCGAAgaa                                  >  1:6580088‑M2/1‑33 (MQ=255)
    tGTCTGACATCTCATGGAACGTTGCAGGCGAAgaag                                 <  1:7348174‑M2/36‑5 (MQ=255)
     gTCTGACATCTCCTGGAACGTTGCAGGCGAAgaaga                                <  1:86322‑M2/36‑6 (MQ=255)
     gTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGc                                <  1:4794591/36‑1 (MQ=255)
     gTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGc                                >  1:5675926/1‑36 (MQ=255)
      tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                               >  1:4985247/1‑36 (MQ=255)
      tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                               >  1:875609/1‑36 (MQ=255)
      tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                               >  1:467556/1‑36 (MQ=255)
      tCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCa                               <  1:389441/36‑1 (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:4926456/1‑36 (MQ=255)
            aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                         >  1:4617156/1‑36 (MQ=255)
            aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                         <  1:1260336/36‑1 (MQ=255)
            aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                         >  1:6550311/1‑36 (MQ=255)
            aTCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGt                         <  1:6934056/36‑1 (MQ=255)
             tctcCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTg                        <  1:3111086/36‑1 (MQ=255)
             tctcCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTg                        <  1:3854859/36‑1 (MQ=255)
             tctcCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTg                        >  1:1508024/1‑36 (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:2244559/1‑36 (MQ=255)
                 cTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATa                    <  1:4984285/36‑1 (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:6179724/36‑1 (MQ=255)
                    gAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACaa                 >  1:2975846/1‑36 (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:6816537/36‑1 (MQ=255)
                     aaCGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaa                >  1:4599821/1‑36 (MQ=255)
                     aaCGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaa                <  1:4333071/36‑1 (MQ=255)
                       cGTTGCAGGCGAAgaagaagcagttcgtgaatacaa              <  1:1958077‑M2/36‑24 (MQ=255)
                       cGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaa              >  1:5071238/1‑36 (MQ=255)
                       cGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaa              <  1:46245/36‑1 (MQ=255)
                       cGTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaa              >  1:6897058/1‑36 (MQ=255)
                        ggtGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa             >  1:171340/3‑36 (MQ=255)
                        gTTGCAGGCGAAgaagaagcagttcgtgaatacaaa             >  1:2409065‑M2/1‑12 (MQ=255)
                        gTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa             >  1:3629766/1‑36 (MQ=255)
                        gTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa             >  1:2438305/1‑36 (MQ=255)
                        gTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa             >  1:5191002/1‑36 (MQ=255)
                        gTTGCAGGCGAAGAAGCAGTTCGTGAATACaaaaaa             >  1:5192115/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:2823115/1‑36 (MQ=255)
                                cGAAGAAGCAGTTCGTGAATACAAAAAAGGcgacga     <  1:1662365/36‑1 (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)
                                  AAGAAGCAGTTCGTGAATACAAAAAAGGCGACGaaa   <  1:7001348‑M1/36‑1 (MQ=255)
                                  AAGAAGCAGTTCGTGAATACAAAAAAGGCGACGaaa   >  1:4863660‑M1/1‑36 (MQ=255)
                                  AAGAAGCAGTTCGTGAATACAAAAAAGGCGACGaaa   <  1:136648‑M1/36‑1 (MQ=255)
                                   AGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAAt  <  1:3263538‑M1/36‑1 (MQ=255)
                                   |                                   
CACCTGTCTGACATCTCCTGGAACGTTGCAGGCGAAGAAGCAGTTCGTGAATACAAAAAAGGCGACGAAAT  >  REL606/980260‑980330

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