Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,786,541 C→T 6.2% 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.2% 50.6 ‑0.6 ‑0.0389A225A (GCC→GCTRCAP_rcc01648ABC transporter permease
Reads supporting (aligned to +/- strand):  new base (10/14):  ref base (182/183):  total (192/197)
Fisher's exact test strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.79e-01

Only 100 of 420 total aligned reads displayed.

TGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGCGTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCT.GGCCGCTCGAACGGGATCGGCCTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGATCC  >  NC_014034/1786444‑1786641
                                                                                                  |                                                                                                    
tGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTtt                                                                                                    >  2:150163/1‑101 (MQ=255)
tGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTtt                                                                                                    >  2:8606835/1‑101 (MQ=255)
 gATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTttc                                                                                                   <  2:7557128/101‑2 (MQ=255)
 gATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTttc                                                                                                   <  2:4477730/101‑2 (MQ=255)
 gATGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGGCTttc                                                                                                   <  1:2389258/101‑2 (MQ=255)
   tGCCGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGt                                                                                                <  1:3112714/101‑1 (MQ=255)
     ccccATCGGATCGCCGCGGGCCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTAt                                                                                              <  2:4231093/98‑1 (MQ=255)
     ccGCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTAt                                                                                              >  1:751465/1‑101 (MQ=255)
       gCATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATcc                                                                                            <  2:3908477/101‑1 (MQ=255)
        cATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCa                                                                                           <  1:9173040/101‑1 (MQ=255)
        cATCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCa                                                                                           >  1:3441147/1‑101 (MQ=255)
          tCTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATc                                                                                         <  2:6562969/101‑1 (MQ=255)
           cTGATCGCCGCGGCGCGGGTGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCt                                                                                        >  1:6481741/1‑101 (MQ=255)
           cTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCt                                                                                        >  1:6966654/1‑101 (MQ=255)
           cTGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCt                                                                                        <  2:5126476/101‑1 (MQ=255)
           cTGATCGCCGCGGCGCGGGCGGGGCGTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCt                                                                                        <  1:3895189/101‑1 (MQ=255)
            tGATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTCTCCCTCTg                                                                                       >  2:533679/1‑101 (MQ=255)
             gATCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGa                                                                                      <  2:2228899/101‑1 (MQ=255)
              aTCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGaa                                                                                     <  1:6693691/101‑1 (MQ=255)
              aTCGCCGCGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGaa                                                                                     >  1:4827691/1‑101 (MQ=255)
                  ccgcGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTc                                                                                 <  2:6408258/101‑1 (MQ=255)
                  ccgcGGCGCGGGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTc                                                                                 <  1:561278/101‑1 (MQ=255)
                    gcggcgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCa                                                                               >  1:6058661/1‑101 (MQ=255)
                        cgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATg                                                                           <  2:3207395/101‑1 (MQ=255)
                        cgcggGCGGGGCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATg                                                                           <  1:474345/101‑1 (MQ=255)
                           gggcggggCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTc                                                                        >  2:4310509/1‑101 (MQ=255)
                           gggcggggCTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTc                                                                        <  2:8905524/101‑1 (MQ=255)
                              cggggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGAc                                                                     >  2:5068084/7‑101 (MQ=255)
                              cggggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGAc                                                                     >  2:1599698/7‑101 (MQ=255)
                              cggggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGAc                                                                     <  1:620988/95‑1 (MQ=255)
                               ggggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACa                                                                    >  1:8747350/6‑101 (MQ=255)
                                gggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAt                                                                   >  2:5578097/5‑101 (MQ=255)
                                gggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAt                                                                   <  2:6957144/97‑1 (MQ=255)
                                gggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAt                                                                   >  2:4833065/5‑101 (MQ=255)
                                 gtctttCGCGGATCTGGAAGGTCGTTCTGGTGTTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  <  2:2928446/98‑1 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  <  1:8607221/98‑1 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  <  2:4125102/98‑1 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  <  2:9152273/98‑1 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  <  2:7438284/98‑1 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  <  2:5841248/98‑1 (MQ=255)
                                 ggctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACAtt                                                                  <  1:7263633/98‑1 (MQ=255)
                                  gctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTg                                                                 <  1:2649531/99‑1 (MQ=255)
                                  gctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTg                                                                 <  1:3662376/99‑1 (MQ=255)
                                  gctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTg                                                                 <  1:7117738/99‑1 (MQ=255)
                                   ctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGc                                                                <  1:7399208/100‑1 (MQ=255)
                                   ctttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGc                                                                >  1:6326640/2‑101 (MQ=255)
                                    tttCGCTGATCTGGAAGGTCGTGCTGGTGTTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGcc                                                               <  1:4205868/101‑1 (MQ=255)
                                     ttCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCg                                                              >  1:8918638/1‑101 (MQ=255)
                                       cGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAg                                                            <  2:8748588/101‑1 (MQ=255)
                                        gCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAgg                                                           <  1:137394/101‑1 (MQ=255)
                                          tGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAGCGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTg                                                         >  2:7821958/1‑101 (MQ=255)
                                          tGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTg                                                         >  1:5204256/1‑101 (MQ=255)
                                          tGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTg                                                         >  2:7660311/1‑101 (MQ=255)
                                            aTCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCt                                                       <  2:6299422/101‑1 (MQ=255)
                                             tCTGGAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTc                                                      <  2:3054053/101‑1 (MQ=255)
                                                ggAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGGGCTCGcc                                                   >  2:1065260/1‑101 (MQ=255)
                                                 gAAGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCt                                                  >  1:1755547/1‑101 (MQ=255)
                                                  aaGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTg                                                 >  2:5721941/1‑101 (MQ=255)
                                                   aGGTCGTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTgg                                                >  2:6164830/1‑101 (MQ=255)
                                                        gTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggcggc                                           >  2:4954413/1‑101 (MQ=255)
                                                        gTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggcggc                                           >  1:4583730/1‑101 (MQ=255)
                                                        gTGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggcggc                                           >  1:4661404/1‑101 (MQ=255)
                                                         tGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggcggcg                                          >  2:1747943/1‑101 (MQ=255)
                                                         tGCTGGTGGTCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGggcggcg                                          <  1:150163/101‑1 (MQ=255)
                                                             ggtggtCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCtt                                      <  2:1389953/101‑1 (MQ=255)
                                                               tggtCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTc                                    <  2:2654367/101‑1 (MQ=255)
                                                                ggtCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCg                                   <  1:3777372/101‑1 (MQ=255)
                                                                 gtCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGt                                  >  2:676210/1‑101 (MQ=255)
                                                                  tCGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTc                                 >  2:7980353/1‑101 (MQ=255)
                                                                   cGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCt                                >  1:2235335/1‑101 (MQ=255)
                                                                   cGAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCt                                >  2:8198823/1‑101 (MQ=255)
                                                                    gAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTg                               >  1:7994878/1‑101 (MQ=255)
                                                                    gAATTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTg                               >  2:2292040/1‑101 (MQ=255)
                                                                     aaTTTCTGGGCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTGCGTCTGc                              >  2:6236713/1‑101 (MQ=255)
                                                                        ttctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTg                           >  2:1284155/6‑101 (MQ=255)
                                                                        ttctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTg                           >  1:3021472/6‑101 (MQ=255)
                                                                        ttctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTg                           <  2:2514043/96‑1 (MQ=255)
                                                                        ttctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTg                           >  2:8426230/6‑101 (MQ=255)
                                                                        ttctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTg                           <  1:1911577/96‑1 (MQ=255)
                                                                         tctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGa                          >  2:2008572/5‑101 (MQ=255)
                                                                          ctgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGAt                         >  1:1065242/4‑101 (MQ=255)
                                                                           tgggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGAGATg                        >  1:2825765/3‑101 (MQ=255)
                                                                            gggCCGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGGTGc                       >  2:7326098/2‑101 (MQ=255)
                                                                               gcGCTGGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGGGGCTTTCGTCTGCGTGATGCTCg                    <  2:3719028/100‑1 (MQ=255)
                                                                                cGCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGc                   <  1:1651099/101‑1 (MQ=255)
                                                                                 gCTCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGcc                  <  1:4430047/101‑1 (MQ=255)
                                                                                   tCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGt                >  1:5511777/1‑101 (MQ=255)
                                                                                   tCGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGt                >  1:5795758/1‑101 (MQ=255)
                                                                                    cGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTc               >  1:127723/1‑101 (MQ=255)
                                                                                    cGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTc               >  1:3547243/1‑101 (MQ=255)
                                                                                    cGAACGGGATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGAGGCGGCTTTCGTCTGCGTGATGCTCGCCGTc               >  1:8148601/1‑101 (MQ=255)
                                                                                         gggATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATg          <  2:4661255/101‑1 (MQ=255)
                                                                                         gggATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATg          <  2:3341873/101‑1 (MQ=255)
                                                                                          ggATCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATgg         >  1:5115009/1‑101 (MQ=255)
                                                                                            aTCGG.CTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCt       >  1:3259514/1‑101 (MQ=255)
                                                                                              tcggcTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAAtagatc      >  2:3441378/5‑96 (MQ=255)
                                                                                              tcggcTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTg      <  1:7816905/97‑1 (MQ=255)
                                                                                               cggcTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGa     >  1:8290208/4‑101 (MQ=255)
                                                                                                ggcTTTTGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGAt    >  2:1839551/3‑101 (MQ=255)
                                                                                                  cTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGATcc  >  2:5950399/1‑101 (MQ=255)
                                                                                                  |                                                                                                    
TGATGCCGCATCTGATCGCCGCGGCGCGGGCGGGCGTTTCGCTGATCTGGAAGGTCGTGCTGGTGGTCGAATTTCT.GGCCGCTCGAACGGGATCGGCCTTTCGTATCCATCTGAATTTCCAGATGTTCGACATTGCCGAGGTGCTCGCCTGGGCGGCGGCTTTCGTCTGCGTGATGCTCGCCGTCGAATGGCTGATCC  >  NC_014034/1786444‑1786641

Base quality scores:  ATCG  < 3 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG