Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,446,527 C→T 10.5% 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→T10.5% 48.1 ‑0.2 ‑0.3181A142A (GCC→GCT) RCAP_rcc01343iron siderophore/cobalamin ABC transporter permease
Reads supporting (aligned to +/- strand):  new base (10/9):  ref base (91/71):  total (101/80)
Fisher's exact test strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.50e-01

Only 100 of 223 total aligned reads displayed.

TCTCGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGA  >  NC_014034/1446427‑1446627
                                                                                                    |                                                                                                    
 ctcGGCGGGGGCCTCGACCGGCGCGGGGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccg                                                                                                     >  2:3246688/1‑101 (MQ=255)
      cggGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttc                                                                                                >  2:583476/1‑101 (MQ=255)
       gggggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttcg                                                                                               >  1:3914356/1‑101 (MQ=255)
       gggggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttcg                                                                                               >  2:1003864/1‑101 (MQ=255)
       gggggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttcg                                                                                               <  2:1455035/101‑1 (MQ=255)
            ccTCGACCGGCGCGGTGGCGGTCATCGTGCTGGTCATCGGCGCGGGCGGGCTGTCGCCGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtg                                                                                          <  1:3806899/101‑1 (MQ=255)
            ccTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtg                                                                                          >  2:741363/1‑101 (MQ=255)
            ccTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtg                                                                                          <  2:3914102/101‑1 (MQ=255)
            ccTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtg                                                                                          <  2:3801652/101‑1 (MQ=255)
             cTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtgg                                                                                         <  2:917845/101‑1 (MQ=255)
             cTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtgg                                                                                         >  1:4236626/1‑101 (MQ=255)
              tCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtggt                                                                                        <  2:240554/101‑1 (MQ=255)
              tCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtggt                                                                                        <  2:2699760/101‑1 (MQ=255)
                   cGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGgcgc                                                                                   >  1:2515391/1‑101 (MQ=255)
                    ggCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCt                                                                                  >  1:839702/1‑101 (MQ=255)
                    ggCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGCCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCt                                                                                  <  1:3103069/101‑1 (MQ=255)
                       gcgGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCtttt                                                                               >  2:3469460/1‑101 (MQ=255)
                        cgGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTg                                                                              >  2:4266303/1‑101 (MQ=255)
                        cgGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTg                                                                              <  1:2197991/101‑1 (MQ=255)
                         ggTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGt                                                                             >  1:1481964/1‑101 (MQ=255)
                           tGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCg                                                                           >  2:1788248/1‑101 (MQ=255)
                           tGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCg                                                                           >  1:4112215/1‑101 (MQ=255)
                           tGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCg                                                                           >  2:489948/1‑101 (MQ=255)
                           gggCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCg                                                                           <  2:4239826/100‑1 (MQ=255)
                             gCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGaa                                                                         <  2:2323006/101‑1 (MQ=255)
                              cGGGCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAAc                                                                        <  2:1992719/101‑1 (MQ=255)
                                 tCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGGGCCTTCGCCGGCTTCGCGCTGGGGGCGCTTTTGTCGAACGgc                                                                     >  2:4145443/1‑101 (MQ=255)
                                  cATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGgcg                                                                    <  1:1295961/101‑1 (MQ=255)
                                  cATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGgcg                                                                    <  1:3556884/101‑1 (MQ=255)
                                      gTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAc                                                                >  1:1584416/1‑101 (MQ=255)
                                       tGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAcc                                                               >  1:3665067/1‑101 (MQ=255)
                                       tGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAcc                                                               <  1:1003949/101‑1 (MQ=255)
                                        gCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCa                                                              >  2:3085345/1‑101 (MQ=255)
                                        gCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCa                                                              <  2:2376901/101‑1 (MQ=255)
                                          tGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGc                                                            >  1:2880146/1‑101 (MQ=255)
                                          tGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGc                                                            >  2:2775126/1‑101 (MQ=255)
                                           gggCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCg                                                           <  2:2987876/101‑1 (MQ=255)
                                            ggCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCgg                                                          <  2:226549/101‑1 (MQ=255)
                                              cATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCTCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGc                                                        <  2:423529/101‑1 (MQ=255)
                                               aTCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGcc                                                       >  1:3336178/1‑101 (MQ=255)
                                               aTCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGcc                                                       >  1:3372662/1‑101 (MQ=255)
                                               aTCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGcc                                                       >  1:1950000/1‑101 (MQ=255)
                                                tCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCg                                                      >  1:2312291/1‑101 (MQ=255)
                                                 cGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGa                                                     <  1:1141888/101‑1 (MQ=255)
                                                    ggcgggcggGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAAcc                                                  <  2:3542335/100‑1 (MQ=39)
                                                    cgcgGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAAcc                                                  <  1:3290526/101‑1 (MQ=40)
                                                    cgcgGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAAcc                                                  >  2:1293978/1‑101 (MQ=40)
                                                     ccgggcgggCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACca                                                 <  1:2264323/100‑1 (MQ=11)
                                                      cgggcgggCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcac                                                >  1:3072631/1‑101 (MQ=14)
                                                       gggcgggcTGTCGCTGTCCATGGGGGCTTTGGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcaca                                               <  2:3133283/101‑1 (MQ=38)
                                                       gggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcaca                                               <  1:1842164/101‑1 (MQ=14)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              >  2:2308815/1‑101 (MQ=14)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              <  1:904212/101‑1 (MQ=38)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              <  2:2199946/101‑1 (MQ=38)
                                                         gcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACg                                             >  1:3120859/1‑101 (MQ=38)
                                                          cgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGa                                            <  1:126061/101‑1 (MQ=38)
                                                           gggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGAt                                           >  2:252041/1‑101 (MQ=38)
                                                           gggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGAt                                           >  1:1723752/1‑101 (MQ=38)
                                                            ggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATc                                          >  1:3000369/1‑101 (MQ=14)
                                                             gctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATcc                                         <  1:3436846/101‑1 (MQ=14)
                                                             gctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATcc                                         >  2:1930780/1‑101 (MQ=38)
                                                              ctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCt                                        <  2:2040488/101‑1 (MQ=14)
                                                              ctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCt                                        <  2:3751815/101‑1 (MQ=38)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       >  1:4032165/1‑101 (MQ=14)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       >  1:96249/1‑101 (MQ=14)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       <  2:2627663/101‑1 (MQ=38)
                                                                gtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcg                                      <  2:509924/101‑1 (MQ=14)
                                                                gtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcg                                      <  1:2459024/101‑1 (MQ=14)
                                                                 tcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcgc                                     >  2:4394023/1‑101 (MQ=38)
                                                                  cgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcgcg                                    <  2:4041394/101‑1 (MQ=38)
                                                                  cgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTcgcg                                    >  1:1961732/1‑101 (MQ=38)
                                                                   gctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCgg                                   >  1:1548026/1‑101 (MQ=38)
                                                                   gctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCgg                                   >  1:442965/1‑101 (MQ=38)
                                                                    ctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCggg                                  >  2:3636721/1‑101 (MQ=37)
                                                                      gtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCggggg                                >  1:2166192/1‑101 (MQ=34)
                                                                        ccATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTc                              <  2:3176142/101‑1 (MQ=32)
                                                                          aTGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgc                            <  2:3927074/101‑1 (MQ=32)
                                                                           tGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgcc                           <  1:2775316/101‑1 (MQ=32)
                                                                           tGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgcc                           <  1:762848/101‑1 (MQ=32)
                                                                           tGGGGGCTTTCGCGGGCCCCTTCGCCCCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgcc                           <  2:395828/101‑1 (MQ=12)
                                                                            gggggCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgccg                          >  1:263873/1‑101 (MQ=33)
                                                                               gcttttccggggCGCCTTCCCCGCCTCCCCGCTGGTGGCCCTTTTGTCGAACGGCGCG.CCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTc                      <  1:2687942/92‑1 (MQ=11)
                                                                                gCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTc                      <  1:528172/101‑1 (MQ=32)
                                                                                 cTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCg                     >  1:1450428/1‑101 (MQ=33)
                                                                                 cTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCg                     >  1:272828/1‑101 (MQ=33)
                                                                                  tttCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAAGGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGc                    >  2:725603/1‑101 (MQ=17)
                                                                                   ttCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCa                   <  1:562674/101‑1 (MQ=32)
                                                                                   ttCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCa                   >  1:4329180/1‑101 (MQ=32)
                                                                                       cgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCCGCtt               >  2:65827/1‑101 (MQ=17)
                                                                                        gggCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCttt              >  1:3482271/1‑101 (MQ=32)
                                                                                        gggCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCttt              >  1:4127984/1‑101 (MQ=32)
                                                                                           cgccttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCa           <  2:4122294/101‑1 (MQ=32)
                                                                                           cgccttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCa           <  1:672257/101‑1 (MQ=32)
                                                                                               ttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgc       >  1:3222614/1‑101 (MQ=32)
                                                                                               ttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgc       >  1:3717128/1‑101 (MQ=32)
                                                                                                tcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgcg      >  2:1569445/1‑101 (MQ=32)
                                                                                                tcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgcg      <  1:589732/101‑1 (MQ=32)
                                                                                                 cgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACgcgc     <  2:2534225/101‑1 (MQ=32)
                                                                                                  gccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCt    <  2:1547921/101‑1 (MQ=32)
                                                                                                    cgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGa  >  2:3266623/1‑101 (MQ=11)
                                                                                                    |                                                                                                    
TCTCGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGA  >  NC_014034/1446427‑1446627

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