Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,015 T→A 12.7% intergenic (+827/‑104) RCAP_rcc01306 → / → RCAP_rcc01307 Cas2 family CRISPR‑associated protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,0150T→A12.7% 85.5 ‑0.4 ‑0.4236intergenic (+827/‑104)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/17):  ref base (98/107):  total (111/125)
Fisher's exact test strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.72e-01

Only 100 of 257 total aligned reads displayed.

GGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGG  >  NC_014034/1395915‑1396115
                                                                                                    |                                                                                                        
gggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCtt                                                                                                          >  1:3462238/1‑101 (MQ=255)
 ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTg                                                                                                         <  2:1596489/101‑1 (MQ=255)
 ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTg                                                                                                         >  1:936474/1‑101 (MQ=255)
 ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTg                                                                                                         <  2:1470601/101‑1 (MQ=255)
 ggCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTAGCTCCCCCAGAGGGGGCATGGATCGAAAATCGCCCTACATCAGCGCTTg                                                                                                         >  1:1532505/1‑101 (MQ=255)
  gCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGc                                                                                                        >  1:1735903/1‑101 (MQ=255)
   cGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGcc                                                                                                       >  1:1802289/1‑101 (MQ=255)
     tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                                                                                     >  1:1210855/1‑101 (MQ=255)
     tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                                                                                     <  1:2701909/101‑1 (MQ=255)
     tGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTc                                                                                                     >  2:3864342/1‑101 (MQ=255)
      ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.a                                                                                                   <  1:858191/101‑1 (MQ=255)
      ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.a                                                                                                   <  2:28980/101‑1 (MQ=255)
      ggATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.a                                                                                                   <  1:782550/101‑1 (MQ=255)
       gATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.Ac                                                                                                  >  2:761121/1‑101 (MQ=255)
        aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACg                                                                                                 >  2:3753561/1‑101 (MQ=255)
        aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACg                                                                                                 <  1:3718506/101‑1 (MQ=255)
        aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACg                                                                                                 <  1:4382238/101‑1 (MQ=255)
           gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCt                                                                                              <  2:2873276/101‑1 (MQ=255)
           gAAAAGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.GCGCCt                                                                                              <  1:3435355/101‑1 (MQ=255)
            aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTg                                                                                             >  2:4469514/1‑101 (MQ=255)
            aaaCGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTg                                                                                             >  2:2609507/1‑101 (MQ=255)
                  tGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..Gt                                                                                    <  2:2335059/101‑1 (MQ=255)
                    cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCg                                                                                  <  2:1699253/101‑1 (MQ=255)
                      ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCt                                                                                <  1:3205247/101‑1 (MQ=255)
                           ggCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCNCTTGCCTC.ACGCCTGC.CGC..GTCGCTccccc                                                                           >  1:551177/1‑101 (MQ=255)
                             ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCcg                                                                         <  2:522928/101‑1 (MQ=255)
                              cTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCcgc                                                                        >  1:4187512/1‑101 (MQ=255)
                               tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCcgcg                                                                       <  1:2550393/101‑1 (MQ=255)
                               tACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCcgcg                                                                       >  2:2652065/1‑101 (MQ=255)
                                aCGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCgg                                                                      >  1:3145741/1‑101 (MQ=255)
                                     gggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCg                                                                 >  1:2970081/1‑101 (MQ=255)
                                      ggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGt                                                                <  2:4563317/101‑1 (MQ=255)
                                      ggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGt                                                                <  2:1007737/101‑1 (MQ=255)
                                        gaGTGGCCTGGTCGCCCCCCCCGCGGGGGCNTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg                                                              <  1:1371432/100‑1 (MQ=255)
                                        cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg                                                              <  2:215747/101‑1 (MQ=255)
                                        cAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTgg                                                              >  2:3148284/1‑101 (MQ=255)
                                           tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                           >  1:2884199/1‑101 (MQ=255)
                                           tGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATc                                                           <  1:1288211/101‑1 (MQ=255)
                                            ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                          >  1:2109726/1‑101 (MQ=255)
                                             gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC..GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                          <  2:1989406/97‑1 (MQ=9)
                                               gcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC..GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                        >  2:935276/3‑101 (MQ=2)
                                               cTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                       >  1:4055892/1‑101 (MQ=255)
                                                tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                      >  2:3553173/1‑101 (MQ=255)
                                                tGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                      >  1:1819565/1‑101 (MQ=255)
                                                 ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC..GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                      <  1:2259283/101‑1 (MQ=2)
                                                 ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                     >  2:3914623/1‑101 (MQ=255)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC..GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                     >  2:4120081/1‑100 (MQ=2)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                    <  1:2840967/101‑1 (MQ=255)
                                                  gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                    <  1:653094/101‑1 (MQ=255)
                                                  gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                    <  1:999350/101‑1 (MQ=255)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcgc                                                >  1:4082608/1‑101 (MQ=255)
                                                       tgcgcccGCGGGGGCGTGTATCGAAACTCGGCCTCGATCCG.GTCAGCCTCGTAGTGCGCGCTG..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                <  2:1095472/97‑7 (MQ=2)
                                                       tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC..GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                <  2:1851575/101‑7 (MQ=2)
                                                       tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                               <  2:1262980/101‑1 (MQ=255)
                                                        ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                              <  2:1917427/101‑1 (MQ=255)
                                                        ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                              >  1:3396542/1‑101 (MQ=255)
                                                          ccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC..GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacc                                             <  2:1370798/101‑10 (MQ=2)
                                                          ccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC..GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacc                                             >  2:3971098/1‑92 (MQ=2)
                                                          ccccGCGGGGGCGTGGATCGAAACTCGCCCTAGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCt                                            >  2:626216/1‑101 (MQ=255)
                                                              gcgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATg                                        <  2:2084258/101‑1 (MQ=255)
                                                              gcgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATg                                        >  2:4370404/1‑101 (MQ=255)
                                                               cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGt                                       >  1:1261568/1‑101 (MQ=255)
                                                                gggggCGTGGATGGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                      >  1:2531607/1‑101 (MQ=255)
                                                                gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                      >  1:4494721/1‑101 (MQ=255)
                                                                 ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                     >  1:1646869/1‑101 (MQ=255)
                                                                  gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..TTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                    <  1:150192/101‑1 (MQ=255)
                                                                  gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                    <  2:2305256/101‑1 (MQ=255)
                                                                  gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                    <  1:4599816/101‑1 (MQ=255)
                                                                  gggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtata                                    <  1:261038/101‑1 (MQ=255)
                                                                   ggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTatat                                   >  1:3816866/1‑101 (MQ=255)
                                                                    gCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATg                                  <  1:2191316/101‑1 (MQ=255)
                                                                     cGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGc                                 >  2:3342787/1‑101 (MQ=255)
                                                                      gTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGcc                                >  2:254650/1‑101 (MQ=255)
                                                                      gTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGcc                                >  2:721165/1‑101 (MQ=255)
                                                                       tGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGccc                               <  2:215449/101‑1 (MQ=255)
                                                                        ggATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCg                              >  2:3393876/1‑101 (MQ=255)
                                                                           tCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTg                           <  1:2005719/101‑1 (MQ=255)
                                                                             gAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                         >  1:2499203/1‑101 (MQ=255)
                                                                             gAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                         <  1:3955797/101‑1 (MQ=255)
                                                                                 cTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACgg                     >  2:3330220/1‑101 (MQ=255)
                                                                                  tCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGc                    <  2:1163960/101‑1 (MQ=255)
                                                                                   cGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCt                   <  2:1646758/101‑1 (MQ=255)
                                                                                   cGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCt                   <  1:3587646/101‑1 (MQ=255)
                                                                                      ccTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTc                >  2:4349471/1‑101 (MQ=255)
                                                                                       cTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCg               <  1:2277716/101‑1 (MQ=255)
                                                                                        tCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGc              <  1:3317705/101‑1 (MQ=255)
                                                                                         cGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt             >  2:3069904/1‑101 (MQ=255)
                                                                                         cGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt             >  1:2161361/1‑101 (MQ=255)
                                                                                         cGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt             >  1:2220908/1‑101 (MQ=255)
                                                                                          gATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTc            >  1:207347/1‑101 (MQ=255)
                                                                                           aTCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcc           <  1:1252636/101‑1 (MQ=255)
                                                                                             ccGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccc         >  2:2749523/1‑101 (MQ=255)
                                                                                              cgcgCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccccc        <  1:4471705/101‑1 (MQ=255)
                                                                                              cgcgCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccccc        >  1:4359143/1‑101 (MQ=255)
                                                                                               gcgcTTGCCTC.CCGCCTAC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc       <  2:4187224/101‑1 (MQ=255)
                                                                                               gcgcTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc       <  2:4055606/101‑1 (MQ=255)
                                                                                               gcgcTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc       <  2:2969852/101‑1 (MQ=255)
                                                                                                 gcTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcgc     <  1:1042640/101‑1 (MQ=255)
                                                                                                   ttGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCgg   >  2:4507021/1‑101 (MQ=255)
                                                                                                   ttGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCgg   <  2:2514235/101‑1 (MQ=255)
                                                                                                    |                                                                                                        
GGGCGTGGATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTC.ACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGG  >  NC_014034/1395915‑1396115

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