Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,958 A→G 10.9% D233D (GAT→GAC)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9580A→G10.9% 36.1 ‑0.2 ‑1.1162D233D (GAT→GAC) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/7):  ref base (79/65):  total (90/72)
Fisher's exact test strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.47e-02

Only 100 of 420 total aligned reads displayed.

CGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221858‑222058
                                                                                                    |                                                                                                    
   cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                                   <  2:2738637/101‑1 (MQ=255)
    ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                                  <  1:3951046/101‑1 (MQ=255)
      tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                                >  1:2066317/1‑101 (MQ=255)
           gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                           <  1:275987/101‑1 (MQ=255)
            acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                          <  2:315849/101‑1 (MQ=255)
            acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                          <  1:635191/101‑1 (MQ=255)
             cacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGAGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcag                                                                                         >  2:4716178/1‑99 (MQ=255)
                cGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCACCGACGGTGGGCGTCGGATGCGTGGGCTCGGGGgggtcggtcggggccgtgg                                                                                      >  1:4600603/1‑101 (MQ=255)
                    gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                                  <  2:1193024/101‑1 (MQ=255)
                      cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcgg                                                                                >  2:306227/1‑101 (MQ=255)
                      cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCACCGGTGGGCGGCGGATGCGGGGGCTCGGTCgggtcgggcgggtccgggggatggg                                                                                >  1:1054892/1‑101 (MQ=255)
                        aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                              >  1:3086837/1‑101 (MQ=255)
                         ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcg                                                                             <  2:3888741/101‑1 (MQ=255)
                           cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                                           >  1:420818/1‑101 (MQ=255)
                           cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                                           <  1:4475755/101‑1 (MQ=255)
                           cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                                           <  1:683070/101‑1 (MQ=255)
                              gggTGTTTGGGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                                        <  1:3391851/100‑1 (MQ=255)
                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                                        <  1:2043251/101‑1 (MQ=255)
                              cGGTGTTGGAGGTCAGCGTCGCCGGATCGCGCTCCTCGACGGTGGGCGTCGGATGCGTGGGCTCGGCCggatcggtcgggtccgtgggatcggtggggtcc                                                                        >  1:380878/1‑101 (MQ=39)
                                    tggtggTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAtcggtcgggtccgtgggatcggtcgggtccggggga                                                                  >  1:3139724/1‑101 (MQ=255)
                                             gCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa                                                         <  1:2998309/101‑1 (MQ=255)
                                              cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                                        <  2:3380427/101‑1 (MQ=255)
                                              cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGAt                                                        >  2:398253/1‑101 (MQ=255)
                                               gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                                       <  1:1144811/101‑1 (MQ=255)
                                                tCGCCGGAGCGCGCTCAGCGACGGTGGGCGTCGGATGCGTGGTCCCGTTCGGATCGTTTGGGTCCGTGGGATCGGTCGGTTCTGTGGGATCGGTCGGATcc                                                      <  2:380834/101‑1 (MQ=25)
                                                 cGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                                                     >  2:1789210/1‑101 (MQ=255)
                                                   ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                   >  2:1845947/1‑101 (MQ=255)
                                                     ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGCGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGCATCGGTCGGATCCGTggg                                                 >  1:951964/1‑101 (MQ=255)
                                                       aTCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                               <  1:2865311/101‑1 (MQ=255)
                                                           cgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                           >  2:4286196/1‑101 (MQ=255)
                                                            gcTCATCGACGGTGGGCGTCGAATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                          >  1:3492552/1‑101 (MQ=255)
                                                             cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                         >  1:1373265/1‑101 (MQ=255)
                                                              tCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                                        >  2:137943/1‑101 (MQ=255)
                                                                  cGAGGGGGGGCGGCGGGTGCGTGGGCTCGGTCGGATGGGTCGGGTCCGGGGGATGGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                    <  2:458265/101‑1 (MQ=18)
                                                                  cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGcc                                    >  2:4147656/1‑101 (MQ=255)
                                                                           gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGtcg                           >  2:2445195/1‑101 (MQ=255)
                                                                            cGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCAGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                          >  1:4019256/1‑101 (MQ=255)
                                                                                ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                      >  1:2106200/1‑101 (MQ=255)
                                                                                  gtGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                    <  1:3649937/100‑1 (MQ=255)
                                                                                  atccGTGGGATAGGTCGGGTCCGTGGGATCGGTGGGATCGGTGGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGAtcggtcgggt                    <  1:4530219/98‑3 (MQ=1)
                                                                                    gtttGGGATCTGTGGGGTCCGTCGGATCCGTGGGGTCCGTAGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                  <  1:1952836/98‑1 (MQ=11)
                                                                                    gCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGGCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                  <  2:1999038/101‑1 (MQ=11)
                                                                                    cggggggATCGGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGGTCCGTCGGATcc                  <  1:2592898/97‑1 (MQ=11)
                                                                                    cggTGGGATCGGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                  <  1:10081/99‑1 (MQ=11)
                                                                                      ttGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCGGACGGATCGGTCCGGTCCGTGGCAtcggccggggccgg                >  1:4672493/2‑95 (MQ=2)
                                                                                      gTGGGGTCCGGCGGAGCGGGCGGTTCGGGGGGACCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                <  1:173033/101‑1 (MQ=11)
                                                                                      gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCgg                >  2:2160633/1‑100 (MQ=255)
                                                                                      gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCgc                >  1:888166/1‑100 (MQ=11)
                                                                                       tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg               >  2:2058830/1‑101 (MQ=255)
                                                                                       tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGGtcggtcggggccggg               >  2:3425017/1‑94 (MQ=14)
                                                                                       tGGGATCGGGCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCgcc               <  2:888092/101‑3 (MQ=11)
                                                                                       gggaatcggtaggTTCGGCGGGGTCCGTAGG.TTGGGCCGGTTCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGtcg              <  2:717757/96‑3 (MQ=2)
                                                                                         ggggcggtgggATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg             <  1:23015/97‑1 (MQ=11)
                                                                                           gtcggtgggATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt           <  1:1297242/100‑1 (MQ=11)
                                                                                           gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGGCCGTCGGATCCGTGGGAt           >  2:1505068/8‑101 (MQ=11)
                                                                                           gtccgtcggATCGTTGGGATCCGTAGTATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt           <  1:1335283/97‑1 (MQ=11)
                                                                                           atcaggcggAGCGCTCGGGCCCGTCGGAGCGGTCGGGCCCGTCGGACCCGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt           <  2:2927323/95‑1 (MQ=2)
                                                                                            tcggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc          >  1:1416748/1‑101 (MQ=255)
                                                                                            tcgctgggaccgggcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtc          <  2:2412796/95‑1 (MQ=14)
                                                                                            tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGAGCCGTGGGATCGGTCGGATCCGTGGGCtc          >  2:1273390/7‑101 (MQ=11)
                                                                                            tccgtcggATCGGTGGTATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc          <  2:168566/98‑1 (MQ=11)
                                                                                            tccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGGTCGGTCGGGCCCGGAGGATCGGTGGGGTCCGTGGGATAGATCGGGTCAGTGGGAta          >  2:2333490/4‑100 (MQ=11)
                                                                                            tccgggggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtc          <  1:3433806/95‑1 (MQ=11)
                                                                                            gccgtcggGTCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc          <  1:475339/98‑1 (MQ=11)
                                                                                             cggtgggATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg         <  1:334272/101‑1 (MQ=11)
                                                                                             cggtcggGTCCGTCGGTTCCTTTGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg         <  2:366152/101‑1 (MQ=11)
                                                                                             cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg         <  1:4147902/101‑1 (MQ=255)
                                                                                             ccgttggatcggtggggtccctgggctcggtcgggtccgtaGGATCTGTCGGGTCCGTCGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg         <  2:3512555/96‑1 (MQ=11)
                                                                                             ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGAGCGGTGGGATCCGTgggggcg         >  2:98262/6‑97 (MQ=11)
                                                                                             ccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg         <  2:1019114/99‑1 (MQ=11)
                                                                                             ccgtcggAGCGGTCGGGTCCGTTGGCTCGGTCGGGTCCGTTGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcc         <  2:4664642/99‑2 (MQ=11)
                                                                                              ggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg        >  2:3254181/1‑101 (MQ=255)
                                                                                              cgtgggatcgttcgggtccgtgggagcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg        <  2:419704/97‑1 (MQ=11)
                                                                                              cgtgggatcggtcgggtccgtgggatcggtggggtccgtgGGATCGGTCGGGTCGGTGGCGTCGGTCGGGTCGGTGGGGGCGGTGGGTGCCGGggggcccg        >  2:4504981/5‑96 (MQ=0)
                                                                                               gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCGGTCGGTTCCGTGGGAtcggt       >  2:2656278/4‑101 (MQ=11)
                                                                                               gtgggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGAtcggg       >  1:2516983/4‑100 (MQ=11)
                                                                                               gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgag       >  1:4305873/1‑99 (MQ=255)
                                                                                               gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAGCGGTCGGGTCCGTGGGTTCGGTCCGGTCCGCCGGAtcggt       >  1:1557745/1‑101 (MQ=14)
                                                                                                tgggatcggtcgggtgcgtgggatcgggcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc      <  1:235564/99‑1 (MQ=11)
                                                                                                tgggatcggtcgggtccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgggtccgtc      <  1:3561128‑M2/99‑10 (MQ=255)
                                                                                                tgggatcggtcgggtccgtgggatcggtcgggtccgtgggATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtg      >  2:403493‑M2/3‑92 (MQ=255)
                                                                                                tgggatcggtcgggtccgtgggatcggtcgggtccgtcGGGTCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTTGGGTCGGTCGGGTCCGTGGGAacagtc      >  2:4320415/3‑101 (MQ=11)
                                                                                                tgggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtc      <  2:4378221/99‑1 (MQ=11)
                                                                                                tgggatcgggcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtc      <  1:1905240/99‑1 (MQ=11)
                                                                                                tcggATCGGTGGGACCCGTGGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc      <  2:2081917/101‑1 (MQ=11)
                                                                                                tcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtc      <  1:572425/101‑1 (MQ=11)
                                                                                                tcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGACCGTGGGGTCGGGGGGGCCCGTCGGAtcggcg      >  2:1587716/1‑99 (MQ=14)
                                                                                                 gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCGGGCGGGTGCGTGGGATCGGTCGGGACCGTGGGAtcggcgg     >  2:604262/2‑98 (MQ=11)
                                                                                                 gggatcggtcgggtccgtcggatcggtcgggtccgtgGGATCGGTCGGGCCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggttg     >  1:15932/2‑99 (MQ=11)
                                                                                                 cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg     >  1:4378468/1‑101 (MQ=255)
                                                                                                 cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggggggtcg     >  2:3621432/1‑93 (MQ=11)
                                                                                                 cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg     <  1:3994935/101‑1 (MQ=11)
                                                                                                  ggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGGGGGATCGGTCGGAGCCGTGGGAtcggtcgg    >  2:658852/1‑101 (MQ=11)
                                                                                                   gatcggtgggatccgtgggatcggtcgggtccgtggGATCGGTCGGGTCCGTCGGGTCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCGGTCcgg   >  2:4359681/1‑101 (MQ=11)
                                                                                                    atcggttgggtccgggggatcggtcgggtccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt  <  1:774913/101‑1 (MQ=11)
                                                                                                    atcggtcgggtccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGATCTGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGat  >  2:2890763‑M2/1‑99 (MQ=255)
                                                                                                    atcggtcgggtccgtgggatcggtcgggtccgtgggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGat  <  2:1051122‑M2/101‑3 (MQ=255)
                                                                                                    atcggtcgggtccgtgggatcggtcgggtccgtgggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGat  >  2:2579453‑M2/1‑99 (MQ=255)
                                                                                                    atcggtcgggtccgtcggatcggtcgggtccgtgggATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGGGGGAGCGGTCGGGt  >  2:4460828/1‑101 (MQ=11)
                                                                                                    atcggtcggatccgtgggatcggtcggatccgtaggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGt  >  1:2250389/1‑101 (MQ=11)
                                                                                                    |                                                                                                    
CGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221858‑222058

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG