Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,066 A→G 2.4% D197D (GAT→GAC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0660A→G 2.4% 19.2 ‑2.0 ‑0.1633D197D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/5):  ref base (298/319):  total (309/324)
Fisher's exact test strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.35e-01

Only 100 of 831 total aligned reads displayed.

GTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGG  >  NC_014034/221962‑222166
                                                                                                         |                                                                                                    
    ggtccgtgggatcggtcgggtccgt.ggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                                                                                      >  1:2451204/1‑101 (MQ=11)
    gctccgtgggatcggtcgggtccgt.cggatccgtgggATCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                                                                                      <  2:3408567/99‑1 (MQ=11)
     gtccgtgggatgggtcgggtccgt.ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                     <  2:2356824/101‑1 (MQ=11)
     gtccgtgggatcggtcgggtccgt.cggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTTGGATCCGTGGGATCGGTCGGGTTCGTGGGAt                                                                                                     >  1:4182538/1‑101 (MQ=14)
      accgtgggatcggtcgggtccgt.gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATTTTTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                                                    <  1:4631262/100‑1 (MQ=11)
       ccgtgggatcggtcggatccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtcc                                                                                                   <  1:2042265/101‑2 (MQ=11)
        cgtgggatcggtcgggtccgt.ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGAtcgg                                                                                                  >  2:2089125/1‑101 (MQ=11)
        cgtgggatcggtcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCGGTGGGAttgg                                                                                                  >  2:3707982/1‑101 (MQ=14)
        cgtgggatcggtcggatccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtccg                                                                                                  <  1:1188791/101‑3 (MQ=11)
        cgggggatcggtcgggtccgt.ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                                  <  1:4862407/98‑1 (MQ=11)
         gtggggtccgtcggatccgt.gggatcggtcggATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGCCGGGTCCGTGGGATCGGTCGGGTCCGAGGGAtcgga                                                                                                 >  2:1399501/1‑100 (MQ=11)
           ggggtccgtcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGAtcgtgcg                                                                                               >  1:1480797/1‑97 (MQ=11)
           gggatcggtcgggtccgt.ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGAtcggtcg                                                                                               >  2:2533318/1‑101 (MQ=11)
           gggatcggtcgggtccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAgcggtcg                                                                                               >  2:4915913/1‑101 (MQ=14)
            ggatcggtcgggtccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGGGGGAtcggtcgg                                                                                              >  1:4586031/1‑101 (MQ=11)
            ggatcggtcgggtccgt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGGCGGGTCCGTGGGGtcgggcgg                                                                                              >  2:3594447/1‑101 (MQ=11)
             gatctgccgggtccgt.gggtgccgtcggGTCCGTGGGATCGGTCGGGTCCGTGCGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggg                                                                                             <  2:1040713/101‑1 (MQ=11)
              gtcggtcgggtccgt.cggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                                                                            <  1:1534433/100‑1 (MQ=11)
              atcggtgggatccgt.gggatcggccggGTCCGGGGGATCGGTCGGGTCCGTGGGACCGGTCAGGGCCGTGGTATCGGTCGGGGCCGTGGGTtgggtcggcc                                                                                            >  2:850949/1‑99 (MQ=2)
               tctgtcggggccgt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                                                                                           <  2:900145/98‑1 (MQ=11)
               tcggtcgggtccgt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGggcggggc                                                                                           >  2:2128977/1‑99 (MQ=11)
               tcggtagggtccgt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                                                                                           <  1:1437389/101‑1 (MQ=11)
                gggtggggtccgtcggggtc.gtgggATCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                                                                          <  1:998666/100‑1 (MQ=11)
                cggtcgggtcggt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                                                                          <  1:3263898/101‑1 (MQ=11)
                cggtcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGGTCGggcggatcc                                                                                          >  1:4128132/1‑101 (MQ=11)
                 ggtcggatccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGGCAGGGGGATCGgtcgggtccg                                                                                         >  2:453301/1‑101 (MQ=11)
                 ggtcggatccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCCgtggggtccg                                                                                         >  2:3177552/1‑101 (MQ=11)
                  gtcgggtccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCGGGGGGATCGATCGGGTCCGTGGGGTCGgccgggtcggt                                                                                        >  2:2193405/1‑101 (MQ=11)
                   tcgggtccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCTTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgccggggccgtg                                                                                       >  1:4531524/1‑100 (MQ=11)
                   tcgggtccgt.cggatcggtgggATCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtg                                                                                       <  1:3797117/101‑2 (MQ=11)
                   tcgggtccgt.cggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGgtcgggtccgtc                                                                                       >  1:900250/1‑101 (MQ=11)
                   tcggatccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtc                                                                                       <  2:2690088/101‑1 (MQ=11)
                   gcgggtccgt.ggggtccgtcggATCCGTGGGATCGGTTGGATCCGTGGGATCGGTCGTATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtg                                                                                       <  1:4312772/100‑2 (MQ=11)
                    cggatccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcg                                                                                      >  1:3735997/1‑101 (MQ=11)
                    cggatccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGCCGGATCCGTGGGAACGGGCGGGTGCGTAGGAACCgtggggcgcgtcg                                                                                      >  1:1089100/1‑101 (MQ=2)
                     gggtccgt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggg                                                                                     >  2:4978312/1‑101 (MQ=11)
                     ggatccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGGTCGGTCGGGGCCGTCGGATCCgtggggtccgtcgg                                                                                     >  2:2940424/4‑101 (MQ=11)
                      gatccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGAAGGATGGGTCGGGTCCGTAGGATCGgtcgggcccgtcgga                                                                                    >  1:4736861/3‑100 (MQ=17)
                       gtgcgt.gggatgggtcggATGCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgggat                                                                                   <  1:747287/98‑3 (MQ=11)
                        tccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtgggggccgtcggggc                                                                                  >  2:124864/1‑99 (MQ=11)
                         ccgt.ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtcc                                                                                 <  2:691258/101‑1 (MQ=11)
                         ccgt.ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGATCGgccgggtccgtggggtcg                                                                                 >  2:5049391/1‑100 (MQ=11)
                         ccgt.aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtgggggccgtcgggtcc                                                                                 >  1:335467/1‑101 (MQ=14)
                          cgt.ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtcgg                                                                                >  1:1753431/1‑99 (MQ=11)
                           gt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCTGGGCCGGGGGATGGgtcgggtcggtcggatcggc                                                                               >  2:659155/1‑98 (MQ=11)
                            t.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGAGATCGGTCGGGGCCGTGGGATCGGAcggggccgcggggtccgtc                                                                              >  1:653321/1‑100 (MQ=11)
                              cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTgg                                                                             <  1:2125042/91‑1 (MQ=11)
                              cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTgg                                                                             >  2:170793/11‑101 (MQ=11)
                              cgggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGGgg                                                                             >  1:2260910/11‑101 (MQ=11)
                               gggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTggg                                                                            >  1:4931036/10‑101 (MQ=11)
                               gggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGTTCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTggg                                                                            >  2:2765094/7‑101 (MQ=11)
                               gggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGGCCGTcgg                                                                            >  2:2028635/7‑101 (MQ=11)
                               ggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtggttcgg                                                                            >  2:18801/1‑94 (MQ=11)
                                ggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGa                                                                           <  1:165614/96‑1 (MQ=11)
                                 gtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGac                                                                          >  2:4042010/8‑100 (MQ=11)
                                 gcgcggaagATCCGTGGGGTAGGTGGGGTGCGTGGGAGGGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAt                                                                          <  1:2608537/93‑1 (MQ=2)
                                  gccgtgggATCGGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtc                                                                         <  1:1934419/95‑1 (MQ=14)
                                   cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcg                                                                        >  1:840959/1‑101 (MQ=11)
                                   cggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCAGTCGGATCGGTCGGAtcg                                                                        >  2:3271322/1‑101 (MQ=11)
                                   ccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg                                                                        <  1:1020335/96‑1 (MQ=14)
                                    cgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccg                                                                       >  2:610933/2‑99 (MQ=11)
                                     gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGGTCCGTGGGAtcggg                                                                      >  1:1271579/1‑100 (MQ=11)
                                      tcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                                                                     <  1:2700319/101‑1 (MQ=11)
                                      tcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                                                                     >  1:2991097/1‑101 (MQ=11)
                                       cggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCTTGGGATCGGTCGGATCCGTGGGGTCCGAGGGAtcggtcg                                                                    >  2:4618690/1‑101 (MQ=11)
                                       cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtgg                                                                    >  1:3600135/1‑99 (MQ=11)
                                        ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtggg                                                                   >  2:4820029/1‑101 (MQ=11)
                                             cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGGtcggtcggggccg                                                              >  2:4695108/1‑96 (MQ=11)
                                              gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGt                                                             >  1:2490084/1‑101 (MQ=14)
                                               tGGGACCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTGCGTGGGGCCCGTCGGAGGGGTGGGGTCGGCGGGATcgttg                                                            >  2:4157416/1‑97 (MQ=2)
                                                  gATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCggcggg                                                         >  2:1554230/1‑96 (MQ=17)
                                                   aTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                                        <  2:3505279/101‑1 (MQ=255)
                                                     ccggggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                      <  1:3967990/96‑1 (MQ=255)
                                                       gtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                    <  1:2798854/101‑1 (MQ=17)
                                                          gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGCtcggtcgg                                                 >  2:1143537/1‑101 (MQ=17)
                                                           ggTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg                                                <  1:5046256/101‑1 (MQ=255)
                                                             tCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTc                                              >  1:3126489/1‑101 (MQ=11)
                                                              ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTcg                                             >  2:1649896/1‑100 (MQ=11)
                                                              ccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                             <  2:4233825/101‑1 (MQ=14)
                                                              ccGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTcc                                             >  1:3879751/1‑101 (MQ=17)
                                                               cGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                            <  2:3887955/101‑1 (MQ=17)
                                                                gTGGGATCTGGCGGATCCGTGGGGTCGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                           <  2:2448758/101‑1 (MQ=17)
                                                                  cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgg                                         >  2:4435467/2‑101 (MQ=17)
                                                                   ggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGggg                                        >  1:3505337/1‑101 (MQ=255)
                                                                    gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTagga                                       >  1:47946/1‑100 (MQ=17)
                                                                     gtcggtcggATCCGGGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                                      <  1:1334935/100‑3 (MQ=14)
                                                                       cggtcggGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTcgggtcc                                    >  2:783493/1‑98 (MQ=14)
                                                                          tgggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtc                                 <  1:604803/99‑8 (MQ=11)
                                                                           cggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcg                                >  2:3987673/1‑93 (MQ=14)
                                                                           cggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcg                                >  1:2381946/1‑93 (MQ=14)
                                                                            gggtCCGTGGGATCGGTCGGGTCCGTAGGGTCGGTTTGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGCGACCGGTCGGGCCCACGGGGCTGGGggg                               >  2:2839076/4‑101 (MQ=18)
                                                                            ggATCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcgg                               <  2:26102/101‑10 (MQ=14)
                                                                            ggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcgg                               >  1:4413193/1‑92 (MQ=14)
                                                                              aTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTggggc                             >  2:132272/1‑99 (MQ=14)
                                                                                   tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCGGTCGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGGCCGGTCGGTtccggc                        >  2:1467129/1‑97 (MQ=14)
                                                                                    gggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcg                       >  2:3796161/1‑101 (MQ=255)
                                                                                    gggATCGGTCGGGTCCGTAGGATCGGCCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcg                       <  1:4241168/101‑1 (MQ=255)
                                                                                        tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCCTGGGATCGGTCGGATCCGTGGGATCGGGCGGGGCCGTGGGGCGCGTTGGATCGGACGCAcc                   >  2:3016220/1‑101 (MQ=18)
                                                                                           gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAACGCt                >  2:4358931/1‑101 (MQ=255)
                                                                                              gggGCCGTAGGACCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTc             <  1:4057574/101‑1 (MQ=255)
                                                                                                         |                                                                                                    
GTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGG  >  NC_014034/221962‑222166

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 40 ≤  ATCG