Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,698 Δ1 bp 7.8% intergenic (+510/‑421) 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,395,6980A→. 7.8% 52.1 ‑1.4 ‑0.0255intergenic (+510/‑421)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (7/13):  ref base (118/117):  total (125/130)
Fisher's exact test strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.74e-01

Only 100 of 272 total aligned reads displayed.

CCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395601‑1395798
                                                                                                   |                                                                                                         
ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTc                                                                                                        <  2:5004547/101‑1 (MQ=255)
   tGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCt                                                                                                     <  2:2908564/101‑1 (MQ=255)
       gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                                 >  2:1759941/1‑101 (MQ=255)
          tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                                              <  2:2587669/101‑1 (MQ=255)
          tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                                              <  2:1076791/101‑1 (MQ=255)
              ccTTTCGCGGGGCGAGGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                                                          <  2:4443916/100‑1 (MQ=255)
               cTTTCGCGGGGCGAGGCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                                                         <  2:4204907/101‑1 (MQ=255)
               cTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                                                         >  2:3282321/1‑101 (MQ=255)
                tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                                        >  1:2062927/1‑101 (MQ=255)
                 ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                                                       >  1:1020281/1‑101 (MQ=255)
                   cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                                                     >  1:4545643/1‑101 (MQ=255)
                    tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                                                     <  1:31130/91‑1 (MQ=2)
                    tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                                                     <  1:293838/91‑1 (MQ=2)
                     gggggCGAAGCGTCGCCCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                                                   <  1:8585/100‑1 (MQ=255)
                     gggggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                                                   <  1:1100674/100‑1 (MQ=255)
                      ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                                  >  1:2585481/1‑101 (MQ=255)
                       gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                                 >  1:4273876/1‑101 (MQ=255)
                        ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                                >  1:3128158/1‑101 (MQ=255)
                        ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCCCGGGGGGGTGGATc                                                                                >  2:3779271/1‑101 (MQ=255)
                         gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                               >  1:3793015/1‑101 (MQ=255)
                         ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                                <  2:1007073/96‑1 (MQ=2)
                         ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                                <  1:2234570/96‑1 (MQ=2)
                          cGAAGCGTCGCTCCCGCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTTGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGGTCGa                                                                              >  1:4260684/1‑101 (MQ=255)
                            aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                            <  1:2306476/101‑1 (MQ=255)
                            aaGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCGGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                            <  1:2088558/101‑1 (MQ=255)
                             gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.TACACCGGGGCTCTGGTG.ACGGTTTGG.GGCGGAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcg                                                                          >  2:3242960/4‑100 (MQ=2)
                             aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                           >  2:3435013/1‑101 (MQ=255)
                               ccgcGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                        <  2:1389154/98‑1 (MQ=255)
                               cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                        >  2:4326350/1‑101 (MQ=255)
                                gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.gg                                                                       >  2:2690386/1‑101 (MQ=255)
                                 ccGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGt                                                                      <  1:3641943/100‑1 (MQ=255)
                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACCCCCGGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcatcga                                                                      <  1:4308045/101‑6 (MQ=2)
                                     cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                  <  2:216899/100‑1 (MQ=255)
                                        ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATc                                                               <  1:340998/101‑1 (MQ=255)
                                         cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                              <  1:2805292/101‑1 (MQ=255)
                                            gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.Ggc                                                          <  1:3567444/101‑1 (MQ=255)
                                            gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.Ggc                                                          >  1:1235156/1‑101 (MQ=255)
                                             cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.Ggcg                                                         >  2:1161170/1‑101 (MQ=255)
                                               ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCa                                                       >  2:3915629/1‑101 (MQ=255)
                                               ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCa                                                       >  1:170856/1‑101 (MQ=255)
                                                 ggggTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCaaa                                                     <  1:3068081/98‑1 (MQ=255)
                                                 ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCaaa                                                     >  1:728828/1‑101 (MQ=255)
                                                  gCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAAc                                                    <  2:2212242/101‑1 (MQ=255)
                                                   cGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACg                                                   >  2:4930133/1‑101 (MQ=255)
                                                    gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACgg                                                  >  1:1530876/1‑101 (MQ=255)
                                                    gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACT.GGTATGGATCT.GGCGCAAACgg                                                  >  2:479927/1‑101 (MQ=255)
                                                     tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..a                                               >  1:1350671/1‑101 (MQ=255)
                                                      ggATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..aa                                              >  2:2285565/1‑101 (MQ=255)
                                                      ggATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..aa                                              <  2:3976818/101‑1 (MQ=255)
                                                      ggATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..aa                                              >  1:4133036/1‑101 (MQ=255)
                                                      gcctCGAAACTCCTGCATGTAGGCGGCGG.GCGGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..aa                                              <  2:678962/98‑1 (MQ=37)
                                                       gATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..aaa                                             >  2:3603901/1‑101 (MQ=255)
                                                       gATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..aaa                                             >  2:3665253/1‑101 (MQ=255)
                                                          cGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTc                                          >  2:2767760/1‑101 (MQ=255)
                                                          cGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTc                                          >  2:4527218/1‑101 (MQ=255)
                                                          cGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTc                                          >  1:516782/1‑101 (MQ=255)
                                                           gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCt                                         >  2:1896613/1‑101 (MQ=255)
                                                            aaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTg                                        <  1:89681/101‑1 (MQ=255)
                                                            aaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTg                                        <  1:220889/101‑1 (MQ=255)
                                                            aaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTg                                        >  1:2734176/1‑101 (MQ=255)
                                                             aaCATCTGCATGTGGTGGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTgt                                       <  1:2106180/101‑1 (MQ=255)
                                                              aCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTgtc                                      >  1:1928921/1‑101 (MQ=255)
                                                               cATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTgtcg                                     <  1:4909450/101‑1 (MQ=255)
                                                                aTCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTgtcgt                                    <  2:1616498/101‑1 (MQ=255)
                                                                  cTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGtcgtcg                                  >  1:3700067/1‑101 (MQ=255)
                                                                  cTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGtcgtcg                                  >  1:2398007/1‑101 (MQ=255)
                                                                   tGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGc                                 >  2:2857642/1‑101 (MQ=255)
                                                                       ggTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AACGTCTGTCGTCGCT.ccc                            <  1:4413511/100‑1 (MQ=255)
                                                                        gTAGTCGGGGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.cccc                           <  1:1161287/101‑1 (MQ=255)
                                                                        gTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.cccc                           >  1:4734729/1‑101 (MQ=255)
                                                                        gTAGTCGGCGG.GAAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.cccc                           <  2:4215985/101‑1 (MQ=255)
                                                                          ggTCGGCGG.GCAGTCCGA.GCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.cccccc                         <  2:516703/100‑1 (MQ=255)
                                                                          cggcggcgg.GCAGTTCGA.CCCGCTGTCGCTGCCCCCGCGGGGGCGTGGATCTAAACT.GGTTTGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.cccccc                         <  1:798152/98‑1 (MQ=38)
                                                                              ggcgg.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCgg                     >  2:2121547/1‑101 (MQ=255)
                                                                                 gg.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCggggg                  >  1:3318497/1‑101 (MQ=255)
                                                                                    gCAGTGCGA.GCCGCAGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCCGCGGGGGc                 <  2:2617497/101‑1 (MQ=38)
                                                                                     cAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGt               <  1:2690571/101‑1 (MQ=255)
                                                                                      aGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTg              >  2:3345035/1‑101 (MQ=255)
                                                                                      aGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTg              <  1:323472/101‑1 (MQ=255)
                                                                                       gTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTgg             <  1:3915921/101‑1 (MQ=255)
                                                                                       gTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTgg             >  2:3049997/1‑101 (MQ=255)
                                                                                         tCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGAt           <  1:3604172/101‑1 (MQ=255)
                                                                                         tCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGAt           >  1:2637204/1‑101 (MQ=255)
                                                                                         ccGA.GCCGCCGTCGCCCCCCCCGCGGGGGCGTTGATCGAAACT.GGTATGGATCT.GGCCCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGAt           <  1:2138346/100‑1 (MQ=38)
                                                                                           gA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGATCg         <  1:4511956/101‑1 (MQ=255)
                                                                                           caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTgg             <  1:4882990/94‑1 (MQ=2)
                                                                                            cagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTGGa            <  2:815732/95‑1 (MQ=2)
                                                                                            aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTGGa            <  1:2524565/95‑1 (MQ=2)
                                                                                             agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTGGAt           >  2:1988707/6‑101 (MQ=2)
                                                                                              gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGGGGTCGCT.CCCCCCGCGGGGGCGTGGATc          >  2:506593/5‑101 (MQ=2)
                                                                                              gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGATCGaa       <  2:4364384/101‑1 (MQ=255)
                                                                                               attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGGCGATCTGGGT.AAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTGGATCg         >  1:2002408/4‑101 (MQ=2)
                                                                                                ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTGGATCGa        <  2:4752545/99‑1 (MQ=2)
                                                                                                cGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAc     <  1:69988/101‑1 (MQ=255)
                                                                                                 tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTGGATCGaa       >  2:3575117/2‑101 (MQ=2)
                                                                                                 gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAcc    >  2:4111666/1‑101 (MQ=255)
                                                                                                  cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTGGATcagat      >  2:4611981/1‑97 (MQ=2)
                                                                                                  cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCt   >  1:3220649/1‑101 (MQ=255)
                                                                                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.CGCCGATCTGGGTAAGGCCGGTGGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAAc     <  1:2866293/101‑1 (MQ=2)
                                                                                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTg  <  2:1488763/101‑1 (MQ=255)
                                                                                                   |                                                                                                         
CCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.GGCGCAAACGG..AAAGTCTGTCGTCGCT.CCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395601‑1395798

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG