Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,446,527 C→T 11.8% A142A (GCC→GCT RCAP_rcc01343 → iron siderophore/cobalamin ABC transporter permease

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,446,5270C→T11.8% 48.6 ‑0.7 ‑0.0170A142A (GCC→GCTRCAP_rcc01343iron siderophore/cobalamin ABC transporter permease
Reads supporting (aligned to +/- strand):  new base (8/12):  ref base (77/73):  total (85/85)
Fisher's exact test strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.75e-01

Only 100 of 211 total aligned reads displayed.

TCTCGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGA  >  NC_014034/1446427‑1446627
                                                                                                    |                                                                                                    
tctcGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgcc                                                                                                      >  1:3577222/1‑101 (MQ=255)
tctcGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgcc                                                                                                      <  2:3344333/101‑1 (MQ=255)
tctcGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcccc                                                                                                      >  1:1934769/1‑101 (MQ=255)
 ctcGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccg                                                                                                     <  1:617120/101‑1 (MQ=255)
  tcGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgc                                                                                                    >  1:2949895/1‑101 (MQ=255)
    ggcggGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgcct                                                                                                  >  1:3279730/1‑101 (MQ=255)
    ggcggGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgcct                                                                                                  <  2:1087354/101‑1 (MQ=255)
      cggGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttc                                                                                                <  1:2629461/101‑1 (MQ=255)
      cggGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttc                                                                                                <  1:2506981/101‑1 (MQ=255)
       gggggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttcg                                                                                               >  1:4652800/1‑101 (MQ=255)
       gggggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttcg                                                                                               >  2:4251826/1‑101 (MQ=255)
        ggggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttcgc                                                                                              >  1:2413177/1‑101 (MQ=255)
          ggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCgcgc                                                                                            >  1:2467609/1‑101 (MQ=255)
           gCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCt                                                                                           <  1:3856418/101‑1 (MQ=255)
              tCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtggt                                                                                        <  2:3124852/101‑1 (MQ=255)
                 gccGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGgc                                                                                     <  2:569081/100‑1 (MQ=255)
                 aCCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGgc                                                                                     <  1:2846044/101‑1 (MQ=255)
                  ccGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGgcg                                                                                    >  1:120228/1‑101 (MQ=255)
                       gcgGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCtttt                                                                               >  1:4418154/1‑101 (MQ=255)
                        cgGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTg                                                                              >  2:1808185/1‑101 (MQ=255)
                         ggTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGt                                                                             >  1:2228778/1‑101 (MQ=255)
                          gTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTc                                                                            >  2:3667880/1‑101 (MQ=255)
                            ggCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGa                                                                          <  2:3529693/101‑1 (MQ=255)
                            ggCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGa                                                                          <  2:1564367/101‑1 (MQ=255)
                            ggCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGa                                                                          >  2:2002715/1‑101 (MQ=255)
                             gCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGaa                                                                         <  1:4134287/101‑1 (MQ=255)
                              cGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGacc                                                                        >  1:2303607/1‑99 (MQ=255)
                              cGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAAc                                                                        <  1:3089959/101‑1 (MQ=255)
                              cGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAAc                                                                        <  1:4734290/101‑1 (MQ=255)
                                gTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACgg                                                                      <  2:2938051/101‑1 (MQ=255)
                                  cATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGgcg                                                                    >  1:168417/1‑101 (MQ=255)
                                  cATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGgcg                                                                    <  1:2928576/101‑1 (MQ=255)
                                   aTCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGgcgc                                                                   >  2:4743725/1‑101 (MQ=255)
                                    cccttcTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGcgcg                                                                  <  1:605442/96‑1 (MQ=255)
                                      gTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAc                                                                >  1:3390585/1‑101 (MQ=255)
                                      gTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAc                                                                >  2:4226832/1‑101 (MQ=255)
                                      gTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAc                                                                >  2:3981806/1‑101 (MQ=255)
                                       tGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAcc                                                               <  2:3602886/101‑1 (MQ=255)
                                        gCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCa                                                              <  2:2200687/101‑1 (MQ=255)
                                          tGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGc                                                            <  1:233963/101‑1 (MQ=255)
                                          tGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGc                                                            >  1:1246567/1‑101 (MQ=255)
                                              cATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGc                                                        >  1:3886690/1‑101 (MQ=255)
                                              cATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGc                                                        >  1:2756514/1‑101 (MQ=255)
                                                tCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCg                                                      <  2:2555056/101‑1 (MQ=255)
                                                 cGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGa                                                     >  2:2955257/1‑101 (MQ=255)
                                                    cgcgGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAAcc                                                  <  1:1808294/101‑1 (MQ=40)
                                                    cgcgGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAAcc                                                  <  2:209142/101‑1 (MQ=40)
                                                     ccgggcgggCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGATGGTGGCGCTTTCGTCGAACGGCGCGACCAGCGGGCCGAACca                                                 >  1:313233/2‑101 (MQ=2)
                                                      cgggcgggCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcac                                                <  1:2264704/101‑1 (MQ=38)
                                                      cgggcgggCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcac                                                <  2:3818148/101‑1 (MQ=38)
                                                       gggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcaca                                               >  1:972064/1‑101 (MQ=38)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              >  2:2473740/1‑101 (MQ=38)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              <  1:3367013/101‑1 (MQ=38)
                                                         gcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACg                                             <  2:3225236/101‑1 (MQ=14)
                                                           gggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGAt                                           <  2:4055614/101‑1 (MQ=38)
                                                            ggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATc                                          <  2:2779149/101‑1 (MQ=14)
                                                            ggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATc                                          >  2:3011118/1‑101 (MQ=38)
                                                            ggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATc                                          >  2:3572384/1‑101 (MQ=38)
                                                            ggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATc                                          <  2:883902/101‑1 (MQ=38)
                                                             gctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATcc                                         <  1:170730/101‑1 (MQ=14)
                                                             gctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATca                                         <  2:1069375/101‑2 (MQ=14)
                                                              gtgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCt                                        <  2:4358658/100‑1 (MQ=14)
                                                              ctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCt                                        >  2:2140591/1‑101 (MQ=38)
                                                              ctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCt                                        <  1:1533169/101‑1 (MQ=38)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       <  2:3377266/101‑1 (MQ=14)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       <  2:2548117/101‑1 (MQ=38)
                                                                gtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcg                                      <  2:4127322/101‑1 (MQ=14)
                                                                gtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcg                                      <  2:757880/101‑1 (MQ=38)
                                                                gtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcg                                      >  2:1041394/1‑101 (MQ=38)
                                                                gtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcg                                      <  2:4116770/101‑1 (MQ=38)
                                                                gtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcg                                      <  1:4252073/101‑1 (MQ=38)
                                                                 tcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcgc                                     <  2:1427697/101‑1 (MQ=14)
                                                                 tcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcgc                                     >  1:3148945/1‑101 (MQ=14)
                                                                   gctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCgg                                   >  2:1897007/1‑101 (MQ=14)
                                                                    ctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCggg                                  <  1:4426128/101‑1 (MQ=11)
                                                                    ctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCggg                                  >  1:1490174/1‑101 (MQ=37)
                                                                    ctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCggg                                  >  1:1061988/1‑101 (MQ=37)
                                                                     tgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCgggg                                 >  2:4440038/1‑101 (MQ=35)
                                                                      gtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCggggg                                <  2:3443566/101‑1 (MQ=34)
                                                                       tcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCgggggc                               >  2:1901984/1‑100 (MQ=18)
                                                                       tcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGt                               <  1:219141/101‑1 (MQ=33)
                                                                        ccATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTc                              <  2:4517660/101‑1 (MQ=32)
                                                                           tGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgcc                           >  1:1408856/1‑101 (MQ=32)
                                                                           tGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgcc                           >  1:4017006/1‑101 (MQ=32)
                                                                            gggggCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGTCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgccg                          >  1:4506227/1‑101 (MQ=18)
                                                                            gggggCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgccg                          >  2:543856/1‑101 (MQ=33)
                                                                                     cgcgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGc                 >  1:567618/1‑101 (MQ=33)
                                                                                     cgcgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGc                 >  1:756585/1‑101 (MQ=33)
                                                                                      gcgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCATACGATCCTCGCGGGGGTCGCCGCCTCGCAGCt                >  2:1359225/1‑101 (MQ=17)
                                                                                      gcgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCt                <  1:4548685/101‑1 (MQ=32)
                                                                                       cgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCtt               <  2:4396691/101‑1 (MQ=32)
                                                                                        gggCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCttt              >  2:1712751/1‑101 (MQ=32)
                                                                                           cgccttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCa           <  1:2067946/101‑1 (MQ=32)
                                                                                            gccttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCaa          >  1:37018/1‑101 (MQ=32)
                                                                                             ccttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTTGCCGCCTCGCAGCTTTTCAAc         >  2:3732072/1‑101 (MQ=17)
                                                                                              cttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcg        >  2:1754853/1‑101 (MQ=32)
                                                                                                tcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgcg      <  2:4652544/101‑1 (MQ=32)
                                                                                                   ccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTg   <  1:82721/101‑1 (MQ=31)
                                                                                                    cgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGa  <  2:2654478/101‑1 (MQ=11)
                                                                                                    cgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGa  <  2:3279529/101‑1 (MQ=11)
                                                                                                    |                                                                                                    
TCTCGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGA  >  NC_014034/1446427‑1446627

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG