Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,942 G→C 5.6% H239D (CAT→GAT)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9420G→C 5.6% 14.7 ‑1.5 ‑0.0184H239D (CAT→GAT) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/3):  ref base (86/87):  total (94/90)
Fisher's exact test strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.47e-01

Only 100 of 222 total aligned reads displayed.

ACGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCC  >  NC_014034/221842‑222042
                                                                                                    |                                                                                                    
acGGGCGAGCGGGTTTCGCCGGTGCCCGACGCCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGGGATCATCGACGGTGGGCGTCGGATg                                                                                                      <  2:4106941/101‑1 (MQ=255)
  ggTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCg                                                                                                    >  2:3797485/1‑101 (MQ=255)
   gTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGt                                                                                                   >  2:4146334/1‑101 (MQ=255)
    tCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTg                                                                                                  <  2:290425/101‑1 (MQ=255)
     cGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTgg                                                                                                 <  1:1248340/101‑1 (MQ=255)
     cGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTgg                                                                                                 >  2:2705006/1‑101 (MQ=255)
      gAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTggg                                                                                                <  1:4278883/101‑1 (MQ=255)
       aGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGc                                                                                               >  2:88906/1‑101 (MQ=255)
          cGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcg                                                                                            <  2:434375/101‑1 (MQ=255)
           ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCCCATCGACGGTGGGCGTCGGGTGGGTGGGCtcgg                                                                                           >  1:4343487/1‑101 (MQ=255)
            gTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggt                                                                                          >  2:1191249/1‑101 (MQ=255)
             tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                         <  2:1024126/101‑1 (MQ=255)
             tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                         <  1:3575717/101‑1 (MQ=255)
             tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                         <  1:2248286/101‑1 (MQ=255)
                cGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCgga                                                                                      <  2:2027219/101‑1 (MQ=255)
                 gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                     <  1:2008377/101‑1 (MQ=255)
                 gCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggat                                                                                     >  2:3217961/1‑101 (MQ=255)
                   cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGGGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                   <  1:4537028/101‑1 (MQ=255)
                    ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                  >  2:4720356/1‑101 (MQ=255)
                    ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                  >  1:1831266/1‑101 (MQ=255)
                     gTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggt                                                                                 >  1:470096/1‑101 (MQ=255)
                      tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                <  2:541756/101‑1 (MQ=255)
                       gCCCGACACCGTCGCCGAGGCTTCGGGGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGCGCCCGGTCggatcggtcg                                                                               >  2:1895466/1‑101 (MQ=255)
                           gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGGGGGCTCGGTCggatcggtcgggtc                                                                           >  2:4701216/1‑101 (MQ=255)
                            acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                          <  2:664049/101‑1 (MQ=255)
                                 gTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                     <  2:1463070/101‑1 (MQ=255)
                                    gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                  <  2:2754476/101‑1 (MQ=255)
                                    gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                  <  2:4545162/101‑1 (MQ=255)
                                       gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                               >  1:2002773/1‑101 (MQ=255)
                                       gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                               <  2:2657097/101‑1 (MQ=255)
                                        gggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                              <  2:4418937/100‑1 (MQ=255)
                                        aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                              <  2:4012253/101‑1 (MQ=255)
                                        aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGGGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggggtcggtc                                                              >  1:3037129/1‑101 (MQ=255)
                                           cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                           >  2:4639133/1‑101 (MQ=255)
                                            ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                          <  1:3797770/101‑1 (MQ=255)
                                             tCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgggcgggtccgtggggtcggtcgggtc                                                         >  1:927272/1‑101 (MQ=255)
                                             tCGGTGTTGGGTGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgggcgggtccgtgggatcggtcgggtc                                                         <  2:156938/101‑1 (MQ=255)
                                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                        >  1:4428982/1‑101 (MQ=255)
                                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                        <  1:3016054/101‑1 (MQ=255)
                                                gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggggccgtgggatcggtcgggtccgt                                                      >  1:3411082/1‑101 (MQ=255)
                                                   ttGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                   >  1:4638630/1‑101 (MQ=255)
                                                   ttGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                   >  1:2319833/1‑101 (MQ=255)
                                                    tggtggTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAtcggtcgggtccgtgggatcggtggggtgcgtggga                                                  >  2:215681/1‑101 (MQ=255)
                                                      gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                <  2:3697826/101‑1 (MQ=255)
                                                       tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                               >  2:2867401/1‑101 (MQ=255)
                                                        ggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGAGGCGTGGGCTCGGTCGGATCGGtcgggtccgtgggatcggtcgggtccgtgggatcgg                                              <  1:2083003/101‑1 (MQ=255)
                                                        ggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGtcgggtccgtgggatcggtcgggtccgtgggatcgg                                              >  2:2214221/1‑101 (MQ=255)
                                                          gcAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                            <  2:1656611/100‑1 (MQ=255)
                                                           cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcggggccgtgggatcggtcg                                           >  2:2802623/1‑101 (MQ=255)
                                                           cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcgggcgggtccgtggggtcggtcg                                           >  2:2374426/1‑101 (MQ=255)
                                                            aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                          <  2:4505861/101‑1 (MQ=255)
                                                              cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                        >  1:1030726/1‑101 (MQ=255)
                                                               gTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                       <  2:819951/101‑1 (MQ=255)
                                                               gTCGCCGGATCGCGCTCATCGACGGTGGGCGCCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTTGGGTCCGTGGGTTCGGTCGGGTc                                       >  1:2679678/1‑101 (MQ=37)
                                                                tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                      <  1:4146651/101‑1 (MQ=255)
                                                                   ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGTg                                   >  2:2190263/1‑101 (MQ=255)
                                                                   ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                   >  1:2768716/1‑101 (MQ=255)
                                                                   ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                   >  2:169117/1‑101 (MQ=255)
                                                                   ccGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                   <  1:2608176/101‑1 (MQ=255)
                                                                     ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                 >  1:4045183/1‑101 (MQ=255)
                                                                     ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                 >  1:687157/1‑101 (MQ=255)
                                                                      gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                <  1:1594880/101‑1 (MQ=255)
                                                                        tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                              <  1:4821728/101‑1 (MQ=255)
                                                                        tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGGGGGAtg                              >  2:3527462/1‑100 (MQ=255)
                                                                        tCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                              <  1:528603/101‑1 (MQ=255)
                                                                         cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                             >  1:3199065/1‑101 (MQ=255)
                                                                         cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                             <  1:719765/101‑1 (MQ=255)
                                                                         cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                             >  2:3295084/1‑101 (MQ=255)
                                                                         cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGCGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCCTGCGGGTCCGTCGGAccg                             >  2:858960/1‑101 (MQ=25)
                                                                          gcgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                            >  2:1598434/1‑101 (MQ=255)
                                                                           cgcTCATGGACGGTGGGCGTCGGATGCTTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                           <  1:215715/101‑1 (MQ=255)
                                                                           cgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                           >  1:2495687/1‑101 (MQ=255)
                                                                            gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                          <  2:2566223/101‑1 (MQ=255)
                                                                             cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                         <  1:1433662/101‑1 (MQ=255)
                                                                                aTCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGt                      <  1:4867085/101‑1 (MQ=255)
                                                                                  ccgtgggatccGTCGGGTCCGGGGGATCCGTCGGATCTGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTcc                    <  2:3954526/91‑1 (MQ=11)
                                                                                  cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                    >  2:4447330/1‑101 (MQ=255)
                                                                                   gACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                   <  2:3061190/101‑1 (MQ=255)
                                                                                   gACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGACGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                   >  2:4204723/1‑101 (MQ=255)
                                                                                    gtggggtcggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCTGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGt                  >  2:4904909/10‑101 (MQ=11)
                                                                                    gcGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                  <  2:427864/100‑1 (MQ=255)
                                                                                     cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                 <  1:288857/101‑1 (MQ=255)
                                                                                      gggggcggTCGGTTCCGTGGGGTCGGGGGGGTCGGTGGGGTCGGTGGGATCGGTGGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                <  2:1831133/94‑3 (MQ=2)
                                                                                      ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                >  2:336155/1‑101 (MQ=255)
                                                                                      ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTgg                >  2:2460259/1‑101 (MQ=255)
                                                                                        tGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTgcat              >  1:1523878/1‑98 (MQ=255)
                                                                                          tcggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCgccgggtc            >  2:3348417/4‑94 (MQ=11)
                                                                                           gCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg           <  1:1895603/101‑1 (MQ=255)
                                                                                           cggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTtggggcc           >  1:843734/3‑97 (MQ=11)
                                                                                            cGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTgcatcaag          <  2:1523748/101‑8 (MQ=255)
                                                                                             gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt         >  1:2181092/1‑101 (MQ=255)
                                                                                              tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc        >  1:3720793/1‑101 (MQ=255)
                                                                                              tCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc        >  1:3150941/1‑101 (MQ=255)
                                                                                              gcGGATCGGTCGGGTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc        <  2:3175577/100‑1 (MQ=11)
                                                                                               cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg       <  1:3350266/101‑1 (MQ=255)
                                                                                                ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGtcggtcgg      >  2:1120016/1‑101 (MQ=255)
                                                                                                 gATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa     >  1:982663/1‑101 (MQ=255)
                                                                                                   tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc   <  1:556772/101‑1 (MQ=255)
                                                                                                    gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc  <  2:3432393/101‑1 (MQ=255)
                                                                                                    gCGTGGGCTCGGGCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc  <  2:1221517/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
ACGGTCGAGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCC  >  NC_014034/221842‑222042

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