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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9580A→G 3.3% 39.9 ‑1.5 ‑1.1758D233D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/10):  ref base (349/384):  total (364/394)
Fisher's exact test strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.33e-02

Only 100 of 1376 total aligned reads displayed.

CGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221858‑222058
                                                                                                    |                                                                                                    
ggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggga                                                                                                      >  1:4198499‑M1/99‑101 (MQ=255)
gatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                                                                      <  1:2645164‑M1/3‑1 (MQ=255)
 gtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggat                                                                                                     <  1:1368218‑M1/4‑1 (MQ=255)
  tccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                                                    >  1:1299333‑M1/97‑101 (MQ=255)
   ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                                   <  1:5430885‑M1/6‑1 (MQ=255)
   ccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                                                                                   <  1:2310487‑M1/6‑1 (MQ=255)
     gtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggt                                                                                                 <  1:19959‑M1/8‑1 (MQ=255)
      tgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                                >  1:3891263‑M1/93‑101 (MQ=255)
        ggatcggtcggttccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtaggatcggtcgg                                                                                              >  2:1368130‑M1/91‑101 (MQ=255)
         gatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtaggatcggtcggg                                                                                             <  1:3975945‑M1/12‑1 (MQ=255)
           tcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtc                                                                                           <  1:2212209‑M1/14‑1 (MQ=255)
            gcacCGTCGCCGGGGCTTTGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                          <  2:2311967/100‑1 (MQ=255)
             cacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccg                                                                                         <  2:4628118/101‑1 (MQ=255)
               tcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggggccgtg                                                                                       <  2:2100456‑M1/18‑1 (MQ=255)
               tcgggtccgggggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtg                                                                                       >  2:3514845‑M1/84‑101 (MQ=255)
                      cgtgggatcggtcgggtccgtggggtcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcgg                                                                                >  1:153310‑M1/77‑101 (MQ=255)
                      cgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcgg                                                                                <  2:1223096‑M1/25‑1 (MQ=255)
                       gtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggc                                                                               <  1:712352‑M1/26‑2 (MQ=255)
                        tgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc                                                                              >  2:5739098‑M1/75‑101 (MQ=255)
                         gggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggacccgtgggatcggtcgggtccgtgggatcggtcg                                                                             <  1:4979941‑M1/28‑1 (MQ=255)
                                gtgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGGGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgt                                                                      >  2:5383788/1‑101 (MQ=255)
                                 tcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtg                                                                     <  1:120866‑M1/36‑1 (MQ=255)
                                  ggggtctgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                    >  1:5204991‑M1/65‑101 (MQ=255)
                                     gtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                                                 <  1:4860614‑M1/40‑1 (MQ=255)
                                     gtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATcggtcgggtccgtgggatcggtcgggtccgtgggat                                                                 >  1:2545213‑M1/62‑101 (MQ=255)
                                      tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtggggtc                                                                <  2:1587510‑M1/41‑4 (MQ=255)
                                      tccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                >  2:1118319‑M1/61‑101 (MQ=255)
                                      tccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                >  1:932886‑M1/61‑101 (MQ=255)
                                       gggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                               <  2:4328694/100‑1 (MQ=255)
                                       cggtgggatcggtcgggtcggtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                               >  1:4305837‑M1/60‑101 (MQ=255)
                                       ccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                               >  1:4105015‑M1/60‑101 (MQ=255)
                                       ccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGgtcggatccgtgggatcggtcgggtccgtgggatcg                                                               >  2:5628527‑M1/60‑101 (MQ=255)
                                       ccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                               >  1:806154‑M1/60‑101 (MQ=255)
                                          tgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                            <  1:4524903‑M1/45‑1 (MQ=255)
                                          tgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGtcgggtccgtgggatcggtcgggtccgtgggatcggtg                                                            <  2:5096728‑M1/45‑2 (MQ=255)
                                          tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                            >  2:2221995/1‑101 (MQ=255)
                                                 ggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                                                     <  2:2452828‑M1/52‑1 (MQ=255)
                                                 cgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCg                                                     <  1:4618855‑M1/52‑1 (MQ=255)
                                                     gggtccgtgggatcggtcgggtccgtgggatctgtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTggg                                                 >  1:3692561‑M1/46‑101 (MQ=255)
                                                      gatccgtgggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGa                                                >  1:3170348‑M1/45‑101 (MQ=255)
                                                       gtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                               >  2:1882333‑M1/44‑101 (MQ=255)
                                                       gtccgtgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                                               >  1:3482350‑M1/44‑101 (MQ=255)
                                                        tccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                              >  2:3966380‑M1/43‑101 (MQ=255)
                                                        tccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                              >  1:867994‑M1/43‑101 (MQ=255)
                                                        tccgtgggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtc                                              >  2:1306129‑M1/43‑101 (MQ=255)
                                                         ccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                             <  2:904314‑M1/60‑1 (MQ=255)
                                                          cgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                            >  1:5299301‑M1/41‑101 (MQ=255)
                                                          cgtgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcgg                                            >  1:472465‑M1/41‑101 (MQ=255)
                                                           gtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                           <  1:5134552‑M1/62‑1 (MQ=255)
                                                            tgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                          <  1:2395505‑M1/63‑1 (MQ=255)
                                                            tgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                          >  1:1587261‑M1/39‑101 (MQ=255)
                                                             cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTAGGATCCGAGGGAtcggtcg                                         >  2:1513776/1‑101 (MQ=255)
                                                              ggatcggtgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                        <  1:4336705‑M1/65‑1 (MQ=255)
                                                               gatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCggg                                       <  1:2921368‑M1/66‑1 (MQ=255)
                                                               gatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGatcggtcgga                                       <  1:3385236‑M1/66‑2 (MQ=255)
                                                               gatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCggg                                       <  2:2708879‑M1/66‑1 (MQ=255)
                                                                    gtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                                  <  1:1682454‑M1/71‑1 (MQ=255)
                                                                     tgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTccggg                                 <  2:3395425‑M1/72‑3 (MQ=255)
                                                                      cggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                <  1:3779940‑M1/73‑1 (MQ=255)
                                                                      cggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCcgtcg                                <  2:1982293‑M1/73‑3 (MQ=255)
                                                                       ggggtgggtggATGCGTTGGGTCCGGTGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                               <  1:5695518/92‑1 (MQ=11)
                                                                         gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                             <  2:484663/101‑1 (MQ=255)
                                                                         gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                             <  2:4066746/101‑1 (MQ=255)
                                                                           ccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcg                           <  2:472415‑M1/78‑1 (MQ=255)
                                                                              tgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                        <  2:5027162‑M1/81‑1 (MQ=255)
                                                                              tCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccgtg                        <  1:3261705/101‑2 (MQ=11)
                                                                               gggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                       >  2:2931762‑M1/20‑101 (MQ=255)
                                                                               cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                       >  2:5494938/1‑101 (MQ=255)
                                                                                ggatcgggcggatccgtgGGATCGGTAGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                      >  1:605412‑M1/19‑101 (MQ=255)
                                                                                 gatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGa                     <  2:4385580‑M1/84‑1 (MQ=255)
                                                                                  aTGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                    >  1:5618691/1‑101 (MQ=255)
                                                                                       tgggagcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtc               <  2:3228432/95‑2 (MQ=11)
                                                                                       tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg               <  1:5498622/101‑1 (MQ=255)
                                                                                         ggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg             <  2:259139‑M1/92‑1 (MQ=255)
                                                                                         ggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTggg             >  1:1008891‑M1/10‑101 (MQ=255)
                                                                                           gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCgtggggt           >  1:2755997‑M1/8‑99 (MQ=255)
                                                                                            tccgtgggatctgtcgggtccgtgggatctgtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc          <  2:1704754/95‑1 (MQ=11)
                                                                                            tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGGCGGATCCGCGGGATCGGCTGGATCCGGGGGacc          >  2:5283964/7‑99 (MQ=11)
                                                                                            tccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGatg          <  2:3679055‑M1/95‑2 (MQ=255)
                                                                                             gcgtgGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATccgtggggtcc         >  2:792939‑M1/6‑97 (MQ=255)
                                                                                             cctttggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg         <  2:2804474/96‑1 (MQ=14)
                                                                                             ccgttggatcggtgggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg         <  2:4142636/96‑1 (MQ=11)
                                                                                             ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg         <  1:4412114‑M1/96‑1 (MQ=255)
                                                                                             ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg         <  2:3866065‑M1/96‑1 (MQ=255)
                                                                                              cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg        >  1:1921158‑M1/5‑101 (MQ=255)
                                                                                              cgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGggataga        <  1:404549‑M1/97‑4 (MQ=255)
                                                                                               gtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGatcggg       <  2:3017370‑M1/98‑2 (MQ=255)
                                                                                               gtcggTTCTGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt       <  2:4511478/101‑1 (MQ=11)
                                                                                               gtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAGCAGGCGGGTCCGTGGGAgcggt       >  2:3312357/1‑101 (MQ=11)
                                                                                                tgggatctgtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc      <  1:3096683/99‑1 (MQ=11)
                                                                                                tgggatcggtcgggtccgtcggatcggtgggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGggggacggcc      >  1:108819/3‑94 (MQ=11)
                                                                                                tgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc      <  2:4245642‑M1/99‑1 (MQ=255)
                                                                                                tcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGGTCCGTGGGATCGGTAGGGTCCGTGGGAtcggtc      >  1:2090636/1‑101 (MQ=255)
                                                                                                 gGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg     >  2:3119045‑M1/2‑101 (MQ=255)
                                                                                                 gGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg     <  1:1447801‑M1/100‑1 (MQ=255)
                                                                                                 cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg     >  1:1812201/1‑101 (MQ=14)
                                                                                                  ggatcggtcgggtccgtggggtccgtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGGGGGAtcggtcgg    >  2:66960/1‑101 (MQ=11)
                                                                                                   gatcggtcgggtccgtgggatcggtcgggtccgtcgGATCGGTGGGATCCGTTGGAGCGGTCGGGTCCGTGGGGTCGGTGGGGTCCGTCGGATCCGCCggg   >  1:3606653/1‑101 (MQ=11)
                                                                                                    atcggtcgggtccgtcggatcggtcgggtccgtgggATCGGACGGGTCCGTGGGCTCGGTCGGGTCCGGGGTATCGGTCGGGGGCGGGGGATCGGTCtgga  >  2:493925/1‑100 (MQ=11)
                                                                                                    atcggtcggatccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGat  >  1:5045279‑M2/1‑99 (MQ=255)
                                                                                                    |                                                                                                    
CGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221858‑222058

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