Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,025 A→T 16.7% T303S (ACG→TCC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,0250A→T16.7% 114.7 ‑2.1 ‑0.5276T303S (ACG→TCG) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (17/22):  ref base (122/90):  total (161/115)
Fisher's exact test strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.97e-01

Only 100 of 290 total aligned reads displayed.

GGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGAT  >  NC_014034/1396925‑1397125
                                                                                                     |                                                                                                        
ggATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCCGCCGGa                                                                                                          >  2:1108903/1‑101 (MQ=255)
ggATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGa                                                                                                          <  2:3776703/101‑1 (MQ=255)
ggATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGa                                                                                                          >  2:1046852/1‑101 (MQ=255)
 gATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAc                                                                                                         <  1:3735938/101‑1 (MQ=255)
  aTCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.g                                                                                                       >  2:4014005/1‑101 (MQ=255)
    cGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GAt                                                                                                     <  2:2494166/101‑1 (MQ=255)
      aaaCAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCTCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGt                                                                                                   >  1:4331037/1‑101 (MQ=255)
         cAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...g                                                                                             <  2:1251902/101‑1 (MQ=255)
           gCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTc                                                                                           >  1:3512471/1‑101 (MQ=255)
             cTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGc                                                                                         >  2:461799/1‑101 (MQ=255)
              tGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCt                                                                                        >  1:809899/1‑101 (MQ=255)
                cgcgGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTcc                                                                                      >  2:1320461/1‑101 (MQ=255)
                  cgGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTcccc                                                                                    <  1:2018958/101‑1 (MQ=255)
                   ggCCGCATCGGGTTCCGCCTGCGTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTccccc                                                                                   <  2:3286297/101‑1 (MQ=255)
                   ggCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTccccc                                                                                   >  2:3539394/1‑101 (MQ=255)
                    gCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTcccccc                                                                                  <  1:2289942/101‑1 (MQ=255)
                    gCCCCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTcccccc                                                                                  <  2:4493936/101‑1 (MQ=255)
                       gCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCcgcg                                                                               >  2:2280191/1‑101 (MQ=255)
                       gCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCcgcg                                                                               <  2:756917/101‑1 (MQ=255)
                        cATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCgg                                                                              >  1:652429/1‑101 (MQ=255)
                         aTCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCggg                                                                             >  1:2390105/1‑101 (MQ=255)
                             gCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCg                                                                         <  2:855754/101‑1 (MQ=255)
                               ttCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTg                                                                       <  2:1676148/101‑1 (MQ=255)
                               ttCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTg                                                                       >  2:931131/1‑101 (MQ=255)
                                tCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGGgg                                                                      >  1:24040/1‑101 (MQ=255)
                                 ggatcccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTg                                                                       >  2:2988851/11‑101 (MQ=2)
                                 ccGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGa                                                                     >  1:2689549/1‑101 (MQ=255)
                                  cGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                    <  1:3148972/101‑1 (MQ=255)
                                   gCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                   >  2:1954706/1‑101 (MQ=255)
                                   gCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                   <  1:1119958/101‑1 (MQ=255)
                                   atcccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGa                                                                     >  1:1601526/9‑101 (MQ=2)
                                     cagcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                   >  2:2125911/8‑101 (MQ=2)
                                     cagcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                   >  2:2212014/8‑101 (MQ=2)
                                      tGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                <  1:562998/101‑1 (MQ=255)
                                      cgCGGGGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                <  2:3641658/100‑1 (MQ=255)
                                      ccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                  <  2:2653565/96‑1 (MQ=2)
                                      agcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                  >  1:1608934/7‑101 (MQ=2)
                                      agcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                  >  1:4130365/7‑101 (MQ=2)
                                        cGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                              <  1:2136283/101‑1 (MQ=255)
                                         ggcgTCGCCCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                               <  2:1396239/98‑1 (MQ=0)
                                         ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                               <  1:1034040/99‑1 (MQ=2)
                                         ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                               <  1:1806798/99‑1 (MQ=2)
                                         ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                               <  2:3071136/99‑1 (MQ=2)
                                         ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                               >  1:3243559/3‑101 (MQ=2)
                                          cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                              >  1:1380751/2‑101 (MQ=2)
                                          cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                              >  2:2026198/2‑101 (MQ=2)
                                          cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                              <  1:1883622/100‑1 (MQ=2)
                                          cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                              <  1:3897068/100‑1 (MQ=2)
                                           ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcg                                                             >  1:394476/1‑100 (MQ=2)
                                           ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GACGCTCCCCCCGCGGGGGCGTGGATCGAAACTGc                                                           >  1:3997235/1‑101 (MQ=255)
                                           cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                                             >  2:581921/2‑100 (MQ=2)
                                            gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaa                                                            >  1:1697012/1‑99 (MQ=2)
                                               gCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgt                                                         >  2:1833530/1‑99 (MQ=2)
                                               gCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTc                                                       <  2:4455475/101‑1 (MQ=255)
                                               gCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaactg                                                         >  1:1441731/1‑96 (MQ=2)
                                               ccTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTc                                                       <  1:1222609/100‑1 (MQ=255)
                                                cTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgta                                                        <  1:3299739/101‑4 (MQ=2)
                                                 tCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCtt                                                     >  2:62121/1‑101 (MQ=255)
                                                 tCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCtt                                                     <  2:4240799/101‑1 (MQ=255)
                                                 tCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaactgac                                                       >  1:4587757/1‑94 (MQ=2)
                                                 cccccccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtag                                                       <  1:103860/100‑5 (MQ=2)
                                                  ccccccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTg                                                    <  2:2735368/101‑1 (MQ=255)
                                                  ccccccGCAGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtaga                                                      >  2:4522955/1‑96 (MQ=2)
                                                   cccccGCGGGGGCGTGGATCGAAACATCAGCACG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGa                                                   >  2:1902119/1‑101 (MQ=255)
                                                     cccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtagatct                                                   <  1:2546810/101‑9 (MQ=2)
                                                       cgcgGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGa                                               <  1:500011/101‑1 (MQ=255)
                                                           ggggCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATggg                                           >  2:706963/1‑101 (MQ=255)
                                                            gggCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGGGATgggg                                          >  1:4436846/1‑101 (MQ=255)
                                                             ggCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGgc                                         <  1:4221397/101‑1 (MQ=255)
                                                              gCGTGGATCGAAACATCAGCTCG.TCGTTGGGCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGgcg                                        <  2:3196683/101‑1 (MQ=255)
                                                               cGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGgcgc                                       >  1:1484964/1‑101 (MQ=255)
                                                                  ggATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAg                                    <  2:2458696/101‑1 (MQ=255)
                                                                   gATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGc                                   >  1:4204726/1‑101 (MQ=255)
                                                                   gATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTAGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGc                                   <  1:4262153/101‑1 (MQ=255)
                                                                    aTCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGca                                  >  2:4348824/1‑101 (MQ=255)
                                                                       gAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACAt                               >  2:1594751/1‑101 (MQ=255)
                                                                       gAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACAt                               >  1:2226515/1‑101 (MQ=255)
                                                                        aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                              >  2:23173/1‑101 (MQ=255)
                                                                        aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCCCGTCTCGCGGGACGGGGCGCCAGCACAtt                              >  2:1780033/1‑100 (MQ=37)
                                                                         aaCACCTGCACG.TCGTTGGTCAGGCGGAC.GATGTGC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgt                             <  2:4536463/101‑1 (MQ=37)
                                                                          aCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtc                            <  2:3698749/101‑1 (MQ=255)
                                                                           cATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtcg                           >  1:875100/1‑101 (MQ=255)
                                                                            aTCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtcgt                          <  1:2175456/101‑1 (MQ=255)
                                                                             tCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGCGGATCGAAACTCCCGTCTTGAGGGATGGGGCGCGAGCACATgtcggc                         >  2:3431528/1‑99 (MQ=255)
                                                                                   cg.tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTccc                   <  2:2940058/101‑1 (MQ=255)
                                                                                   cg.tcgtTGGTCAGCCGGAC.GAGGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTccc                   >  1:2646522/1‑101 (MQ=255)
                                                                                    g.tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccc                  >  2:4145463/1‑101 (MQ=255)
                                                                                       cgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccccc                >  1:1977417/1‑101 (MQ=255)
                                                                                          tGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcgcg             >  1:3586265/1‑101 (MQ=255)
                                                                                          tGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcgcg             >  2:4094792/1‑101 (MQ=255)
                                                                                          tGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGCGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcgcg             >  1:3385921/1‑101 (MQ=255)
                                                                                           ggTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCgg            >  1:4508524/1‑101 (MQ=255)
                                                                                             tCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCgggg          <  1:2856993/101‑1 (MQ=255)
                                                                                              cAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCCCggggg         >  1:4061413/1‑101 (MQ=255)
                                                                                                gCCGGAC.GATGGCC...GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCg       <  2:24035/101‑1 (MQ=255)
                                                                                                 ccGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGt      <  1:1320554/101‑1 (MQ=255)
                                                                                                  cGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTg     <  1:919609/101‑1 (MQ=255)
                                                                                                   ggAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTgg    >  1:2983789/1‑101 (MQ=255)
                                                                                                    gAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGa   >  1:2432077/1‑101 (MQ=255)
                                                                                                     aC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGAt  <  1:2076375/101‑1 (MQ=255)
                                                                                                     |                                                                                                        
GGATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGAT  >  NC_014034/1396925‑1397125

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG