Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,066 A→G 2.0% 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.0% 10.5 ‑2.0 0.0633D197D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (10/6):  ref base (265/348):  total (277/356)
Fisher's exact test strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 862 total aligned reads displayed.

GGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGG  >  NC_014034/221966‑222166
                                                                                                    |                                                                                                    
 gtccctgggctcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGAt                                                                                                     <  2:1036155/101‑1 (MQ=11)
 atccgtgggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAt                                                                                                     >  1:2716335/2‑101 (MQ=11)
  tccgtgggatctgtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                                    <  2:3732944/101‑1 (MQ=14)
  tccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGCAtc                                                                                                    <  2:3955194/101‑1 (MQ=11)
  tccgggggatcgggcgggtccgggggatcggtcggCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                                    <  1:4209479/101‑1 (MQ=14)
  gccgggggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                                    <  1:615824/100‑1 (MQ=11)
  gcccggggatcgggcgggtccgggggatcggtcggATCCGTGGGATCGGGCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGAtc                                                                                                    <  1:105116/96‑1 (MQ=11)
   tcgtgggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAtcc                                                                                                   <  2:1331365/100‑2 (MQ=11)
   ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTACGGGGGGGCGGTCGGATGCGGGGGATCGGTCGGATCAGTGGGAtag                                                                                                   >  2:3515860/1‑99 (MQ=11)
   ccgtgggatcggtcggatccgtgggatcggtcggGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGTTCGTTGGGAGCGCTCAGATCCGTGGGAtcg                                                                                                   >  2:3845194/1‑101 (MQ=11)
    cgtgggatgggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGAATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGAtcgg                                                                                                  >  1:3407624/1‑101 (MQ=11)
     gtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTGGGAtcggg                                                                                                 >  2:3896356/1‑100 (MQ=11)
     gtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGCTCGGTCGGGTCCGTGGGATCGGTCGGCCCCGTGGGATCGGTCGGGTCCCTGGGAccggc                                                                                                 >  1:4403440/1‑100 (MQ=11)
     gtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                                                 >  1:2817207/1‑101 (MQ=14)
     gtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTCGGAtgggt                                                                                                 >  1:1493980/1‑101 (MQ=11)
      ttgggtcggtcggagccgtgggatcggtcggCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGAtcggtc                                                                                                <  2:4574783/96‑1 (MQ=11)
      tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                                >  1:924531/1‑101 (MQ=14)
      tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                                >  1:2667794/1‑101 (MQ=14)
       gggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGGTGCGTGGGGTCGGTCGGATGCATGGGATCGGCCGGGTCCCCGGGAtcggtcg                                                                                               >  1:3219/1‑101 (MQ=14)
       gggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGATCCGGGGGATCGGTCGGATCGGTGGGAtcggtcg                                                                                               >  2:2149933/1‑101 (MQ=11)
       gggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGTTCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGAtcgggcg                                                                                               >  2:4150676/1‑101 (MQ=14)
        ggagctgtcggctctgtgggatcggtcggCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                              <  2:4091724/95‑1 (MQ=11)
            cgggcggctccggaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtcc                                                                                          <  2:1848073/101‑1 (MQ=11)
             ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGCTCGgtcgggtcgg                                                                                         >  2:908137/1‑99 (MQ=11)
             ggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccg                                                                                         <  2:4487683/101‑1 (MQ=255)
              gtggtgtccgtggggtccgtcggATCGGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgt                                                                                        <  2:4004638/96‑1 (MQ=14)
              gtcgggtccgtggggtccgtcggATCCGTGGGATCGGCCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                                                                                        <  1:3936528/101‑1 (MQ=14)
               tcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggggccgtg                                                                                       >  1:3873746/1‑100 (MQ=11)
               tcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtc                                                                                       <  1:125622/101‑1 (MQ=11)
                cgggtccgtgggatcgggcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcg                                                                                      <  1:3283602/101‑1 (MQ=11)
                cggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGCATCGGTCGGGGCCGTGGGATCGgtgggatccgtgg                                                                                      >  2:799619/1‑99 (MQ=11)
                cggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcg                                                                                      <  1:3065660/101‑1 (MQ=11)
                 gggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgacgggtgcggggg                                                                                     >  2:2601970/1‑97 (MQ=11)
                 gggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgtggg                                                                                     <  1:173160/101‑1 (MQ=11)
                 ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                                     >  1:1830145/4‑101 (MQ=255)
                 ggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgg                                                                                     >  1:86043/4‑101 (MQ=11)
                  gatccgtaggatgggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                                                    <  1:1003866/99‑1 (MQ=11)
                   gtgcgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgtgggat                                                                                   <  1:799693/98‑3 (MQ=11)
                   agcggtcgggtcggtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgggt                                                                                   <  2:4657572/91‑1 (MQ=11)
                    tccgtggggtccgtcggATCCGTGGGATCGGTCGGTTCCGTGGGATCGGTCGGATCCGTGGGACCGGTCGGGTCTGTGGGGTCGgtggggtcgggcggggc                                                                                  >  2:3133564/1‑99 (MQ=11)
                       gtgggatcggtcggATCCGTGGTATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGCCGGGACCGTAGGATCGgtcgggtccgtgggatcggc                                                                               >  1:3133218/1‑98 (MQ=11)
                        gcgggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTg                                                                              <  2:3198379/93‑1 (MQ=11)
                            gtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATAGGTCGGGTCCGTCGGGCTCGTGGGAt                                                                          >  2:722339/8‑101 (MQ=11)
                            gtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGag                                                                          >  2:269094/5‑100 (MQ=11)
                            ggcgggcgcATGCGTGGGGGCGGGCGGCTCTGTGGGATCGGTCGCCTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAt                                                                          <  2:3058182/99‑1 (MQ=2)
                             tccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGGtc                                                                         >  2:580725/7‑101 (MQ=11)
                              cggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtcc                                                                        <  2:1702268/101‑2 (MQ=11)
                              ccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg                                                                        <  1:706315/96‑1 (MQ=14)
                               ggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccg                                                                       >  1:2860732/1‑99 (MQ=11)
                               cgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcgg                                                                       >  2:115367/5‑101 (MQ=14)
                               cgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcgg                                                                       >  1:15768/5‑101 (MQ=14)
                                 tcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc                                                                     <  1:2497332/101‑1 (MQ=11)
                                 tcggATCGGTCGGATCCTTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                     <  2:612909/101‑1 (MQ=17)
                                 tcggATCCGTGGGATCGGTCGGATCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGAGCGGTCGGGGCCCGGGGGGCCGTCGGAtccgtg                                                                     >  2:2924872/1‑100 (MQ=11)
                                  cggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTGCGGGGGGTCGGTGGGAtccgtgg                                                                    >  2:521604/1‑99 (MQ=11)
                                    ggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcgggcggg                                                                  >  1:506346/3‑100 (MQ=14)
                                      tCCGGGGGATCGGTCGGGTCCGTGGGATCGGACGGATCCGTGGGATCGTTCGGGTCCGCCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTc                                                                <  1:1185119/101‑1 (MQ=14)
                                        ggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCg                                                              >  2:1417577/5‑101 (MQ=17)
                                        cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcggggccg                                                              >  1:2385721/1‑96 (MQ=14)
                                         gtcggATCCGTGGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGt                                                             <  1:4501247/98‑1 (MQ=17)
                                         gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGt                                                             <  2:2816991/101‑1 (MQ=14)
                                          tcggATCCGTGGGGTCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                            <  2:3645592/99‑1 (MQ=18)
                                          tcggATCCGTGGGGTCCGTGTGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                            <  2:2755818/99‑1 (MQ=18)
                                           gggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGGgg                                                           >  1:3499769/1‑101 (MQ=14)
                                           cggATCTGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                           <  2:1225581/100‑1 (MQ=255)
                                           cggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGCCCGTGGGATCGGTCGGATCCGTgg                                                           >  2:2463002/2‑101 (MQ=18)
                                           aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                                           >  1:256266/2‑99 (MQ=14)
                                            ggAGCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTcgg                                                          <  2:1655779/101‑1 (MQ=14)
                                              aTCGGTCGGATCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                                        <  2:2330595/101‑1 (MQ=11)
                                                   tggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                   >  1:216157/3‑101 (MQ=255)
                                                    ggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                  >  1:2266865/2‑101 (MQ=255)
                                                        tCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                              <  2:2385559/101‑1 (MQ=17)
                                                         ccGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGcc                                             >  2:4411427/1‑101 (MQ=17)
                                                          cGTGGGATCTGGCGGACCCGTGGGATCGGTGGGGTCCGTGGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCg                                            <  2:3319341/101‑1 (MQ=2)
                                                           gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                                           >  1:357352/1‑101 (MQ=255)
                                                            tGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                          >  1:2934503/1‑101 (MQ=255)
                                                             gggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCggag                                         >  2:2882408/1‑98 (MQ=17)
                                                             cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGGTCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGgg                                         >  2:1968857/2‑101 (MQ=17)
                                                              ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGACCGGTGGGGTCCGTGGGACCGGTCGGGTCGCTGggggcgcggcg                                        >  2:3356394/1‑93 (MQ=11)
                                                               gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggga                                       >  1:2040465/1‑100 (MQ=17)
                                                                gtcggtccgATCCGTGGGCTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggat                                      <  2:4295817/100‑3 (MQ=14)
                                                                aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggggt                                      >  1:4309263/1‑100 (MQ=18)
                                                                aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                                      <  2:3446796/101‑3 (MQ=14)
                                                                  cggtcggATCCTTGCGATCGGTGGGGTCCGTGGCATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcg                                    <  1:4023107/101‑5 (MQ=14)
                                                                    gtcggATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGGCCGGTCGGGTCCGTaggatcggt                                  >  1:1673454/1‑95 (MQ=14)
                                                                     tcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtg                                 >  2:300212/1‑101 (MQ=255)
                                                                     tcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtc                                 <  2:2522783/101‑8 (MQ=11)
                                                                     gcggATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtc                                 <  2:2034848/100‑8 (MQ=14)
                                                                      cggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcg                                >  1:4102190/1‑93 (MQ=11)
                                                                       ggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGAgggatcggtcgg                               >  1:138866/1‑92 (MQ=11)
                                                                        gggccGTGGGGTCGGGCGGGGCCGTGGGATCGGTCGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGa                              <  1:4595249/98‑1 (MQ=18)
                                                                        gATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGGCGGGTCCGTgggatcgggcggg                              >  1:3674423/1‑91 (MQ=11)
                                                                           ccGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcg                           >  1:3775256/1‑101 (MQ=255)
                                                                            cgggggAGCGGGGGGGTCCGTAGGATCGGGCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcgg                          <  1:1979261/98‑1 (MQ=21)
                                                                               gggGTCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcg                       <  1:3010780/101‑1 (MQ=255)
                                                                                            tCCGTAGGATCGGTCGGGTCCGTCGGGTCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgcc          <  1:4162808/101‑1 (MQ=255)
                                                                                                  ggATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCgccgccgcc    <  1:580786/101‑1 (MQ=255)
                                                                                                  ggATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCgccgccgcc    >  2:2237881/1‑101 (MQ=255)
                                                                                                   gATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGccgccgccg   >  2:4233628/1‑101 (MQ=255)
                                                                                                    aTCGGTCGGCTCCGTCGGGGCCGTGGGATCGGTCGGGTCCGCGGGATCGGTCTGGGCCGTGGGGCTGGTGGGATCGGTCGGACCGGTGGCGCCGCCGCCgg  >  2:1550786/1‑101 (MQ=21)
                                                                                                    |                                                                                                    
GGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGG  >  NC_014034/221966‑222166

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