Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,940 A→G 4.1% H239H (CAT→CAC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9400A→G 4.1% 9.3 ‑0.7 ‑0.2192H239H (CAT→CACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/3):  ref base (87/97):  total (92/100)
Fisher's exact test strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.65e-01

Only 100 of 216 total aligned reads displayed.

CACGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGA  >  NC_014034/221841‑222039
                                                                                                   |                                                                                                   
     tCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTg                                                                                               <  1:4331900/101‑1 (MQ=255)
      cGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTgg                                                                                              <  1:2914902/101‑1 (MQ=255)
        aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGc                                                                                            <  2:4380411/101‑1 (MQ=255)
        aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGc                                                                                            <  2:1411292/101‑1 (MQ=255)
          ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                          <  1:2800200/101‑1 (MQ=255)
          ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                          <  2:3022950/101‑1 (MQ=255)
          ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGGGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGCCGGATGCGGGGGCtc                                                                                          >  1:1634603/1‑101 (MQ=255)
           cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcg                                                                                         <  1:508946/101‑1 (MQ=255)
           cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcg                                                                                         <  1:2161191/101‑1 (MQ=255)
            ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                        >  2:3382063/1‑101 (MQ=255)
             gTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggt                                                                                       <  2:1021567/101‑1 (MQ=255)
                tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtcgg                                                                                    <  2:3868947/101‑1 (MQ=255)
                tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCAACGACGGTGGGCGTCGGATGCGTGGGCtcggtcgg                                                                                    <  1:4098587/101‑1 (MQ=255)
                 cGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCgga                                                                                   <  2:422840/101‑1 (MQ=255)
                 cGCCGGTGCCCGACACCGTCGCCGAGGCTGCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCgga                                                                                   <  1:3169516/101‑1 (MQ=255)
                  gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                  <  2:3095386/101‑1 (MQ=255)
                  gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGCGTTGGGGGTCAGCGTCGCCGGATCCCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCtcggtcggtc                                                                                  >  1:3734544/1‑99 (MQ=255)
                   ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                 >  2:3399520/1‑101 (MQ=255)
                   ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                 <  2:2418473/101‑1 (MQ=255)
                     ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                               <  1:3406758/101‑1 (MQ=255)
                         cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgg                                                                           >  1:2323019/1‑101 (MQ=255)
                            gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggcc                                                                        >  1:2047417/1‑99 (MQ=255)
                             acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                       >  2:1717864/1‑101 (MQ=255)
                              cacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccg                                                                      <  2:3812482/101‑1 (MQ=255)
                               acCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgg                                                                     >  2:3571311/1‑100 (MQ=255)
                                 cGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgg                                                                   >  2:4114083/1‑101 (MQ=255)
                                   tCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggga                                                                 <  1:478697/101‑1 (MQ=255)
                                   gcGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggga                                                                 <  2:885941/100‑1 (MQ=255)
                                     gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                               <  2:4041019/101‑1 (MQ=255)
                                      ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                              <  2:2391339/101‑1 (MQ=255)
                                       cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcgg                                                             >  2:860885/1‑101 (MQ=255)
                                        gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                            >  2:1160995/1‑101 (MQ=255)
                                          ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcg                                                          <  2:1576670/101‑1 (MQ=255)
                                           gCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgg                                                         >  1:3469868/1‑101 (MQ=255)
                                             ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcggtgggatcggtcgggt                                                       >  1:756821/1‑101 (MQ=255)
                                              tCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtc                                                      >  1:221183/1‑101 (MQ=255)
                                               cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGGGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                     >  1:1344715/1‑101 (MQ=255)
                                                ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                    >  1:4228742/1‑101 (MQ=255)
                                                ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcggtcgggcccg                                                    >  1:1323444/1‑101 (MQ=255)
                                                 gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                   >  2:3614997/1‑101 (MQ=255)
                                                  tgtTGGTGGTCAGCGTCGCCTGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGGCggatcggtcgggtccgtgggatcggtcggatccgtg                                                  <  2:146380/101‑1 (MQ=255)
                                                  tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                  >  2:3148172/1‑101 (MQ=255)
                                                  tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                  <  1:3684990/101‑1 (MQ=255)
                                                   gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                 <  1:4236293/101‑1 (MQ=255)
                                                   gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                 >  1:2435846/1‑101 (MQ=255)
                                                    ttGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                <  2:1816860/101‑1 (MQ=255)
                                                    ttGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                <  1:2250892/101‑1 (MQ=255)
                                                     tggtggTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAtcggtcgggtccgtgggatcggtcgggtccgtggga                                               <  1:534991/101‑1 (MQ=255)
                                                      ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                              <  1:1735420/101‑1 (MQ=255)
                                                       ggggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                             <  2:3416810/99‑1 (MQ=255)
                                                        tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                            <  2:2328442/101‑1 (MQ=255)
                                                          gtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGAGGAGTGGGCTCGGGCGGATCGGTcggggccgggggattggtgtagttcgtggggtcggt                                          >  1:1991489/1‑101 (MQ=255)
                                                           tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                         >  2:1115933/1‑101 (MQ=255)
                                                           tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                         <  2:1043157/101‑1 (MQ=255)
                                                            cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                        <  1:4065570/101‑1 (MQ=255)
                                                             aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGccgggtccgtgggctcggtcgggtccgcgggatcggcggg                                       >  2:2652984/1‑97 (MQ=255)
                                                             aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtggg                                       >  2:3927896/1‑101 (MQ=255)
                                                             aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                       <  2:1567353/101‑1 (MQ=255)
                                                             aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcgggcggggccgtgggatcggtcgg                                       >  2:800779/1‑101 (MQ=255)
                                                              gCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa                                      <  2:416416/101‑1 (MQ=255)
                                                               cGTCGCCGGGTCGCGCTCATCGGCGGTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                     <  1:1304708/101‑1 (MQ=255)
                                                               cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGAt                                     >  2:3217367/1‑101 (MQ=255)
                                                                 tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTTGGATCGGTCGGATcc                                   <  1:3754101/101‑1 (MQ=255)
                                                                    ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                <  1:662528/101‑1 (MQ=255)
                                                                       ggtCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                             <  1:1053284/99‑1 (MQ=255)
                                                                         tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                           <  2:3778276/101‑1 (MQ=255)
                                                                           gcgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                         <  2:2550190/101‑1 (MQ=255)
                                                                             gcTCATCGGCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                       <  1:3515409/101‑1 (MQ=255)
                                                                             gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                       >  2:1788020/1‑101 (MQ=255)
                                                                             gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                       >  2:2771580/1‑101 (MQ=255)
                                                                                 accGACGGTTGCCCTCGGATGCGTGGCCTCGTCCGGTCCGGTCGGGCCCTTGGAATCGGCCGGG.CCGGTGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                  <  1:2653156/99‑1 (MQ=9)
                                                                                  tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                  <  1:2122477/101‑1 (MQ=255)
                                                                                  tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                  <  1:2576961/101‑1 (MQ=255)
                                                                                   cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                 <  1:825498/101‑1 (MQ=255)
                                                                                     aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt               <  1:3874441/101‑1 (MQ=255)
                                                                                     aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt               <  2:2936354/101‑1 (MQ=255)
                                                                                      tgggctctgTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGtc              <  2:3171346/93‑2 (MQ=11)
                                                                                      cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGGtcgggcggcgcccgg              >  2:784879/1‑94 (MQ=21)
                                                                                       ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGtcg             >  2:3999026/1‑99 (MQ=255)
                                                                                       ggTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg             <  1:3243246/101‑1 (MQ=255)
                                                                                        gTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg            <  2:2322869/101‑1 (MQ=255)
                                                                                         gatcggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGCGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTcggg           >  2:1273889/6‑100 (MQ=11)
                                                                                          gggCGTCGGAGGGGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt          <  2:1830900/101‑1 (MQ=255)
                                                                                           tcggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtc         <  1:4164851/98‑1 (MQ=1)
                                                                                           ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc         <  2:1344618/101‑1 (MQ=255)
                                                                                           ggCGTCGGATGCGTGGGCTCGGTCGGATCGGGCGGGTCCGTGGGATAGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtg         >  1:3971212/1‑100 (MQ=255)
                                                                                            tggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtcc        <  1:4118884/99‑2 (MQ=11)
                                                                                            gCGTGGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg        <  2:3272893/101‑1 (MQ=255)
                                                                                            gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg        >  1:237556/1‑101 (MQ=255)
                                                                                            gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg        <  1:3382300/101‑1 (MQ=255)
                                                                                             cGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTAGGGGCCGTGGGAtcgg       >  2:3301221/1‑101 (MQ=255)
                                                                                              gTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccgt      >  2:203224/1‑101 (MQ=11)
                                                                                              gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt      <  2:2047280/101‑1 (MQ=255)
                                                                                               tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc     <  1:1116024/101‑1 (MQ=255)
                                                                                                 ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg   >  1:1863547/1‑101 (MQ=255)
                                                                                                 ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg   <  1:2771770/101‑1 (MQ=255)
                                                                                                 ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg   >  2:2812592/1‑101 (MQ=255)
                                                                                                 ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCGCGAACGCACGGATCCGAGGGATCGGTCGGGTCCGTGGGAtcggtcgg   >  2:3681979/1‑101 (MQ=18)
                                                                                                  gATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGCCGGGTCCGTGGGATCGGGCGGa  >  1:3115316/1‑101 (MQ=255)
                                                                                                  gATCCGTGGGATCGGTGGGTTCCGTCGGATCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcggtcggg  <  2:4422411/101‑2 (MQ=11)
                                                                                                   |                                                                                                   
CACGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGA  >  NC_014034/221841‑222039

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG