Predicted mutation
evidence seq id position mutation freq annotation gene description
JC NC_014034 3,112,237 +GGC 100% coding (382/438 nt) RCAP_rcc02916 ← hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 31122374 (0.020)255 (1.120) 146/192 0.0 0.98 coding (382/438 nt) RCAP_rcc02916 hypothetical protein
?NC_014034 3112238 = 4 (0.020)coding (381/438 nt) RCAP_rcc02916 hypothetical protein

Only 100 of 269 total aligned reads displayed.

GGGTAGGGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAA.......................................................................................................  >  NC_014034/3112138‑3112237
.......................................................................................................GCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATGGCGGCGATC  >  NC_014034/3112238‑3112337
                                                                                                    |||                                                                                                    
GGGTAGGGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAG                                                                                                        >  2:53913/1‑101 (MQ=1)
 GGTAGGGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGG                                                                                                       >  2:2556142/1‑101 (MQ=11)
 GGTAGGGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGG                                                                                                       >  1:3299140/1‑101 (MQ=11)
  GTAGGGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGC                                                                                                      >  1:990392/1‑101 (MQ=31)
    AGGGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGC                                                                                                    <  1:4489226/101‑1 (MQ=32)
    AGGGTTGGACGGACGGCGCGGACGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGC                                                                                                    >  1:2698572/1‑101 (MQ=17)
     GGGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCC                                                                                                   <  1:2094463/101‑1 (MQ=32)
      GGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCT                                                                                                  <  1:2481183/101‑1 (MQ=32)
      GGTTGGACGGACGGCGCCGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCT                                                                                                  <  2:4262384/101‑1 (MQ=17)
       GTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTG                                                                                                 >  2:638683/1‑101 (MQ=32)
        TTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGC                                                                                                >  2:4117238/1‑101 (MQ=32)
          GGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGA                                                                                              >  1:3027838/1‑101 (MQ=32)
              GGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATATGGTCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGC                                                                                          <  1:4054016/101‑1 (MQ=16)
                  GGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGG                                                                                      >  1:4636015/1‑101 (MQ=32)
                  GGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGG                                                                                      >  2:4419439/1‑101 (MQ=32)
                    CGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCT                                                                                    <  1:3914443/101‑1 (MQ=32)
                      CGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCA                                                                                  <  1:3516305/101‑1 (MQ=33)
                      CGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCA                                                                                  >  2:3405275/1‑101 (MQ=32)
                      CGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCA                                                                                  <  1:1208541/101‑1 (MQ=33)
                       GGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAG                                                                                 <  1:766415/101‑1 (MQ=32)
                       GGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAG                                                                                 >  2:238816/1‑101 (MQ=32)
                          CGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGg                                                                              >  2:966887/1‑100 (MQ=17)
                          CGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGT                                                                              <  1:3198670/101‑1 (MQ=32)
                           GTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGgg                                                                             >  1:3234537/1‑99 (MQ=17)
                           GTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTG                                                                             >  1:36611/1‑101 (MQ=32)
                            TGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCATGGTGC                                                                            <  2:2474204/101‑1 (MQ=18)
                            TGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGC                                                                            >  2:3191708/1‑101 (MQ=32)
                            TGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGC                                                                            <  1:2032196/101‑1 (MQ=32)
                            TGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGC                                                                            <  2:4346262/101‑1 (MQ=32)
                              GTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAG                                                                          <  1:225554/101‑1 (MQ=33)
                              GTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAG                                                                          <  2:2302650/101‑1 (MQ=32)
                              GTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCAAAAACGAGAGGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAG                                                                          >  2:3948456/1‑101 (MQ=12)
                                TTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCC                                                                        >  2:1628199/1‑101 (MQ=32)
                                 TTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCG                                                                       >  1:413633/1‑101 (MQ=32)
                                 TTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCG                                                                       >  1:1501958/1‑101 (MQ=32)
                                 TTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCG                                                                       <  1:3105343/101‑1 (MQ=33)
                                  TGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGC                                                                      >  1:2663555/1‑101 (MQ=32)
                                   GCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCG                                                                     <  2:2747500/101‑1 (MQ=32)
                                     TCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGGTCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGA                                                                   >  1:1683464/1‑101 (MQ=17)
                                     TCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGA                                                                   <  1:2400212/101‑1 (MQ=32)
                                     TCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGA                                                                   >  2:1749565/1‑101 (MQ=32)
                                      CTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAA                                                                  >  1:3136408/1‑101 (MQ=32)
                                       TTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAG                                                                 >  1:610036/1‑101 (MQ=32)
                                        TGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGA                                                                <  2:2941502/101‑1 (MQ=32)
                                        TGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGA                                                                <  1:509144/101‑1 (MQ=32)
                                        TGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGA                                                                >  1:4564426/1‑101 (MQ=32)
                                         GCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGAC                                                               >  2:2138907/1‑101 (MQ=32)
                                         GCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGAC                                                               <  1:2838941/101‑1 (MQ=32)
                                          CTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACC                                                              <  1:421936/101‑1 (MQ=32)
                                            CATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCAT                                                            <  1:191130/101‑1 (MQ=32)
                                            CATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCAT                                                            >  1:2229027/1‑101 (MQ=32)
                                             ATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATC                                                           >  2:3478285/1‑101 (MQ=32)
                                             ATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATC                                                           >  1:1515882/1‑101 (MQ=32)
                                             ATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATC                                                           >  2:4563439/1‑101 (MQ=32)
                                              TGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCA                                                          <  2:1198923/101‑1 (MQ=32)
                                               ggcGGCGTGAACCTCAGCTGATGGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATAAC                                                         <  1:1562020/98‑1 (MQ=255)
                                                TGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACG                                                        <  2:1300961/101‑1 (MQ=32)
                                                   GCGTCAATCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTACGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCG                                                     <  2:3142104/101‑1 (MQ=255)
                                                     GTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGAC                                                   <  1:2722441/101‑1 (MQ=32)
                                                     GTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGAC                                                   <  1:2555604/101‑1 (MQ=32)
                                                      TGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACG                                                  <  1:1628299/101‑1 (MQ=32)
                                                      TGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACG                                                  >  2:1816609/1‑101 (MQ=32)
                                                       GAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGC                                                 >  1:1874019/1‑101 (MQ=33)
                                                          CCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGC                                              >  1:478428/1‑101 (MQ=33)
                                                           CTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCC                                             <  2:1538423/101‑1 (MQ=32)
                                                           CTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCCACGCGGCC                                             >  2:4160765/1‑101 (MQ=17)
                                                                CTGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCC                                        <  2:863311/101‑1 (MQ=32)
                                                                 TGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCA                                       >  2:1080012/1‑101 (MQ=33)
                                                                 TGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCA                                       <  2:3808120/101‑1 (MQ=32)
                                                                 TGATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCA                                       <  1:2813073/101‑1 (MQ=32)
                                                                  GATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAG                                      >  1:3603236/1‑101 (MQ=33)
                                                                  GATCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAG                                      >  1:2078516/1‑101 (MQ=33)
                                                                    TCGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGACGACCATCACGCCGACGCGGCCGCCCCAGTC                                    >  2:70667/1‑101 (MQ=17)
                                                                     CGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCC                                   <  1:3769875/101‑1 (MQ=32)
                                                                     CGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCC                                   >  1:581196/1‑101 (MQ=33)
                                                                     CGGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCC                                   >  2:1491757/1‑101 (MQ=32)
                                                                      GGTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCC                                  <  1:1106524/101‑1 (MQ=32)
                                                                       GTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCA                                 <  1:3948698/101‑1 (MQ=32)
                                                                       GTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCA                                 >  2:2225056/1‑101 (MQ=33)
                                                                       GTGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCA                                 <  2:3027650/101‑1 (MQ=32)
                                                                        TGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAG                                >  1:886489/1‑101 (MQ=33)
                                                                        TGATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAG                                >  1:865907/1‑101 (MQ=33)
                                                                          ATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCC                              >  1:1979924/1‑101 (MQ=33)
                                                                          ATCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCC                              <  2:1918307/101‑1 (MQ=32)
                                                                           TCGAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCC                             >  2:2895169/1‑101 (MQ=33)
                                                                             GAAAACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGC                           >  1:639847/1‑101 (MQ=32)
                                                                                 ACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAA                       <  1:2139028/101‑1 (MQ=32)
                                                                                 ACGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCAAGTCCCAGCCCGCCAAA                       >  2:2073588/1‑101 (MQ=17)
                                                                                  CGAGATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAG                      <  1:291976/101‑1 (MQ=32)
                                                                                   GAGATGGCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGC                     <  2:3965709/101‑1 (MQ=17)
                                                                                     GATGTCGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGG                   <  2:4645341/101‑1 (MQ=32)
                                                                                          CGTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATG              >  1:2290717/1‑101 (MQ=32)
                                                                                           GTAGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGA             >  2:38022/1‑101 (MQ=32)
                                                                                             AGACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATG           <  1:4184396/101‑1 (MQ=32)
                                                                                              GACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATGG          >  1:3707265/1‑101 (MQ=32)
                                                                                               ACGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATGGC         <  1:1384396/101‑1 (MQ=32)
                                                                                                CGAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATGGCG        >  2:4213823/1‑101 (MQ=32)
                                                                                                 GAAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATGGCGG       <  2:3603017/101‑1 (MQ=32)
                                                                                                  AAGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATGGCGGC      >  2:3466273/1‑101 (MQ=32)
                                                                                                   AGGCGCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATGGCGGCG     <  1:1816718/101‑1 (MQ=32)
                                                                                                    |||                                                                                                    
GGGTAGGGTTGGACGGACGGCGCGGGCGTGGTTTTGCTCTTGCTCATGTGGGCGTGAACCTCAGCTGATCGGTGATCGAAAACGAGATGTCGTAGACGAA.......................................................................................................  >  NC_014034/3112138‑3112237
.......................................................................................................GCCTGCGACGGCGCGGCTCAGGGTGCAGCCGCGGAAGACCATCACGCCGACGCGGCCGCCCCAGTCCCAGCCCGCCAAAGCGGCGATGATGGCGGCGATC  >  NC_014034/3112238‑3112337

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG