Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,949 C→A 6.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 6.1% 43.6 ‑2.3 ‑1.6408E236D (GAG→GAT) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (18/7):  ref base (203/178):  total (221/187)
Fisher's exact test strand distribution p-value = 9.59e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.54e-02

Only 100 of 469 total aligned reads displayed.

AGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221849‑222049
                                                                                                    |                                                                                                    
aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGc                                                                                                      >  2:4100359/1‑101 (MQ=255)
aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGGGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGCGGGCGTCGGATGCGTGGGc                                                                                                      >  1:2918801/1‑101 (MQ=255)
  ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                                    >  2:8350659/1‑101 (MQ=255)
   cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcg                                                                                                   <  1:7151452/101‑1 (MQ=255)
    ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                                  >  2:3183407/1‑101 (MQ=255)
    ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                                  >  2:3249726/1‑101 (MQ=255)
    ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                                  >  1:1687854/1‑101 (MQ=255)
     gTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggt                                                                                                 <  2:1682818/101‑1 (MQ=255)
     gTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGAGGTCAGCGTCGCCGGGTCGCGCTCATCGACGTGGGGAGGCGGATGCGTGGGCttggc                                                                                                 >  1:1014011/1‑100 (MQ=255)
      tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                                <  2:7009963/101‑1 (MQ=255)
       ttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtcg                                                                                               >  1:312522/1‑101 (MQ=255)
        tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtcgg                                                                                              <  1:5548821/101‑1 (MQ=255)
        tCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtcgg                                                                                              <  2:3686521/101‑1 (MQ=255)
          gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                            >  1:5797805/1‑101 (MQ=255)
          gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                            >  2:778923/1‑101 (MQ=255)
           ccGGGGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                           <  1:6226729/101‑1 (MQ=255)
               tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                       <  2:1313197/101‑1 (MQ=255)
                 cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggccgg                                                                                     >  2:8871877/1‑101 (MQ=255)
                     gcacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  2:2483917/100‑1 (MQ=255)
                     acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  1:491290/101‑1 (MQ=255)
                     acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGATCATCGACGGGGGGCGTCGGATGCGTGGGCTCGGTCggatgggtcgggtcc                                                                                 >  2:6860260/1‑101 (MQ=255)
                      cacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccg                                                                                >  1:2755938/1‑101 (MQ=255)
                         cGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGGCggatcggtcggggccgtgg                                                                             >  1:2173257/1‑101 (MQ=255)
                           tCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggga                                                                           <  1:812633/101‑1 (MQ=255)
                            cGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggat                                                                          >  2:146536/1‑101 (MQ=255)
                              ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                                        <  1:5352247/101‑1 (MQ=255)
                              ccGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                                        <  1:6383203/101‑1 (MQ=255)
                               cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcgg                                                                       >  1:6079805/1‑101 (MQ=255)
                                gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                                      >  1:6478024/1‑101 (MQ=255)
                                gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcggt                                                                      >  2:6633572/1‑101 (MQ=255)
                                  ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgggggatcggtcg                                                                    >  1:2241553/1‑101 (MQ=255)
                                     ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                 <  2:633877/101‑1 (MQ=255)
                                      tCGGTGTTGGGGGTCAGCGTCGCCGGATCGCGCTCATCGGCGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtc                                                                <  1:4100478/101‑1 (MQ=255)
                                       cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatgggtcgggtcc                                                               >  2:505912/1‑101 (MQ=255)
                                        ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                              <  2:5359993/101‑1 (MQ=255)
                                        ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                              >  1:2579475/1‑101 (MQ=255)
                                        ggTGTGGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                              >  1:5455665/1‑101 (MQ=255)
                                         gtggTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                             <  2:3812909/101‑1 (MQ=255)
                                           gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGTGCACGGTCggatcggtcgggtccgttggatcggtcgggtcctggg                                                           >  2:2269917/1‑97 (MQ=38)
                                               gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGGTCggtcgggtccgtgggatcggtcggggccgtgggatc                                                       >  2:6093915/1‑101 (MQ=255)
                                               gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                       >  2:9195370/1‑101 (MQ=255)
                                               gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                       >  1:5612723/1‑101 (MQ=255)
                                                tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTTGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                      <  1:729517/101‑1 (MQ=255)
                                                 ggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                     <  1:768105/101‑1 (MQ=255)
                                                  gtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGtcgggtccgtggggtcggtcgggtccgtcggatcggc                                                    >  1:8601625/1‑100 (MQ=255)
                                                   tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGGCgggtccctgggatcggtcgggtccctgggatcggtc                                                   >  2:808754/1‑101 (MQ=38)
                                                    cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                                  >  2:5789967/1‑101 (MQ=255)
                                                     aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                                 <  2:4339297/101‑1 (MQ=255)
                                                       cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGAt                                               >  1:9003921/1‑101 (MQ=255)
                                                        gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                              <  2:5455480/101‑1 (MQ=255)
                                                        gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                              >  1:1178911/1‑101 (MQ=255)
                                                             cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGGGGGGTCGGTCGGGTTCGTGGGCTCGGTCGGGTCCGGgg                                         >  2:3631121/1‑101 (MQ=255)
                                                              ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                        <  2:3893058/101‑1 (MQ=255)
                                                               gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                       <  2:7500611/101‑1 (MQ=255)
                                                                  cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                    <  2:6503855/101‑1 (MQ=255)
                                                                     gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                 >  1:3665841/1‑101 (MQ=255)
                                                                     gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                 >  1:4429603/1‑101 (MQ=255)
                                                                      cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                >  1:1869642/1‑101 (MQ=255)
                                                                       tcagcggtggTGGGCGTCGGATGCGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcggtcgg                               <  2:6711918/93‑1 (MQ=14)
                                                                       tCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAgcggtcgg                               >  2:1066431/1‑101 (MQ=255)
                                                                       tCAGCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                               >  1:8811806/1‑101 (MQ=255)
                                                                        cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCggg                              >  1:5207715/1‑101 (MQ=255)
                                                                         aTCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCgggc                             >  1:7831206/1‑100 (MQ=255)
                                                                          tgggacggttgCGTCGGATGAATCGGGTCCATCCGATCGCTCGGGTCCGTAGGAGCGGTCCGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                            <  2:1550486/91‑1 (MQ=2)
                                                                           cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                           >  2:4492195/1‑101 (MQ=255)
                                                                            gACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                          >  2:4905161/1‑101 (MQ=255)
                                                                            gACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGGGGGATCGGGCGGGTCCGTGGGAGCGGTCGGGGCCg                          >  1:4218842/1‑101 (MQ=21)
                                                                             aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCgc                         >  1:5061624/1‑100 (MQ=255)
                                                                              cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                        >  1:5889622/1‑101 (MQ=255)
                                                                                   ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                   <  1:3930986/101‑1 (MQ=255)
                                                                                      gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGAtcggt                >  2:3459553/1‑101 (MQ=255)
                                                                                      gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                >  1:4779354/1‑101 (MQ=255)
                                                                                       tgcgATGCGTCGCCTCCCTGGCATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  2:8684869/98‑1 (MQ=11)
                                                                                       tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTGGGAtcggtc               >  2:1749216/1‑101 (MQ=255)
                                                                                       cttccgatctGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  1:1622490/92‑1 (MQ=1)
                                                                                        cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg              >  2:8550603/1‑101 (MQ=255)
                                                                                        cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg              >  1:4683012/1‑101 (MQ=255)
                                                                                        cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg              >  2:2849384/1‑101 (MQ=255)
                                                                                         ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             >  2:4664734/1‑101 (MQ=255)
                                                                                         ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGCGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg             <  1:149785/101‑1 (MQ=255)
                                                                                           gtGCGTGGGGTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           <  1:1346800/100‑1 (MQ=255)
                                                                                           atccGTGGGATCGGTCGGGTGCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcggtcgggg           >  1:8221704/4‑99 (MQ=11)
                                                                                           agcggTGTGCTCGGTGGGATCGGTCGGATCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGAt           <  2:8270357/97‑1 (MQ=11)
                                                                                           aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           >  1:6957447/1‑101 (MQ=255)
                                                                                           aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           >  1:905254/1‑101 (MQ=255)
                                                                                           aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt           >  1:8328160/1‑101 (MQ=255)
                                                                                            tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc          <  2:5889387/101‑1 (MQ=255)
                                                                                            tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc          <  1:6864874/101‑1 (MQ=255)
                                                                                             gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc         >  2:1975441/1‑101 (MQ=255)
                                                                                             cggTGGGATCCGTCGGATCGGTGGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTcc         <  2:3123418/99‑1 (MQ=11)
                                                                                              cGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg        >  1:965321/1‑101 (MQ=255)
                                                                                              cGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg        <  2:212049/101‑1 (MQ=255)
                                                                                              cGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg        >  1:8463941/1‑101 (MQ=11)
                                                                                               gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt       <  1:793077/101‑1 (MQ=255)
                                                                                               gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt       >  1:3549314/1‑101 (MQ=11)
                                                                                                tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtc      <  1:6712656/101‑2 (MQ=11)
                                                                                                ggggGTCGGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      <  2:7256576/100‑1 (MQ=255)
                                                                                                ggggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      <  2:8990829/100‑1 (MQ=255)
                                                                                                 gggCGCGGTCGGAGCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGg     >  2:2698007‑M2/1‑100 (MQ=255)
                                                                                                  ggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgg    <  1:319546‑M2/101‑3 (MQ=255)
                                                                                                    |                                                                                                    
AGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221849‑222049

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG