Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,446,527 C→T 9.7% 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→T 9.7% 84.1 ‑1.6 ‑0.1350A142A (GCC→GCTRCAP_rcc01343iron siderophore/cobalamin ABC transporter permease
Reads supporting (aligned to +/- strand):  new base (11/23):  ref base (142/174):  total (153/197)
Fisher's exact test strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.19e-01

Only 100 of 438 total aligned reads displayed.

TCTCGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGA  >  NC_014034/1446427‑1446627
                                                                                                    |                                                                                                    
 ctcGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCCGTCCATGGGGGCTTTCGCGGGcgccttcgccg                                                                                                     >  1:4455593/1‑101 (MQ=255)
    ggcggGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgcct                                                                                                  <  1:1005234/101‑1 (MQ=255)
    ggcggGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgcct                                                                                                  >  1:8864860/1‑101 (MQ=255)
     gcggGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgcctt                                                                                                 >  1:8538746/1‑101 (MQ=255)
       gggggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttcg                                                                                               >  2:2280448/1‑101 (MQ=255)
       gggggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGcgccttcgccgccttcg                                                                                               >  2:5277824/1‑101 (MQ=255)
         gggCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTcgcg                                                                                             <  2:7898107/101‑1 (MQ=255)
           gCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCt                                                                                           >  1:9024636/1‑101 (MQ=255)
           gCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCt                                                                                           <  1:4080963/101‑1 (MQ=255)
           gCCGCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCt                                                                                           <  2:4211572/101‑1 (MQ=255)
            ccTCGACNGGCGCGNTGGCGGTGATCNNNNNGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCCCCGCCTTCGCGCtg                                                                                          <  1:6940499/101‑1 (MQ=14)
            ccTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtg                                                                                          >  2:6197781/1‑101 (MQ=255)
             cTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtgg                                                                                         >  2:172873/1‑101 (MQ=255)
             cTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtgg                                                                                         >  2:1939065/1‑101 (MQ=255)
               cGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtggtg                                                                                       >  1:8782993/1‑101 (MQ=255)
               cGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtggtg                                                                                       <  2:446026/101‑1 (MQ=255)
                gACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCtggtgg                                                                                      >  1:3766888/1‑101 (MQ=255)
                 aCCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGGTGGTGgc                                                                                     >  2:6237372/1‑101 (MQ=255)
                   cGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGgcgc                                                                                   >  2:3517445/1‑101 (MQ=255)
                       gcgGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCtttt                                                                               <  2:757286/101‑1 (MQ=255)
                         ggTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGt                                                                             >  1:2333452/1‑101 (MQ=255)
                          gTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTc                                                                            >  2:6903557/1‑101 (MQ=255)
                          gTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTc                                                                            <  2:1560305/101‑1 (MQ=255)
                           tGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCg                                                                           >  2:508760/1‑101 (MQ=255)
                             gCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGaa                                                                         <  2:3729995/101‑1 (MQ=255)
                                gTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACgg                                                                      <  2:219571/101‑1 (MQ=255)
                                gTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACgg                                                                      >  1:465181/1‑101 (MQ=255)
                                gTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACgg                                                                      >  1:658038/1‑101 (MQ=255)
                                gTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACgg                                                                      <  2:2749909/101‑1 (MQ=255)
                                   aTCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGgcgc                                                                   <  2:8046139/101‑1 (MQ=255)
                                    tCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGcgcg                                                                  <  1:734896/101‑1 (MQ=255)
                                     cGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGa                                                                 >  1:7843738/1‑101 (MQ=255)
                                     cGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGa                                                                 >  2:7730177/1‑101 (MQ=255)
                                      gTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAc                                                                >  2:499723/1‑101 (MQ=255)
                                       tGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAcc                                                               <  2:8538361/101‑1 (MQ=255)
                                       tGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGAcc                                                               <  1:391667/101‑1 (MQ=255)
                                       tGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACCGCGCGAcc                                                               >  2:7641972/1‑101 (MQ=255)
                                        gCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCa                                                              <  2:5429815/101‑1 (MQ=255)
                                         cTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAg                                                             <  2:6354040/101‑1 (MQ=255)
                                         cTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAg                                                             >  2:1256719/1‑101 (MQ=255)
                                         cTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAg                                                             >  1:1395197/1‑101 (MQ=255)
                                             gCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCggg                                                         >  2:2663160/1‑101 (MQ=255)
                                             gCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTATTGTCGAACGGCGCGACCAGCggg                                                         >  2:9063425/1‑101 (MQ=255)
                                              cATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGCCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGc                                                        <  2:5902147/101‑1 (MQ=255)
                                                tCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCg                                                      <  1:4169146/101‑1 (MQ=255)
                                                  ggCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGac                                                    >  2:2511694/1‑100 (MQ=255)
                                                  ggCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGaa                                                    <  1:2620545/101‑1 (MQ=255)
                                                      cgggcgggCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcac                                                <  1:6267925/101‑1 (MQ=38)
                                                      cgggcgggCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcac                                                <  1:6863336/101‑1 (MQ=38)
                                                      cgggcgggCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcac                                                >  1:6018848/1‑101 (MQ=38)
                                                       gggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACcaca                                               <  1:6365029/101‑1 (MQ=38)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              <  1:1548816/101‑1 (MQ=14)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              <  1:8271268/101‑1 (MQ=14)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              <  1:1736263/101‑1 (MQ=38)
                                                        ggcgggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCacac                                              <  1:1222675/101‑1 (MQ=38)
                                                           gggcTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGAt                                           <  2:8350048/101‑1 (MQ=14)
                                                             gctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATcc                                         <  2:1369472/101‑1 (MQ=38)
                                                             gctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATcc                                         <  2:3066234/101‑1 (MQ=38)
                                                              ctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCGGCGGGCCGAACCACACGATCCt                                        <  2:1018300/101‑1 (MQ=9)
                                                              ctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCt                                        >  1:2185620/1‑101 (MQ=14)
                                                              ctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCTCGACCAGCGGGCCGAACCACACGATCCt                                        >  1:3298717/1‑101 (MQ=38)
                                                              ctgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCt                                        <  2:6571084/101‑1 (MQ=38)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       <  1:2918390/101‑1 (MQ=14)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCTGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       <  2:1182058/101‑1 (MQ=14)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       <  1:3517552/101‑1 (MQ=38)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       <  1:6474/101‑1 (MQ=38)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTc                                       >  2:3285672/1‑101 (MQ=38)
                                                               tgtcgctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGATGGTGGCGCTTTTGTCGAACGGCGCGCCCAGCGGGCCGCACCACCCGAGcccc                                       >  2:8948582/1‑99 (MQ=14)
                                                                    ctgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCggg                                  >  2:4039228/1‑101 (MQ=37)
                                                                     tgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCgggg                                 >  1:4901260/1‑101 (MQ=35)
                                                                     tgtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCgggg                                 <  1:5218780/101‑1 (MQ=35)
                                                                      gtcCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCggggg                                >  2:869388/1‑101 (MQ=34)
                                                                         cATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcg                             <  2:7025334/101‑1 (MQ=32)
                                                                          aTGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgc                            >  1:3672441/1‑101 (MQ=32)
                                                                           tGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgcc                           <  1:2663242/101‑1 (MQ=32)
                                                                           tGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTcgcc                           >  1:7653096/1‑101 (MQ=32)
                                                                              gggCTTTCGCGGGGGCCTTCGCCGCCTCCGCGCTGGTGGCGCTTTTGTGGAACGGCGCGACCACCGGGCCGAACCACCCGATCCTCGCGGGGGTCgccgcc                        <  2:595438/101‑1 (MQ=11)
                                                                               ggCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCgccgccc                       <  1:6118100/101‑2 (MQ=17)
                                                                                 cTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCg                     <  1:8646637/101‑1 (MQ=32)
                                                                                  tttCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGc                    >  1:7948643/1‑101 (MQ=32)
                                                                                   tccgcgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCa                   <  2:4418059/99‑1 (MQ=17)
                                                                                     cgcgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGc                 <  1:706232/101‑1 (MQ=32)
                                                                                       cgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCtt               <  1:1466239/101‑1 (MQ=32)
                                                                                       cgGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCtt               >  1:5518716/1‑101 (MQ=33)
                                                                                        gggCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCttt              <  2:4828582/101‑1 (MQ=32)
                                                                                        gggCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCttt              >  1:2779429/1‑101 (MQ=32)
                                                                                         ggCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCtttt             <  2:3943928/101‑1 (MQ=32)
                                                                                          gcgcCTTCGCCGCCTTCGCGTTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTCTTc            <  1:7204703/101‑1 (MQ=12)
                                                                                            gccttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCaa          <  2:465173/101‑1 (MQ=32)
                                                                                            gccttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCaa          >  2:3198157/1‑101 (MQ=32)
                                                                                            gccttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGCGCCGAACCACACGATCCNNNNNCGGGTCGCCGCNTCNCAGCTTTTCaa          >  1:9137862/1‑101 (MQ=11)
                                                                                              gatctccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcg        <  1:8097955/96‑1 (MQ=14)
                                                                                              cttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcg        <  2:3590004/101‑1 (MQ=32)
                                                                                              cttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcg        >  1:7483448/1‑101 (MQ=32)
                                                                                               ttcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgc       <  2:7352753/101‑1 (MQ=32)
                                                                                                tcgccgccttcgcGCTGGTGGCGCTTTTGTCGTACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgcg      <  2:6761765/101‑1 (MQ=17)
                                                                                                tcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgcg      >  1:571045/1‑101 (MQ=32)
                                                                                                tcgccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAAcgcg      <  1:8407531/101‑1 (MQ=32)
                                                                                                  gccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCt    <  2:8162857/101‑1 (MQ=32)
                                                                                                   ccgccttcgcGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTg   >  2:7488753/1‑101 (MQ=31)
                                                                                                    |                                                                                                    
TCTCGGCGGGGGCCTCGACCGGCGCGGTGGCGGTCATCGTGCTGGGCATCGGCGCGGGCGGGCTGTCGCTGTCCATGGGGGCTTTCGCGGGCGCCTTCGCCGCCTTCGCGCTGGTGGCGCTTTTGTCGAACGGCGCGACCAGCGGGCCGAACCACACGATCCTCGCGGGGGTCGCCGCCTCGCAGCTTTTCAACGCGCTGA  >  NC_014034/1446427‑1446627

Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG