Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,559 G→C 15.6% Q147H (CAG→CGC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,5590G→C15.6% 84.6 ‑1.9 ‑0.1230Q147H (CAG→CAC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/22):  ref base (102/92):  total (116/114)
Fisher's exact test strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.51e-01

Only 100 of 255 total aligned reads displayed.

TGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCG  >  NC_014034/1396459‑1396659
                                                                                                     |                                                                                                     
tGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAg                                                                                                       <  2:3748515/101‑1 (MQ=255)
 ggATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgc                                                                                                      >  2:3730966/1‑101 (MQ=255)
 ggATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgc                                                                                                      <  2:2016854/101‑1 (MQ=255)
 ggATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgc                                                                                                      >  1:2849742/1‑101 (MQ=255)
  gATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgcc                                                                                                     >  2:3863858/1‑101 (MQ=255)
      gAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGGCGCTCCCCCCGCGGGGGCGAGGACCGAAACTGCTTTC.GTGGGGCTGGCTAGccgccc                                                                                                 >  1:3613982/1‑101 (MQ=38)
       aaaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGccgccc.g                                                                                               >  1:520350/1‑101 (MQ=255)
       aaaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGccgccc.g                                                                                               >  2:1173472/1‑101 (MQ=255)
         aCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCc.gcg                                                                                             >  1:2879154/1‑101 (MQ=255)
         aCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCc.gcg                                                                                             >  2:475075/1‑101 (MQ=255)
          cGGCCCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.gcgc                                                                                            <  2:1745843/101‑1 (MQ=255)
          cGGACCTGAAATTCCGCGGGAAACTCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.gcgc                                                                                            <  1:2574129/101‑1 (MQ=255)
           ggACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGcc                                                                                           >  2:93132/1‑101 (MQ=255)
            gACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCa                                                                                          <  1:950825/101‑1 (MQ=255)
             aCCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCaa                                                                                         <  2:3131673/101‑1 (MQ=255)
               cTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGt                                                                                       >  1:2112474/1‑101 (MQ=255)
                  aaaTTCCGCGGGAAACGCGGTGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGc                                                                                    >  1:1104376/1‑101 (MQ=255)
                  aaaTTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGc                                                                                    >  1:2022688/1‑101 (MQ=255)
                   aaTTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCt                                                                                   <  1:475132/101‑1 (MQ=255)
                    aTTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTc                                                                                  <  1:3911975/101‑1 (MQ=255)
                     ttCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCAGCGccgacctcggc                                                                                  <  1:789723/101‑10 (MQ=37)
                     ttCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTcc                                                                                 >  2:363685/1‑101 (MQ=255)
                          cgGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCCCTccccccc                                                                            >  2:3151363/1‑100 (MQ=255)
                           gggAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCcgc                                                                           >  2:3115237/1‑101 (MQ=255)
                            ggAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTGCCCCcgcg                                                                          >  2:3206959/1‑101 (MQ=255)
                             gAAACGCGGCGGCGGGCGCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCTGTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCgg                                                                         <  2:3995428/101‑1 (MQ=38)
                             gAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCgg                                                                         <  2:2924059/101‑1 (MQ=255)
                                aCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCCAGTCGCTCCCCCCGCggggg                                                                      >  1:1679214/1‑101 (MQ=255)
                                 cgcgGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGc                                                                     >  1:758223/1‑101 (MQ=255)
                                    ggcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTg                                                                  >  2:147744/1‑101 (MQ=255)
                                     gcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTgg                                                                 <  1:593468/101‑1 (MQ=255)
                                     gcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTgg                                                                 <  2:2702541/101‑1 (MQ=255)
                                      ttgggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTgg                                                                 >  1:964861/11‑101 (MQ=2)
                                      cggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                <  1:3500837/101‑1 (MQ=255)
                                      cggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                >  1:1921959/1‑101 (MQ=255)
                                        gggacggggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGAt                                                               <  1:3070347/93‑1 (MQ=2)
                                         ggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATc                                                              <  2:1531513/94‑1 (MQ=2)
                                         ggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATc                                                              >  1:710791/8‑101 (MQ=2)
                                          gacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCg                                                             <  2:777493/95‑1 (MQ=2)
                                          gAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                            >  1:3479219/1‑101 (MQ=255)
                                           acggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGa                                                            >  2:1872176/6‑101 (MQ=2)
                                           aGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                           >  2:1671555/1‑101 (MQ=255)
                                            gCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                          <  2:4044722/101‑1 (MQ=255)
                                            gCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCCCCCCCCGCGGGGGCGTGGGTCGaaa                                                          >  1:2124528/1‑101 (MQ=255)
                                             cACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                         >  1:2625086/1‑101 (MQ=255)
                                             cACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                         <  1:3176931/101‑1 (MQ=255)
                                              gcGTCGCGCCCCCCGCGTGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg                                                        <  1:1887697/100‑1 (MQ=255)
                                               agTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAAct                                                        <  1:3098174/100‑2 (MQ=2)
                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt                                                      <  2:4192719/101‑1 (MQ=255)
                                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActg                                                       <  1:2323237/101‑3 (MQ=2)
                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                                    <  2:131230/101‑1 (MQ=255)
                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtg                                                     <  1:4122007/101‑5 (MQ=2)
                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtc                                                   <  1:152770/101‑1 (MQ=255)
                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgc                                                    >  2:164189/1‑96 (MQ=2)
                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgc                                                    <  2:4108099/101‑6 (MQ=2)
                                                     tCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgt                                                 >  2:609393/1‑101 (MQ=255)
                                                     cccccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgt                                                 <  2:3257949/100‑1 (MQ=255)
                                                      ccccccGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccgg                                                 >  1:1183267/1‑93 (MQ=2)
                                                       cccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgtcg                                               >  2:3949492/1‑101 (MQ=255)
                                                       cccccGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccggt                                                <  1:2389571/101‑10 (MQ=2)
                                                       cccccGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccggt                                                >  2:3029435/1‑92 (MQ=2)
                                                        ccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTAGc                                              >  2:3921989/1‑101 (MQ=255)
                                                        ccccGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCC.AACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccggtt                                               <  1:381193/101‑11 (MQ=2)
                                                          ccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCtt                                            >  2:509971/1‑101 (MQ=255)
                                                          ccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCtt                                            <  1:1475602/101‑1 (MQ=255)
                                                           cgcgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTa                                           <  2:1021764/101‑1 (MQ=255)
                                                            gcgGGGGTGTGGATCGAAACTGCTTTC.GTGGGAGTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTac                                          <  1:3665318/101‑1 (MQ=255)
                                                            gcgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTACCCCCGCGGGGGCGAGGATAGAAAAAAATTCGTCGCTTac                                          >  1:445185/1‑101 (MQ=37)
                                                              gggggCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacga                                        >  2:4212282/1‑101 (MQ=255)
                                                               ggggCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacgac                                       >  1:1037435/1‑101 (MQ=255)
                                                                  gCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGcc                                    <  1:1607830/101‑1 (MQ=255)
                                                                    gTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTa                                  >  2:3663005/1‑101 (MQ=255)
                                                                     tGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTAt                                 >  2:2488559/1‑101 (MQ=255)
                                                                      ggATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATc                                >  1:11561/1‑101 (MQ=255)
                                                                      ggATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATc                                >  2:2229406/1‑101 (MQ=255)
                                                                       gATCGAAACTGGTTAC.GTGGGACTGGCCGGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCt                               <  2:3481605/101‑1 (MQ=255)
                                                                        aTCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTg                              <  2:1104282/101‑1 (MQ=255)
                                                                          cGAAACTGCTTTC.GTGGGACTGGCCAGCAGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAt                            >  2:3056537/1‑101 (MQ=255)
                                                                           gAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAtt                           <  1:2164356/101‑1 (MQ=255)
                                                                           gAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAtt                           >  1:4056389/1‑101 (MQ=255)
                                                                           gAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAtt                           <  2:2386246/101‑1 (MQ=255)
                                                                            aaaCTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAttt                          <  1:1966419/101‑1 (MQ=255)
                                                                               cTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCgg                       <  2:2120533/101‑1 (MQ=255)
                                                                                tGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCggg                      >  1:3304975/1‑101 (MQ=255)
                                                                                 gCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGt                     <  1:871722/101‑1 (MQ=255)
                                                                                    ttC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGc                  >  2:2605145/1‑101 (MQ=255)
                                                                                        gTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTcc               >  2:1912821/1‑101 (MQ=255)
                                                                                        gTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTcc               >  1:1368723/1‑101 (MQ=255)
                                                                                        gTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTcc               >  1:174106/1‑101 (MQ=255)
                                                                                         tGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTccc              <  2:3530219/101‑1 (MQ=255)
                                                                                         tGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTccc              <  2:1403651/101‑1 (MQ=255)
                                                                                          gggACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTcccc             >  1:2764017/1‑101 (MQ=255)
                                                                                             gcTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCcg          <  1:2488743/100‑1 (MQ=255)
                                                                                               tGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCcgcg        >  2:3830573/1‑101 (MQ=255)
                                                                                                 gCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCggg      >  1:3209246/1‑101 (MQ=255)
                                                                                                 gCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCggg      <  2:2252992/101‑1 (MQ=255)
                                                                                                  ccAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCgggg     <  2:2542618/101‑1 (MQ=255)
                                                                                                  ccAGCCGCCC.GCGCCAAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGATCCCCCCGCgggg     <  1:3805148/101‑1 (MQ=255)
                                                                                                    aGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGc   >  2:2340562/1‑101 (MQ=255)
                                                                                                     gccgccC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCg  >  1:2803709/1‑101 (MQ=255)
                                                                                                     |                                                                                                     
TGGATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCC.GCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCG  >  NC_014034/1396459‑1396659

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG