Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,563 C→G 18.5% P149A (CCC→GCC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,5630C→G18.5% 106.5 ‑0.9 ‑0.4227P149A (CCC→GCC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (19/23):  ref base (98/87):  total (117/110)
Fisher's exact test strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.09e-01

Only 100 of 252 total aligned reads displayed.

CGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCGTGGA  >  NC_014034/1396464‑1396663
                                                                                                     |                                                                                                    
cGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAgccgcc                                                                                                     >  1:3589218/1‑101 (MQ=255)
cGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAgccgcc                                                                                                     <  2:2114405/101‑1 (MQ=255)
   aaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGccgcccgc                                                                                                  <  2:2877478/101‑1 (MQ=255)
   aaCGGACCTGAAATTCCGCGGGAAAAGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGccgcccgc                                                                                                  >  2:567031/1‑101 (MQ=255)
      ggACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGcc                                                                                               >  2:2577149/1‑101 (MQ=255)
       gACCCGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCa                                                                                              <  2:1240579/101‑1 (MQ=255)
        aCCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCaa                                                                                             <  1:3174943/101‑1 (MQ=255)
             aaaTTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGc                                                                                        <  2:699844/101‑1 (MQ=255)
                  ccGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTcccc                                                                                   >  1:451277/1‑101 (MQ=255)
                     cgGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCcg                                                                                <  1:3940841/101‑1 (MQ=255)
                       ggAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCcgcg                                                                              >  2:3667426/1‑101 (MQ=255)
                        gAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCgg                                                                             >  2:2498898/1‑101 (MQ=255)
                         aaaCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCggg                                                                            >  1:338772/1‑101 (MQ=255)
                         aaaCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCggg                                                                            >  1:3715078/1‑101 (MQ=255)
                          aaCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCgggg                                                                           >  2:3370717/1‑101 (MQ=255)
                          aaCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCCCGCCCAGCCGCTCCCCCCGCgggg                                                                           >  2:3954264/1‑101 (MQ=255)
                            cgcgGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGc                                                                         <  2:1512149/101‑1 (MQ=255)
                              cggcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCgg                                                                       >  2:820662/1‑100 (MQ=255)
                               ggcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTg                                                                      <  1:1539591/101‑1 (MQ=255)
                                gcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTgg                                                                     <  1:1187025/101‑1 (MQ=255)
                                 ttgggccggagTCGCGCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTgg                                                                     <  1:4168031/91‑1 (MQ=2)
                                 ttgggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCCCCCCGCGGGGGCGCgg                                                                     >  2:3603160/11‑101 (MQ=2)
                                 ttgggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTgg                                                                     <  1:4227880/91‑1 (MQ=2)
                                  tgggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGa                                                                    >  1:1615433/10‑101 (MQ=2)
                                  tgggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGa                                                                    >  1:4365740/10‑101 (MQ=2)
                                  ggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   >  2:1652462/1‑101 (MQ=255)
                                   gggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGAt                                                                   >  2:946460/9‑101 (MQ=2)
                                    cgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 <  2:3239442/101‑1 (MQ=255)
                                     gacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCg                                                                 >  2:1155358/7‑101 (MQ=2)
                                     gAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                >  1:1590526/1‑101 (MQ=255)
                                     gAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                <  1:364332/101‑1 (MQ=255)
                                       gggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGaa                                                               <  2:3570191/97‑1 (MQ=2)
                                       gCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.GGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                              <  1:1374645/101‑1 (MQ=255)
                                       gCACGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGGGGATCGaaa                                                              >  1:735965/1‑101 (MQ=255)
                                       cggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGaa                                                               <  1:2480993/97‑1 (MQ=2)
                                        cccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                             <  1:3293299/99‑1 (MQ=255)
                                         gagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGGCAACCTATGTCCCTCACCCCGCGGGGGCGTGGATCGACAc                                                             >  2:1866597/3‑101 (MQ=2)
                                          ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAAct                                                            <  2:4385113/100‑2 (MQ=2)
                                          cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGa                                                           >  1:1749200/1‑101 (MQ=255)
                                          cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGa                                                           >  2:1991920/1‑101 (MQ=255)
                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt                                                          >  1:2379599/1‑101 (MQ=255)
                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActg                                                           <  2:1205994/101‑3 (MQ=2)
                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                                        >  2:1598266/1‑101 (MQ=255)
                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                                        <  1:2025649/101‑1 (MQ=255)
                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                                        >  1:2868377/1‑101 (MQ=255)
                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtg                                                         >  2:734698/1‑97 (MQ=2)
                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtg                                                         <  2:3831300/101‑5 (MQ=2)
                                             cGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                                        <  1:212933/101‑1 (MQ=255)
                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtc                                                       >  2:2259416/1‑101 (MQ=255)
                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtc                                                       <  2:1749081/101‑1 (MQ=255)
                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAG.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgc                                                        <  1:826396/101‑6 (MQ=2)
                                               cTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgcc                                                       <  1:3503333/101‑7 (MQ=2)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccg                                                      >  1:3123024/1‑94 (MQ=2)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccg                                                      >  1:280896/1‑94 (MQ=2)
                                                 ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTTGCCAGCCTCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGTGTCGAAACGATTtagtc                                                    >  2:1000415/1‑101 (MQ=37)
                                                 ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgtc                                                    <  1:661760/101‑1 (MQ=255)
                                                 ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgtc                                                    <  1:399589/101‑1 (MQ=255)
                                                 ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCATTtcggc                                                    >  2:1336621/1‑99 (MQ=255)
                                                 ccccccGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccgg                                                     >  1:861629/1‑93 (MQ=2)
                                                 ccccccGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccgg                                                     >  1:123665/1‑93 (MQ=2)
                                                 ccccccGAGGGGGCGTGGATCGAAACAGCTTTC.GATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccgg                                                     <  1:1042252/101‑9 (MQ=2)
                                                  cccccGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccggt                                                    >  1:19679/1‑92 (MQ=2)
                                                  cccccGCGGGGGCGTGGATCGAAACAGCTTTCAA.TGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccggt                                                    <  1:3050161/101‑10 (MQ=2)
                                                   ccccGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGc                                                  >  1:522058/1‑101 (MQ=255)
                                                    cccGCGGGGGCGTCGATCGAAACTGCTTTC.G.TGGGCCTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCt                                                 <  2:1417937/101‑1 (MQ=255)
                                                     ccGGGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCtt                                                <  2:4274068/101‑1 (MQ=255)
                                                      cgcgGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTa                                               <  1:3415596/101‑1 (MQ=255)
                                                      cgcgGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTa                                               <  2:2676829/101‑1 (MQ=255)
                                                       gcgGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTac                                              <  1:2695631/101‑1 (MQ=255)
                                                       gcgGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTac                                              >  1:974813/1‑101 (MQ=255)
                                                        ggggggCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacg                                             <  1:2772321/100‑1 (MQ=255)
                                                           gggCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacgacc                                          >  2:2365714/1‑100 (MQ=255)
                                                            ggCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGc                                         >  1:2157769/1‑101 (MQ=255)
                                                              cGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCCTGGATCGAAACGATTTCGTCGCTTACGACGCCt                                       >  1:3843048/1‑101 (MQ=255)
                                                               gTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTa                                      >  2:272369/1‑101 (MQ=255)
                                                                 ggATCGAAACTGCTTTC.G.TGGGACTGGCCAGGCGCCCGCGCCAAGCCGCTCCCCCCGCGGGGGCGGGGGTCCGAACCACTTCGTCGGTTCCGCCGGCTATc                                    >  2:546779/1‑101 (MQ=255)
                                                                  gATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCt                                   <  1:3542329/101‑1 (MQ=255)
                                                                    tCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGa                                 <  2:451237/101‑1 (MQ=255)
                                                                     cGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAt                                >  1:3913835/1‑101 (MQ=255)
                                                                       aaaCTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGACACGATTTCGTCGCTTACGACGCCTATCTGAttt                              >  1:3819549/1‑101 (MQ=255)
                                                                       aaaCTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAttt                              >  2:4045724/1‑101 (MQ=255)
                                                                       aaaCTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAttt                              >  2:2043152/1‑101 (MQ=255)
                                                                         aCTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCg                            >  1:1278164/1‑101 (MQ=255)
                                                                          cTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCgg                           >  2:173185/1‑101 (MQ=255)
                                                                            gCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGt                         <  1:1723591/101‑1 (MQ=255)
                                                                            gCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGt                         >  2:1116845/1‑101 (MQ=255)
                                                                            gCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATATGATTTCGGGt                         >  2:1933661/1‑101 (MQ=255)
                                                                             cTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTc                        <  1:2259561/101‑1 (MQ=255)
                                                                              tttC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCg                       >  1:2616436/1‑101 (MQ=255)
                                                                               ttC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGc                      >  1:788801/1‑101 (MQ=255)
                                                                               ttC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGc                      >  2:2351128/1‑101 (MQ=255)
                                                                                   g.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTcc                   >  2:3441424/1‑101 (MQ=255)
                                                                                      gggACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTcccc                 >  2:2966089/1‑101 (MQ=255)
                                                                                           tGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCcgcg            >  2:1627515/1‑101 (MQ=255)
                                                                                            ggCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTTGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCgg           <  1:1089474/101‑1 (MQ=255)
                                                                                               cAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCggggg        >  2:3185802/1‑101 (MQ=255)
                                                                                               cAGCCGCCCGCGCCAAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCggggg        <  2:3413865/101‑1 (MQ=255)
                                                                                                aGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGc       <  2:788719/101‑1 (MQ=255)
                                                                                                   cgcccgcGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCGTg    <  2:2512771/101‑1 (MQ=255)
                                                                                                    gcccgcGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCGTgg   <  1:1652568/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
CGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.G.TGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCGTGGA  >  NC_014034/1396464‑1396663

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG