Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,066 A→G 2.6% 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.6% 23.9 ‑0.2 0.0795D197D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/13):  ref base (343/424):  total (355/440)
Fisher's exact test strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 1047 total aligned reads displayed.

GGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGG  >  NC_014034/221966‑222166
                                                                                                    |                                                                                                    
 tcttccgatctcggtcggagccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                     <  1:4113566/91‑1 (MQ=14)
 tctcgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                     <  2:4665892/98‑1 (MQ=11)
 ggccgtggggtcgggcgggtccgtggggtcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                     <  1:2756345/99‑1 (MQ=14)
 gaccgtgtgatccgtcgggtccgtgggaccggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                     <  2:5394203/99‑1 (MQ=11)
  tccttgggttcggtctgatccttgggttcggtcgtATCCGTGGTATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGAtc                                                                                                    <  1:1513869/101‑1 (MQ=11)
  tccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                                    >  1:4433911/1‑101 (MQ=14)
   ccgtggggtcggtcgggtccgtgggatcggtcggGGCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                                                                   <  1:2181278/101‑1 (MQ=11)
   ccgtgggatcggtcggatccgtgggagcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGAtcg                                                                                                   <  2:899862/101‑1 (MQ=11)
    cgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCCGTCGTGTCCGTGGGATCGGTCGGGTCCCTGGGATCGGTCGCGTCCTTGGGAtcgg                                                                                                  >  1:5755369/1‑101 (MQ=11)
    cgtgggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtccg                                                                                                  >  2:1332928/1‑99 (MQ=11)
     gtggggtccgtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                                 <  2:1011482/101‑1 (MQ=14)
     gtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCGGCGGGAtcggt                                                                                                 >  2:838623/1‑101 (MQ=14)
      tgggatcggtcggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGAtcgggc                                                                                                >  1:3952069/1‑99 (MQ=11)
       gggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGTTCGGGTCCGGGGGATCGGTCGGGACCGGTGGAtttttag                                                                                               >  2:2957013/1‑95 (MQ=1)
        ggagtggacgggtccgtggactggtccggGTCCGAGGGATCGGTCGGGTCCGTGGGATAGGTCGGATCCGTGCGATAGGTCGGGTCCGTGGGGtcggtggg                                                                                              >  2:180883/9‑101 (MQ=2)
          tccgtccgggtccgtgtgatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGTGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt                                                                                            <  2:5532299/95‑1 (MQ=11)
            cggtcgggtgcgtggggtcggtgggGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                                                                          <  1:5472893/101‑1 (MQ=11)
            cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGGCCGTGGGATCGgtcgggtcc                                                                                          >  1:1125334/1‑101 (MQ=11)
            cgggcgggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatcc                                                                                          <  1:1362309/101‑1 (MQ=11)
             ggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTTGGGGCCGTGGGATCGgtcgggaccg                                                                                         >  2:5705180/1‑101 (MQ=11)
             ggtcgggtccgcggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccg                                                                                         >  2:3212298/1‑101 (MQ=11)
             ggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccg                                                                                         >  1:310377/1‑101 (MQ=11)
             ggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccg                                                                                         >  1:502489/1‑101 (MQ=11)
              gtcgggtccgtcggatcggccggGTCCGTGGGATCGGTCGGGTCCGTGCGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAACGatcgggtccgt                                                                                        >  1:2862388/1‑101 (MQ=11)
               tcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtc                                                                                       <  2:3114884/101‑1 (MQ=11)
               tcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGAACGgacgggtccggg                                                                                       >  2:3603995/1‑99 (MQ=14)
               tcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtagatcggaag                                                                                       >  1:4955135/1‑91 (MQ=14)
                cggatccgttggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcg                                                                                      <  1:4301436/101‑1 (MQ=11)
                  ggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGGCCGTGGGATCGgtcgggtccgtggga                                                                                    >  2:4548153/1‑100 (MQ=11)
                    tccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggctccgtgggatc                                                                                  >  1:1627629/1‑101 (MQ=11)
                     ccgtgggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                                 <  1:100380/101‑1 (MQ=14)
                     ccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGGGGGATCGgtcggctccggggggtgc                                                                                 >  2:1094962/1‑99 (MQ=11)
                     ccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTCGGATCCgcggggtccgtcgggtcc                                                                                 >  2:5237618/1‑101 (MQ=14)
                      cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCTCGgtcgggtccgtgggagatc                                                                                >  2:1865217/1‑96 (MQ=11)
                      cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGtcgggtccgtggggtccg                                                                                >  1:5161733/1‑101 (MQ=14)
                      cgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggatcgg                                                                                <  1:1996529/101‑3 (MQ=255)
                      cgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCGGtggggtccggcgggtccg                                                                                >  2:118550/1‑101 (MQ=11)
                       gtggggtacgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgt                                                                               <  2:5161447/92‑1 (MQ=11)
                       gtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggt                                                                               >  1:1229353/1‑101 (MQ=11)
                       gtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggt                                                                               >  2:707847/1‑101 (MQ=11)
                       gtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggt                                                                               >  1:4631241/1‑101 (MQ=11)
                       gtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggggtccgtcgggtccgt                                                                               >  1:4224916/4‑101 (MQ=14)
                        tggggtccgtcggATCCGGGGGATCGGCCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGCGGGGTCGgccgggaccgtggggacggtc                                                                              >  2:2474915/1‑95 (MQ=2)
                         gggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCGGTgg                                                                             >  2:3323535/1‑101 (MQ=11)
                         gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgg                                                                             <  1:3227381/101‑1 (MQ=11)
                         cggttcggttggGTTCGTGGGATCGGTCTGGTCCGTCGGCTCGGGGGGATCCGTGGGATCTGTCGGGTCAGTGGGATCGgtcgggtccgtcggatcgggcg                                                                             <  1:3304972/97‑6 (MQ=2)
                          gggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTggg                                                                            >  2:271014/10‑101 (MQ=11)
                          gggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTcgg                                                                            <  1:35248/95‑1 (MQ=11)
                          ggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTggg                                                                            >  1:4537971/1‑101 (MQ=11)
                           ggtcggtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGa                                                                           <  2:2561160/99‑1 (MQ=14)
                           gatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGa                                                                           <  1:2119136/101‑1 (MQ=14)
                            gtccgggggATCCGTGGGATCGGTCGGTTCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAt                                                                          <  1:4855131/94‑1 (MQ=11)
                            atcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAt                                                                          >  1:2924899/1‑101 (MQ=11)
                             tccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtc                                                                         <  1:1843725/95‑1 (MQ=14)
                             tccgggggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtc                                                                         <  2:2345082/95‑1 (MQ=14)
                              cggtcggGTCCGTGGGATCGTGTGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcg                                                                        <  1:2506777/101‑1 (MQ=11)
                              cggccggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtcc                                                                        <  2:4290622/101‑2 (MQ=11)
                               cgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGATCGAATCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGGACCGTCGGAtccg                                                                       <  2:5782228/100‑3 (MQ=11)
                                gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTGGGATCCGTGGGGGCCGTGGGAtcggt                                                                      >  2:2076077/1‑101 (MQ=1)
                                gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggt                                                                      <  1:1712459/101‑1 (MQ=14)
                                 cgatctCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                                                                     <  2:1667551/96‑1 (MQ=11)
                                  cggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcg                                                                    <  1:3837058/101‑1 (MQ=14)
                                   gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGTTCGGTTGGCTCGGTGGGGCCCCTGGGGtcggccgg                                                                   >  1:3415605/4‑101 (MQ=11)
                                   gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgg                                                                   <  2:4355420/98‑1 (MQ=11)
                                   gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgg                                                                   >  2:4063857/4‑101 (MQ=11)
                                   ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg                                                                   <  1:1302496/101‑1 (MQ=11)
                                     gtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgggt                                                                 <  2:1584611/100‑3 (MQ=11)
                                     aTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCGGTGGGATCCGTGGGAt                                                                 >  1:5648149/1‑101 (MQ=11)
                                        cGTGGGATCGGTCGGGTCCGCGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGGTCCGCGGGGtcggtcggggccg                                                              >  1:3054419/1‑96 (MQ=14)
                                         gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCgg                                                             >  1:5400310/1‑100 (MQ=14)
                                          taggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGtc                                                            <  1:476235/99‑2 (MQ=14)
                                           cggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                           >  2:285289/2‑101 (MQ=255)
                                           cggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                           >  2:3906945/2‑101 (MQ=255)
                                            ggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTggg                                                          <  2:5038015/101‑1 (MQ=14)
                                             gATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                                         >  1:4027943/1‑101 (MQ=17)
                                              aTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAGCGGTCGGATCCGTGGGAt                                                        >  2:212497/1‑101 (MQ=18)
                                                cggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                      <  2:2230143/101‑1 (MQ=14)
                                                   tcgttttcGTTGGATCTATCGGATCGGTGCGAGCGGTCGGGTCCTTGGGATCGCTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                   <  2:3493044/94‑1 (MQ=9)
                                                    cggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg                                                  >  1:2330857/1‑101 (MQ=17)
                                                    cggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg                                                  >  1:718301/1‑101 (MQ=17)
                                                     gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGGCCGTGGGAtcggtcgg                                                 >  2:1275622/1‑101 (MQ=17)
                                                     ggatCCGTAGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg                                                 >  2:3065465/4‑101 (MQ=14)
                                                         ccGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                             <  1:907104/101‑1 (MQ=17)
                                                           gTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGt                                           >  1:1949722/1‑101 (MQ=17)
                                                            tGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGta                                          <  1:963172/101‑2 (MQ=17)
                                                             cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgg                                         >  2:193980/2‑101 (MQ=17)
                                                              gggtctgtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTagg                                        <  2:533811/95‑1 (MQ=14)
                                                              ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggg                                        <  1:3685992/101‑1 (MQ=17)
                                                               gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggga                                       <  2:3992019/101‑2 (MQ=17)
                                                                aTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCGGTGGGGCCCGCGGCGTCGGTGGCGGGCGTggggt                                      >  2:4200935/1‑100 (MQ=11)
                                                                 tcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatc                                     <  1:5340555/101‑4 (MQ=14)
                                                                   ggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcgg                                   >  1:2681783/1‑96 (MQ=14)
                                                                     gcggATCCGTGGGTTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGTTCCGTGGGATCGGTCGGGTCCGTaggatcggtc                                 <  2:4469939/100‑8 (MQ=14)
                                                                          tCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc                            <  2:1480128/101‑1 (MQ=255)
                                                                              tgtgATCGGTCGTGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtc                        <  2:4941753/98‑1 (MQ=255)
                                                                                ggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcgg                      <  2:770555/101‑1 (MQ=255)
                                                                                     ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCGGTAGGACCGGTGGGGcccgtcggggcgg                 >  2:4813857/1‑96 (MQ=2)
                                                                                      gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCt                >  1:622402/1‑101 (MQ=255)
                                                                                      gtcggGTCCGTAGGATCGGGCGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGGCCGCt                <  2:5607446/101‑1 (MQ=255)
                                                                                        cggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGt              <  2:380805/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGG  >  NC_014034/221966‑222166

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG