Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 448,780 A→C 2.6% intergenic (+73/+147) RCAP_rcc00419 → / ← pykA1 NADH:flavin oxidoreductase/NADH oxidase/pyruvate kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034448,7800A→C 2.6% 7.1 ‑0.3 ‑0.5239intergenic (+73/+147)RCAP_rcc00419/pykA1NADH:flavin oxidoreductase/NADH oxidase/pyruvate kinase
Reads supporting (aligned to +/- strand):  new base (6/3):  ref base (166/64):  total (172/67)
Fisher's exact test strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.42e-01

Only 100 of 337 total aligned reads displayed.

TGCGCTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGTCGTGCCGCAGCG  >  NC_014034/448680‑448879
                                                                                                       |                                                                                                   
tGCGCTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.Ca                                                                                                     >  2:6331248/1‑101 (MQ=255)
 gcgcTGCTGCTGATGAAGGACGTTTGACGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CAc                                                                                                    <  2:4200173/101‑1 (MQ=255)
  cgcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACt                                                                                                   >  2:1807816/1‑101 (MQ=255)
  cgcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACt                                                                                                   >  2:386224/1‑101 (MQ=255)
  cgcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACt                                                                                                   <  1:5103053/101‑1 (MQ=255)
   gcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTc                                                                                                  <  1:1430728/101‑1 (MQ=255)
   gcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTc                                                                                                  <  2:193297/101‑1 (MQ=255)
   gcTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTc                                                                                                  <  1:4894526/101‑1 (MQ=255)
      gCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGc                                                                                               >  2:2192965/1‑101 (MQ=255)
      gCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGc                                                                                               >  1:2796469/1‑101 (MQ=255)
        ggCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCt                                                                                             <  2:2502544/101‑1 (MQ=255)
        ggCTGATGAAGCACGTTTGGGGCTGACTTGGCCTCGGGCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCt                                                                                             <  2:7221136/101‑1 (MQ=255)
         gCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTg                                                                                            <  2:4528502/101‑1 (MQ=255)
         gCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTg                                                                                            <  1:3539725/101‑1 (MQ=255)
         gCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTg                                                                                            >  1:5613063/1‑101 (MQ=255)
          cTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTgc                                                                                           <  1:5484526/101‑1 (MQ=255)
           tgatgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTgcg                                                                                          >  2:3193799/1‑101 (MQ=255)
            gatgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTgcgc                                                                                         >  1:6828937/1‑101 (MQ=255)
            gatgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACAGCCGAGGGCCCAAA.GG.CACTCGGACTgtgc                                                                                         >  1:4123804/1‑101 (MQ=255)
             atgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGcgcc                                                                                        >  1:320951/1‑101 (MQ=255)
              tgaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGcgccc                                                                                       <  1:8533797/101‑1 (MQ=255)
               gaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGcgcccg                                                                                      >  2:4726896/1‑101 (MQ=255)
               gaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGcgcccg                                                                                      >  1:8769090/1‑101 (MQ=255)
               gaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGcgcccg                                                                                      >  1:8161203/1‑101 (MQ=255)
               gaAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGcgcccg                                                                                      <  1:5353127/101‑1 (MQ=255)
                aaGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGcgcccgc                                                                                     <  2:7984769/101‑1 (MQ=255)
                 aGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGcgcccgcc                                                                                    >  1:6177612/1‑101 (MQ=255)
                  ggACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCgcccgccc                                                                                   <  1:169532/101‑1 (MQ=255)
                   gACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGcccgcccc                                                                                  <  1:7797126/101‑1 (MQ=255)
                   gACGTTTGAGGCTGACTTGGCCCCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGcccgcccc                                                                                  <  1:3942766/101‑1 (MQ=255)
                    aCGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCCcccgccccc                                                                                 >  1:4036632/1‑100 (MQ=255)
                      gTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCccgccccgcc                                                                               >  2:3397/1‑101 (MQ=255)
                      gTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCccgccccgcc                                                                               >  1:3261969/1‑101 (MQ=255)
                      gTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCccgccccgcc                                                                               >  2:3000802/1‑101 (MQ=255)
                       tttGAGGCTGACTTGGCCTCGGCCATTGCCGAGACC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCa                                                                              >  1:4771494/1‑101 (MQ=255)
                        ttGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCAt                                                                             >  2:8198017/1‑101 (MQ=255)
                        ttGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCAt                                                                             >  1:7429062/1‑101 (MQ=255)
                        tggAGCCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCAt                                                                             <  2:6413/99‑1 (MQ=255)
                         tGAGGCTGCCTTGGCCGCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATc                                                                            <  1:6648465/101‑1 (MQ=255)
                         tGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGGCTGGGCCCGCCCCGCCacg                                                                            >  2:7346306/1‑99 (MQ=255)
                         tGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCCTc                                                                            >  2:6722800/1‑101 (MQ=255)
                         tGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATc                                                                            >  1:1174549/1‑101 (MQ=255)
                         tGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATc                                                                            >  1:7508143/1‑101 (MQ=255)
                          gAGGCTGACTGGGCCTCGGCCATTGCCGAGGCC.AATCGCAACGGAGGTCCTGTGCCCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATcg                                                                           <  2:1496331/101‑1 (MQ=255)
                            ggCTGCCTTGGCCTGGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CCCTCGGCCTGCGCCCGCCCCGCCATcggc                                                                         <  1:6588807/101‑1 (MQ=255)
                            ggCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATcggc                                                                         >  2:4593817/1‑101 (MQ=255)
                              cTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATcggcgg                                                                       >  2:2822104/1‑101 (MQ=255)
                                gACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCggcgggg                                                                     >  1:8657985/1‑100 (MQ=255)
                                 aCTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCggcggggg                                                                    >  2:3269307/1‑99 (MQ=255)
                                  cTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGGGGTCCGGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGgcgc                                                                   <  1:925937/101‑1 (MQ=255)
                                  cTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCCGCGGgagc                                                                   >  2:240165/1‑101 (MQ=255)
                                   ttGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCggcggggcgc                                                                  >  2:6949666/1‑97 (MQ=255)
                                   ttGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGgggcg                                                                  >  2:7767899/1‑100 (MQ=255)
                                     ggCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCggcggggggggg                                                                >  1:3430373/1‑95 (MQ=255)
                                     ggCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCggcggggccagg                                                                >  2:6619196/1‑95 (MQ=255)
                                     ggCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGgggcgcg                                                                >  2:848053/1‑98 (MQ=255)
                                      gCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCggcgggggggggg                                                               >  1:3571238/1‑94 (MQ=255)
                                      gCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCggcggggggcggc                                                               >  1:4156119/1‑94 (MQ=255)
                                         tCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCcgcggggggcgggggg                                                            >  1:421866/1‑91 (MQ=21)
                                           ggCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGGAGGCGcgcc                                                          >  2:3466739/1‑101 (MQ=255)
                                             ccATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGCCGGCGcgccgc                                                        >  1:2347699/1‑101 (MQ=255)
                                              cATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGGGGCCCCgccgcc                                                       >  1:2478071/1‑101 (MQ=21)
                                              cATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGCGGGGGCgccgcc                                                       >  2:6130160/1‑101 (MQ=25)
                                                ttGCCGAGGCC.AATCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAG.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCgg                                                     <  2:7831143/101‑1 (MQ=255)
                                                  gcggAGGCC.AATCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCgggg                                                   <  1:6900498/98‑1 (MQ=255)
                                                   ccGAGCCC.AATCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGc                                                  <  1:1788178/101‑1 (MQ=255)
                                                    ccgggccaaaTCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGc                                                  <  2:9053643/93‑1 (MQ=255)
                                                    cGAGGCC.AACCGAAACAGAGGTCCTGTGCCCGCGGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCg                                                 <  1:8314140/101‑1 (MQ=18)
                                                     caGGCCAAATCCCAACAGAGTTCCTGTGCCCGCCGAGGCCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCg                                                 <  2:7441184/100‑1 (MQ=17)
                                                      gggCC.AGCCGCAACAGAGGTCCTGTGCCCGCCGGGGCCCCAAA.GG.CGCTCGGCTTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGa                                               <  2:5752030/100‑1 (MQ=14)
                                                      gggCC.AAGCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGa                                               <  1:9109447/100‑1 (MQ=255)
                                                      agccc.AACCGCAACAGAGGTCCTGTGCCCGCCGGGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGa                                               <  1:7857477/98‑1 (MQ=17)
                                                      aGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGGCGGCGCGCCGCCGCGGCggg                                               >  2:5978660/1‑100 (MQ=18)
                                                       ggCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGCGGGCGCCGCGCCGGGGGgggg                                              >  1:6013307/1‑99 (MQ=18)
                                                         cc.AATCGCAACAGAGGTCCTGTGCCCGCCGAGGGCCCAAA.GG.CCCTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAgcg                                            <  1:7545771/101‑1 (MQ=255)
                                                          c.AATCGCAACAGAGGTCCTGTGCCCGCCGGGGGCCCAAG.GG.CCCTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAgcgc                                           <  1:650384/101‑1 (MQ=17)
                                                          c.AACCGCAACAGAGGTCCGGTGCCCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAgcgc                                           <  2:2669991/101‑1 (MQ=255)
                                                              tCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGGGGCCGCGCCGCCGCGGGGGGGCGCGgg                                        >  2:4814415/1‑101 (MQ=18)
                                                               cGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGGCGGCCCGCGGCCCGGGGGGGGGGCAGGc                                       >  2:5353805/1‑101 (MQ=17)
                                                                 cAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGGCGGCGGGCCGGCCGGGCGGGGGGCAGGccc                                     >  2:7087957/1‑100 (MQ=17)
                                                                  acgggagGTCATGGGCCCCCCGAGGGCCTAGC.GC.CACTCGGTGTGGGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCctc                                    <  2:1113335/97‑1 (MQ=14)
                                                                     ggggttCCTGTGCCCGCCGGGGCCCCAAG.GG.CCCTCGGCCGGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCtgt                                 <  2:5579838/96‑1 (MQ=11)
                                                                     agagGCCTTGTGCCCGCCGGGGGCCCAAG.GG.CACTCGCCCGGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCtgt                                 <  1:2151172/101‑1 (MQ=14)
                                                                       ggttCCTGTGCCCG.CGGGGGCCCAAA.GGCCCCCCGGCCGGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCtgttg                               <  2:6289717/98‑1 (MQ=11)
                                                                        ggTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGGCGGCCCGCCGCCGGGGGGGGGCGCGGGCCGCcgtttt                              >  1:5812626/1‑99 (MQ=14)
                                                                         gcctgttgcccgccgAGGGCCCAAA.GG.CCCTCGCCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTtgtg                             <  2:8314043/92‑1 (MQ=11)
                                                                           ccggggcccgcggAGGGCCCAAG.GC.CCCTCGGCCTGCCCCCGCCCCGCCATCGGCGGGCGCAGTCGAGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCt                           <  1:1990043/94‑1 (MQ=11)
                                                                            cTGGGACCGCCGAGGGCCCAAG.GG.CCCTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGTGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCtt                          <  2:1655751/101‑1 (MQ=11)
                                                                                gcccgcggGGGGCCAAAG.GG.CACCCGCCCGGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTgggg                      <  1:104812/99‑1 (MQ=14)
                                                                                gACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGGCGGCGCGCCCCCCGGGCGGGGGGCGGGCCTCTGTTGTGCTTgtgg                      >  2:4056251/1‑101 (MQ=14)
                                                                                 ccggccgA.GGCCCAAA.GGCCCCCCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGc                     <  2:866460/98‑1 (MQ=11)
                                                                                  ccgccgAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGCAGGCGCGCCGCCCGGGGGGGGCGCCGGGCGCTGTTGTGCTTgtgggg                    >  2:5654617/1‑99 (MQ=14)
                                                                                  cccgcgggggCCCAAAGGG.CCCTCGGCCTGCCCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGc                     <  1:3353172/94‑1 (MQ=11)
                                                                                   cgccgAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGGCGGCGCGCCCCCCGGGCGGGGGGCCGGGCCCCGGTGTGCTTggggggg                   >  1:8027534/1‑98 (MQ=9)
                                                                                     ccgGGGGCCAAAG.GG.CACTCGGCCGGCCCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGc                 <  2:3571130/101‑1 (MQ=17)
                                                                                      cgggggCCCAAG.GG.CCCTCGCCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCg                <  1:1594644/98‑1 (MQ=14)
                                                                                           gCCCAAAGGG.CCCTGGGCCTGCCCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGTCg            <  1:4458905/101‑1 (MQ=14)
                                                                                            cccAAG.GG.CCCTCGGCCTGCCCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGTCGTg          <  1:4594129/101‑1 (MQ=255)
                                                                                             ccAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGGCGGCGCGCCGCCCGGGGGGGGGGCAGGGCCCCGTTTTGCTTTGGGGGGGCgggggggc         >  2:4890433/1‑95 (MQ=2)
                                                                                                      cACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGGGCGGGCGCGGCGCCCGGGCGGGGGGGAGGGCCCTGTTTTGCGTGGGGGGGGGGGGCGTGCCGCCgcg  >  2:1550502/1‑101 (MQ=9)
                                                                                                       |                                                                                                   
TGCGCTGCGGCTGATGAAGGACGTTTGAGGCTGACTTGGCCTCGGCCATTGCCGAGGCC.AATCGCAACAGAGGTCCTGTGACCGCCGAGGGCCCAAA.GG.CACTCGGCCTGCGCCCGCCCCGCCATCGGCGGGCGCAGTCGCGCCGCCGGGGCGGAGCGCAGGCCTCTGTTGTGCTTGGGGCAGGCGGTCGTGCCGCAGCG  >  NC_014034/448680‑448879

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG