Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,949 C→A 5.1% E236D (GAG→GAT)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9490C→A 5.1% 14.5 ‑3.8 ‑0.2183E236D (GAG→GAT) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/2):  ref base (72/101):  total (80/103)
Fisher's exact test strand distribution p-value = 2.23e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.93e-01

Only 100 of 228 total aligned reads displayed.

GGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221853‑222049
                                                                                                |                                                                                                    
ggTTTCGCCGGTGCCCGCCACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                                  <  2:4262564/101‑1 (MQ=255)
ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                                  >  1:679324/1‑101 (MQ=255)
ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                                  <  2:4535348/101‑1 (MQ=255)
ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCccgg                                                                                                  >  2:4060438/1‑101 (MQ=255)
  tttCGCTGGGGCCCGACACCGTCGGCCAGGCTGCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                                <  2:1394260/101‑1 (MQ=255)
  tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGGGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                                >  1:688902/1‑101 (MQ=255)
     cGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCgga                                                                                             <  2:3753074/101‑1 (MQ=255)
       ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                           <  2:4151906/101‑1 (MQ=255)
       ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                           <  2:4422593/101‑1 (MQ=255)
        cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                          <  2:2738637/101‑1 (MQ=255)
         ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                         >  2:2090479/1‑101 (MQ=255)
         ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                         <  1:3951046/101‑1 (MQ=255)
            tcccGACGCCGTCGCCGAGGGTTCGGTGGTGGTGGTCAGCGTCGCCGGATCGCGCTGATCGACGGTGGGCGTCGGATGCGTGGGCTGGGTCggatcggtcg                                                                                      <  2:993931/100‑1 (MQ=255)
              ccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggg                                                                                    <  2:1028536/101‑1 (MQ=255)
              ccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggg                                                                                    <  2:3461923/101‑1 (MQ=255)
                gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                  >  1:3731733/1‑101 (MQ=255)
                gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                  <  1:275987/101‑1 (MQ=255)
                 acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  1:4561306/101‑1 (MQ=255)
                 acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  2:315849/101‑1 (MQ=255)
                  cacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGAGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcag                                                                                >  2:4716178/1‑99 (MQ=255)
                    ccGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtg                                                                              <  1:583452/101‑1 (MQ=255)
                      gTCGCCGAGGCTTCGGTGTTGGCGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCggatcggtcgggtccggggg                                                                            >  1:2890939/1‑101 (MQ=255)
                        cGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggat                                                                          >  1:653551/1‑101 (MQ=255)
                        cGCCGAGGCTTCGGTGATGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggat                                                                          >  1:4147742/1‑101 (MQ=255)
                         gCCGAGCGTTCGGTGTTTGTGGTCAGCGTCGCCGGATCGCGCTGATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                         <  2:2663570/101‑1 (MQ=39)
                          ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                                        <  1:3832768/101‑1 (MQ=255)
                             aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                     >  1:3086837/1‑101 (MQ=255)
                             aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                     <  1:2147139/101‑1 (MQ=255)
                              ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcg                                                                    <  2:3888741/101‑1 (MQ=255)
                                   cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  1:4203716/101‑1 (MQ=255)
                                   cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                               <  1:2043251/101‑1 (MQ=255)
                                     gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                             <  1:4493210/101‑1 (MQ=255)
                                         tggtggTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAtcggtcgggtccgtgggctcggtcggggccgtggga                                                         >  1:3620910/1‑101 (MQ=255)
                                          ggtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                                        <  2:1620623/101‑1 (MQ=255)
                                           gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                       <  2:2066196/101‑1 (MQ=255)
                                               tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                   <  2:653487/101‑1 (MQ=255)
                                                cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                                  <  2:4075373/101‑1 (MQ=255)
                                                cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                                  >  1:4245183/1‑101 (MQ=255)
                                                 aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                                 <  1:306261/101‑1 (MQ=255)
                                                   cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                               <  2:3380427/101‑1 (MQ=255)
                                                   cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGAt                                               >  2:398253/1‑101 (MQ=255)
                                                    gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                              <  1:1162896/101‑1 (MQ=255)
                                                    gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                              >  1:1537472/1‑101 (MQ=255)
                                                      cGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                                            >  2:1789210/1‑101 (MQ=255)
                                                      cGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                                            >  1:2294911/1‑101 (MQ=255)
                                                          ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                        <  1:3248745/101‑1 (MQ=255)
                                                          ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGCGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGCATCGGTCGGATCCGTggg                                        >  1:951964/1‑101 (MQ=255)
                                                           gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                       >  1:1303145/1‑101 (MQ=255)
                                                            gtcccgcTTATCGGCGGTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                      <  2:3583826/97‑1 (MQ=21)
                                                            aTCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGTGATCGGTCGGATCCGTGGGAt                                      <  1:2688410/101‑1 (MQ=255)
                                                              tgggcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                    <  2:645868/98‑1 (MQ=255)
                                                              cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                    <  1:3475307/101‑1 (MQ=255)
                                                               gcgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                   <  1:2090600/101‑1 (MQ=255)
                                                                cgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                  <  1:4603905/101‑1 (MQ=255)
                                                                 gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                 <  2:1228745/101‑1 (MQ=255)
                                                                  cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                >  1:1373265/1‑101 (MQ=255)
                                                                  cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                >  2:4367339/1‑101 (MQ=255)
                                                                      tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                            <  2:1099600/101‑1 (MQ=255)
                                                                       cGAGGGGGGGCGGCGGGTGCGTGGGCTCGGTCGGATGGGTCGGGTCCGGGGGATGGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                           <  2:458265/101‑1 (MQ=18)
                                                                       cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                           >  1:3524447/1‑101 (MQ=255)
                                                                          cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                        >  2:573512/1‑101 (MQ=255)
                                                                          cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                        <  1:3359561/101‑1 (MQ=255)
                                                                             tGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa                     >  1:3406263/1‑101 (MQ=255)
                                                                              gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                    <  1:1789319/101‑1 (MQ=255)
                                                                              gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                    <  1:4367586/101‑1 (MQ=255)
                                                                               tccggggtaTGTGTCGGATCCGTGGGATCCGTCGGGTCCGGGGGATCCGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                   <  1:4708120/93‑1 (MQ=11)
                                                                                  gttgtgtccGTGAGATCGGTTGTATCCGTGGGGTCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccgt                <  2:1363514/93‑1 (MQ=14)
                                                                                  gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                >  2:2356309/1‑101 (MQ=255)
                                                                                  gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                <  2:3083794/101‑1 (MQ=255)
                                                                                   tggggttcGTGGGATCCGTGGGATCGGTCGGATCCGGGGGATTGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  1:2789855/94‑1 (MQ=11)
                                                                                   gggggtGCGTTGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  1:2347087/96‑1 (MQ=11)
                                                                                   ggggggcggggggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc               <  1:914173/91‑1 (MQ=14)
                                                                                     ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             >  1:1686012/1‑101 (MQ=255)
                                                                                     ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             <  1:1133834/101‑1 (MQ=255)
                                                                                       gtccGTCGGATCGGTGGGACCGGGCCGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           <  2:3216736/98‑1 (MQ=11)
                                                                                       gtGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           <  1:3649937/100‑1 (MQ=255)
                                                                                       atccGTGGGATAGGTCGGGTCCGTGGGATCGGTGGGATCGGTGGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGAtcggtcgggt           <  1:4530219/98‑3 (MQ=1)
                                                                                       aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           >  1:996551/1‑101 (MQ=255)
                                                                                       aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAATGGTGGGac           >  2:2906720/1‑100 (MQ=255)
                                                                                        tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc          <  1:1315873/101‑1 (MQ=255)
                                                                                        tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGCGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc          <  1:4017661/101‑1 (MQ=255)
                                                                                         gCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGGCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc         <  2:1999038/101‑1 (MQ=11)
                                                                                         cggtcggCTCCGTCGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc         <  2:944253/96‑1 (MQ=11)
                                                                                         cggTGGGATCGGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc         <  1:10081/99‑1 (MQ=11)
                                                                                           gggggCTCGGTCGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGTATCCGt       <  2:4600350/99‑1 (MQ=255)
                                                                                           gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCgg       >  2:2160633/1‑100 (MQ=255)
                                                                                           gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt       <  1:137959/101‑1 (MQ=255)
                                                                                            tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      >  2:1747828/1‑101 (MQ=255)
                                                                                            tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      <  1:1036279/101‑1 (MQ=255)
                                                                                            tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      >  2:2887691/1‑101 (MQ=255)
                                                                                            tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      >  2:2058830/1‑101 (MQ=255)
                                                                                            tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGGtcggtcggggccggg      >  2:3425017/1‑94 (MQ=14)
                                                                                            tGGGATCGGGCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCgcc      <  2:888092/101‑3 (MQ=11)
                                                                                             gggGTCGGGTGGGTCCGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg     <  1:2052917/101‑3 (MQ=11)
                                                                                             gggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGg     <  1:3133588‑M2/101‑2 (MQ=255)
                                                                                             gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg     >  1:2581653/1‑99 (MQ=1)
                                                                                             gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg     >  1:527460/1‑99 (MQ=11)
                                                                                              ggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgg    >  2:3809729‑M2/1‑99 (MQ=255)
                                                                                              ggCGCGATCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGATCCGTGgg    <  2:1905128‑M2/101‑3 (MQ=255)
                                                                                               tcTCCGTGGGATCGGTCGTGTCCGGGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa   <  2:4301037/100‑1 (MQ=11)
                                                                                                |                                                                                                    
GGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221853‑222049

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