Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,027 G→C 16.3% T303T (ACG→TCC RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,0270G→C16.3% 97.4 ‑2.1 ‑0.4275T303T (ACG→ACCRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (17/22):  ref base (122/90):  total (161/114)
Fisher's exact test strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.43e-01

Only 100 of 291 total aligned reads displayed.

ATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCG  >  NC_014034/1396927‑1397127
                                                                                                      |                                                                                                       
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)
    aaaCAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGt                                                                                                     >  1:3480855/1‑101 (MQ=255)
     aaCAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAATCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTc                                                                                                    >  1:4036891/1‑101 (MQ=255)
     aaCAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTc                                                                                                    >  2:132336/1‑101 (MQ=255)
       cAGCGTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...g                                                                                               <  2:2985757/101‑1 (MQ=255)
           cTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGc                                                                                           >  2:461799/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:1638261/1‑101 (MQ=255)
                     gCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCcgcg                                                                                 <  2:756917/101‑1 (MQ=255)
                     gCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCcgcg                                                                                 >  2:2280191/1‑101 (MQ=255)
                     gCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCcgcg                                                                                 >  1:2214320/1‑101 (MQ=255)
                       aTCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCggg                                                                               >  2:3046191/1‑101 (MQ=255)
                        tCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCgggg                                                                              <  1:4161544/101‑1 (MQ=255)
                         cgccTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCggggg                                                                             <  2:1688783/98‑1 (MQ=255)
                          ggCTTCCGCCTGCGTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGc                                                                            <  2:2023416/101‑1 (MQ=255)
                           gCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCg                                                                           <  2:855754/101‑1 (MQ=255)
                             ttCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTg                                                                         >  2:931131/1‑101 (MQ=255)
                              tCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGGgg                                                                        >  1:4574773/1‑101 (MQ=255)
                                gatcccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTgg                                                                        >  1:2302473/10‑101 (MQ=2)
                                cGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      <  1:3148972/101‑1 (MQ=255)
                                 gCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                     <  1:1119958/101‑1 (MQ=255)
                                 atcccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGa                                                                       >  1:1601526/9‑101 (MQ=2)
                                  tcccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                      <  1:2154770/94‑1 (MQ=2)
                                    ccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                    <  2:2653565/96‑1 (MQ=2)
                                    agcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                    >  1:1608934/7‑101 (MQ=2)
                                    agcggcgTCGCCCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                    <  1:642706/95‑1 (MQ=2)
                                     gCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                 <  2:375360/101‑1 (MQ=255)
                                       ggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                 >  1:467300/4‑101 (MQ=0)
                                       gTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                               >  1:536849/1‑101 (MQ=255)
                                       ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                 <  1:1034040/99‑1 (MQ=2)
                                       ccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                 >  1:3243559/3‑101 (MQ=2)
                                        tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                              >  1:3686622/1‑101 (MQ=255)
                                        cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                <  1:1883622/100‑1 (MQ=2)
                                        cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                >  1:1380751/2‑101 (MQ=2)
                                         ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcg                                                               >  1:1834991/1‑100 (MQ=2)
                                         cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                                               >  1:80515/2‑100 (MQ=2)
                                         cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                                               >  2:581921/2‑100 (MQ=2)
                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaa                                                              >  1:1697012/1‑99 (MQ=2)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgt                                                           >  1:4491855/1‑99 (MQ=2)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaactg                                                           >  1:1441731/1‑96 (MQ=2)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAAC.CCGCCACG.TCTTGGAACAGCGTGCC.ATTGGCCCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaacta                                                           >  2:1279441/1‑96 (MQ=2)
                                             ccTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTc                                                         <  1:1222609/100‑1 (MQ=255)
                                              cTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgta                                                          <  1:3299739/101‑4 (MQ=2)
                                              cTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCt                                                        >  2:2844597/1‑101 (MQ=255)
                                               tCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCtt                                                       <  2:4240799/101‑1 (MQ=255)
                                               tCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCtt                                                       >  2:62121/1‑101 (MQ=255)
                                               cccccccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtag                                                         <  1:103860/100‑5 (MQ=2)
                                                ccccccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtaga                                                        >  1:3931295/1‑96 (MQ=2)
                                                ccccccGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtaga                                                        <  1:2719672/101‑6 (MQ=2)
                                                ccccccGCAGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGcgtaga                                                        >  2:4522955/1‑96 (MQ=2)
                                                  cgccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAg                                                    >  2:1880572/3‑101 (MQ=255)
                                                   cccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGt                                                   >  1:2566625/1‑101 (MQ=255)
                                                    ccGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTg                                                  >  1:2249181/1‑101 (MQ=255)
                                                         ggggCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATggg                                             >  2:2418484/1‑101 (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)
                                                             ggTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGgcgc                                         <  2:3602927/100‑1 (MQ=255)
                                                             cGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGgcgc                                         >  1:1484964/1‑101 (MQ=255)
                                                              gTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGcgcg                                        >  2:1250105/1‑101 (MQ=255)
                                                                 gATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGc                                     <  1:3851425/101‑1 (MQ=255)
                                                                 gATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTAGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGc                                     <  1:4262153/101‑1 (MQ=255)
                                                                  aTCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGca                                    >  2:4348824/1‑101 (MQ=255)
                                                                      aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                                >  1:2658734/1‑101 (MQ=255)
                                                                      aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                                <  1:2418648/101‑1 (MQ=255)
                                                                      aaaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATg                                >  2:23173/1‑101 (MQ=255)
                                                                       aaCATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgt                               >  2:2725228/1‑101 (MQ=255)
                                                                       aaCACCTGCACG.TCGTTGGTCAGGCGGAC.GATGTGC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgt                               <  2:4536463/101‑1 (MQ=37)
                                                                         cATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtcg                             >  1:875100/1‑101 (MQ=255)
                                                                          aTCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATgtcgc                            >  1:2092834/1‑100 (MQ=255)
                                                                           tCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGCGGATCGAAACTCCCGTCTTGAGGGATGGGGCGCGAGCACATgtcggc                           >  2:3431528/1‑99 (MQ=255)
                                                                            cAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGGGATGGGGCGCGAGCACATgtcgtgc                          >  2:314530/1‑99 (MQ=255)
                                                                              gCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCt                        <  2:1083881/101‑1 (MQ=255)
                                                                                tcg.tcgTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcc                      <  1:2393973/101‑1 (MQ=255)
                                                                                 cg.tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTccc                     >  2:3694417/1‑101 (MQ=255)
                                                                                  g.tcgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccc                    >  2:4145463/1‑101 (MQ=255)
                                                                                     cgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccccc                  <  2:4448888/101‑1 (MQ=255)
                                                                                     cgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccccc                  >  1:1977417/1‑101 (MQ=255)
                                                                                     cgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccccc                  <  2:180699/101‑1 (MQ=255)
                                                                                     cgtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTcccccc                  <  1:2538913/101‑1 (MQ=255)
                                                                                      gtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcg                 >  2:4468523/1‑101 (MQ=255)
                                                                                      gtTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcg                 <  1:4062906/101‑1 (MQ=255)
                                                                                        tGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCcgcg               >  2:4094792/1‑101 (MQ=255)
                                                                                           tCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCgggg            <  1:2856993/101‑1 (MQ=255)
                                                                                            cAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCCCcgggg           >  2:1939638/1‑101 (MQ=255)
                                                                                               ccGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGt        <  1:1320554/101‑1 (MQ=255)
                                                                                                cGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTg       <  2:2214158/101‑1 (MQ=255)
                                                                                                cGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTg       <  1:919609/101‑1 (MQ=255)
                                                                                                cGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGAGATGGGGCGCGAGCGCATGTCGTCGTTCCCCCCGCCGGGGCgcg       >  2:2096918/1‑99 (MQ=37)
                                                                                                 ggAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTgg      >  1:2983789/1‑101 (MQ=255)
                                                                                                 ggAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTgg      <  1:3539652/101‑1 (MQ=255)
                                                                                                 ggAC.GATCTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTgg      >  2:144544/1‑101 (MQ=255)
                                                                                                  gAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGa     >  1:3967068/1‑101 (MQ=255)
                                                                                                    c.GATGTCC...GTCGCTCCCNCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATc   <  2:2169907/101‑1 (MQ=255)
                                                                                                    c.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATc   >  1:4188968/1‑101 (MQ=255)
                                                                                                      gATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg  <  2:2389943/101‑1 (MQ=255)
                                                                                                      |                                                                                                       
ATCGAAACAGCCTGCGCGGCCGCATCGGCTTCCGCCTGCGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCC...GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCG  >  NC_014034/1396927‑1397127

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG