Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,786,541 C→T 6.0% A225A (GCC→GCT RCAP_rcc01648 → ABC transporter permease

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,786,5410C→T 6.0% 26.4 ‑0.6 ‑0.1217A225A (GCC→GCTRCAP_rcc01648ABC transporter permease
Reads supporting (aligned to +/- strand):  new base (4/9):  ref base (87/117):  total (91/126)
Fisher's exact test strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.76e-01

Only 100 of 229 total aligned reads displayed.

TGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGCGTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCT.GGCCGCTCGAACGGGATCGGCCTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGATCC  >  NC_014034/1786444‑1786641
                                                                                                  |                                                                                                    
tGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTtt                                                                                                    <  1:2586514/101‑1 (MQ=255)
tGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTtt                                                                                                    >  2:2554746/1‑101 (MQ=255)
tGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTtt                                                                                                    <  1:1452474/101‑1 (MQ=255)
tGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTtt                                                                                                    >  1:2756415/1‑101 (MQ=255)
 gATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTttc                                                                                                   <  1:796166/101‑2 (MQ=255)
 gATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTttc                                                                                                   <  2:227058/101‑2 (MQ=255)
 gATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGGCGAATTTCTGGGCCGATCGAACGGGATCGGCTttc                                                                                                   <  2:2308356/101‑2 (MQ=255)
  aTGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCg                                                                                                 <  2:2820813/101‑1 (MQ=255)
  aTGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCg                                                                                                 >  2:4378216/1‑101 (MQ=255)
  aTGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCg                                                                                                 <  1:393980/101‑1 (MQ=255)
   tGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGt                                                                                                <  1:875822/101‑1 (MQ=255)
   tGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGt                                                                                                <  2:3243508/101‑1 (MQ=255)
      cGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATc                                                                                             <  1:2061616/101‑1 (MQ=255)
       gCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATcc                                                                                            <  1:3450019/101‑1 (MQ=255)
         gtCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCAt                                                                                          <  1:3031805/100‑1 (MQ=255)
         aTCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCAt                                                                                          <  1:3688177/101‑1 (MQ=255)
           cTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCt                                                                                        <  1:1489474/101‑1 (MQ=255)
           cTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCt                                                                                        <  1:446300/101‑1 (MQ=255)
            cgATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTg                                                                                       <  1:1774206/100‑1 (MQ=255)
             gATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGa                                                                                      >  1:475134/1‑101 (MQ=255)
             gATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGa                                                                                      <  2:2799584/101‑1 (MQ=255)
             gATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGa                                                                                      <  1:101710/101‑1 (MQ=255)
             actCGCCGCGGCGCGGGCGGGGCTTTCGGTGATCCGGAACGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGa                                                                                      <  2:1117475/99‑1 (MQ=255)
               tCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAAt                                                                                    <  1:468016/101‑1 (MQ=255)
                cgcctcGGCGCGGGCGGGGCTTTCGCTGATCTTGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAAtt                                                                                   <  2:2237723/101‑1 (MQ=255)
                   cgcggcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTcc                                                                                <  2:4050267/101‑1 (MQ=255)
                   cgcggcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTcc                                                                                <  2:4278120/101‑1 (MQ=255)
                   cgcggcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTcc                                                                                >  1:1216800/1‑101 (MQ=255)
                    gcggcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCa                                                                               >  1:3533245/1‑101 (MQ=255)
                     cggcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAg                                                                              >  2:2090131/1‑101 (MQ=255)
                     cggcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAg                                                                              >  1:3840092/1‑101 (MQ=255)
                      ggcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGa                                                                             >  2:1474888/1‑101 (MQ=255)
                       gcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGAt                                                                            <  2:576839/101‑1 (MQ=255)
                         gcgggcggGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGt                                                                          <  1:2860005/101‑1 (MQ=255)
                          cgggcgggGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGtt                                                                         <  2:1689981/101‑1 (MQ=255)
                          cgggcgggGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGtt                                                                         <  2:4254752/101‑1 (MQ=255)
                          cgggcgggGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGtt                                                                         <  1:4415196/101‑1 (MQ=255)
                           gggcggggCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTc                                                                        >  2:611100/1‑101 (MQ=255)
                           gggcggggCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTc                                                                        >  2:1626298/1‑101 (MQ=255)
                           gggcggggCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTc                                                                        >  1:2475945/1‑101 (MQ=255)
                            ggcggggcTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCg                                                                       <  2:478542/101‑1 (MQ=255)
                              cggggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGAc                                                                     >  2:1646040/7‑101 (MQ=255)
                              cggggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGAc                                                                     <  2:1496053/95‑1 (MQ=255)
                               ggggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACa                                                                    >  1:3945183/6‑101 (MQ=255)
                               ggggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACa                                                                    <  1:3728771/96‑1 (MQ=255)
                                gggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAt                                                                   >  2:362915/5‑101 (MQ=255)
                                gggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAt                                                                   >  2:4429858/5‑101 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACTtt                                                                  <  1:2711284/98‑1 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  <  1:4595032/98‑1 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  >  1:1630334/4‑101 (MQ=255)
                                  gttttCGCTGGTTTTGAATGTTGTACTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.TTTTCGTATCAATCTGAATTTCCAGATGTTCGACATTg                                                                 <  2:2163953/99‑1 (MQ=255)
                                   ctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGc                                                                <  2:4333689/100‑1 (MQ=255)
                                     ttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCg                                                              <  2:2174438/101‑1 (MQ=255)
                                       cGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAg                                                            <  2:2729589/101‑1 (MQ=255)
                                       cGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAg                                                            <  1:4265333/101‑1 (MQ=255)
                                        gCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAgg                                                           >  2:3686278/1‑101 (MQ=255)
                                              cTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCg                                                     <  2:4345485/101‑1 (MQ=255)
                                              cTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCg                                                     <  2:347458/101‑1 (MQ=255)
                                               tGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGc                                                    <  2:2872557/101‑1 (MQ=255)
                                                 gAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCt                                                  >  1:65833/1‑101 (MQ=255)
                                                  aaGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTg                                                 >  2:3151158/1‑101 (MQ=255)
                                                  aaGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTg                                                 <  1:4104552/101‑1 (MQ=255)
                                                  aaGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTCCCGAGGTGCTCGCCTg                                                 <  1:4145943/101‑1 (MQ=255)
                                                    ggTCGTGGTGGTGGTCGAATTTCTGGGCCGCCCGAACGGGATCGG.CTTTAGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTggg                                               <  1:2092726/101‑1 (MQ=255)
                                                     gTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggc                                              <  2:1966925/101‑1 (MQ=255)
                                                     gTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggc                                              <  2:3850111/101‑1 (MQ=255)
                                                       cGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggcgg                                            <  1:3979597/101‑1 (MQ=255)
                                                        gTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGGGCTCGCCTGggcggc                                           >  2:4637788/1‑101 (MQ=255)
                                                         tGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggcggcg                                          >  1:2232475/1‑101 (MQ=255)
                                                           cTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggcggcggc                                        >  2:1834161/1‑101 (MQ=255)
                                                            tggtggTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCTGCt                                       >  2:3168294/1‑101 (MQ=255)
                                                            tggtggTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCt                                       <  2:414057/101‑1 (MQ=255)
                                                               tggtCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTc                                    >  2:2404766/1‑101 (MQ=255)
                                                                ggtCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCg                                   <  2:4496158/101‑1 (MQ=255)
                                                                ggtCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCg                                   >  1:1346589/1‑101 (MQ=255)
                                                                       tttCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTGGCCCGGGGGGGGGCgttttttttctg                            >  2:3942711/1‑92 (MQ=255)
                                                                       tttCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGt                            >  2:729363/1‑101 (MQ=255)
                                                                         tctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGa                          <  1:1999402/97‑1 (MQ=255)
                                                                          ctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGAt                         >  2:2040650/4‑101 (MQ=255)
                                                                          ctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGAt                         <  1:611164/98‑1 (MQ=255)
                                                                           tgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATg                        <  2:1748098/99‑1 (MQ=255)
                                                                           tgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATg                        <  1:2554898/99‑1 (MQ=255)
                                                                            gggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGc                       >  1:3386691/2‑101 (MQ=255)
                                                                             ggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCt                      >  1:1693278/1‑101 (MQ=255)
                                                                              gCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTc                     >  2:4612094/1‑101 (MQ=255)
                                                                              gCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTc                     <  2:1555230/101‑1 (MQ=255)
                                                                                    cGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTc               >  1:2271719/1‑101 (MQ=255)
                                                                                    cGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCCCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTc               <  1:3151358/101‑1 (MQ=255)
                                                                                     gAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCg              >  1:231799/1‑101 (MQ=255)
                                                                                     gAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCg              >  1:2404956/1‑101 (MQ=255)
                                                                                     gAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCg              >  2:3330183/1‑101 (MQ=255)
                                                                                     gAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCg              <  2:2658880/101‑1 (MQ=255)
                                                                                        cGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAAt           <  1:1974259/101‑1 (MQ=255)
                                                                                        cGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAAt           <  1:1187408/101‑1 (MQ=255)
                                                                                          ggATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCACGCCGTCGAATgg         >  2:1022572/1‑101 (MQ=255)
                                                                                              tcggcTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTg      <  2:2140959/97‑1 (MQ=255)
                                                                                              tcggcTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTg      <  1:113979/97‑1 (MQ=255)
                                                                                                cgcTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGAt    <  2:1371141/99‑1 (MQ=255)
                                                                                                 gcTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGATc   <  2:1236484/100‑1 (MQ=255)
                                                                                                  cTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGATcc  >  1:4550477/1‑101 (MQ=255)
                                                                                                  |                                                                                                    
TGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGCGTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCT.GGCCGCTCGAACGGGATCGGCCTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGATCC  >  NC_014034/1786444‑1786641

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG