Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,165 C→G 4.9% R349R (CGC→CGG RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,1650C→G 4.9% 18.9 ‑4.0 ‑0.1244R349R (CGC→CGGRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/10):  ref base (121/111):  total (123/121)
Fisher's exact test strand distribution p-value = 1.87e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.80e-01

Only 100 of 270 total aligned reads displayed.

CTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC  >  NC_014034/1397065‑1397265
                                                                                                     |                                                                                                    
cTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTgcgc                                                                                                      <  2:3568065/101‑1 (MQ=255)
 tGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGcc                                                                                                     <  1:2568881/101‑1 (MQ=255)
  gCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCg                                                                                                    >  2:4424280/1‑101 (MQ=255)
       cTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCt                                                                                               >  1:5205190/1‑101 (MQ=255)
       cTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GGGCGCCGTCGcc                                                                                               >  2:4705124/1‑100 (MQ=255)
        ttGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTc                                                                                              >  1:4793348/1‑101 (MQ=255)
        ttGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTc                                                                                              >  2:3757710/1‑101 (MQ=255)
            gTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTccccc                                                                                          <  1:5017990/101‑1 (MQ=255)
             tGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTcccccc                                                                                         >  1:3154092/1‑101 (MQ=255)
             tGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTcccccc                                                                                         <  1:5726106/101‑1 (MQ=255)
              gATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCcg                                                                                        <  2:5773229/101‑1 (MQ=255)
              gATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCcg                                                                                        >  1:5169680/1‑101 (MQ=255)
              gATGGGGCGCGAGCACATGTCGACGCCGCCCCCGGGGGGGCGTGGATCGAAACGGCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCcg                                                                                        <  2:1268944/101‑1 (MQ=38)
                    gcgcGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCggggg                                                                                  <  2:3329736/101‑1 (MQ=255)
                    gcgcGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCggggg                                                                                  >  1:4774396/1‑101 (MQ=255)
                       ggAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGt                                                                               <  2:2193607/100‑1 (MQ=255)
                       cgAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCgg                                                                               >  1:2624732/1‑100 (MQ=255)
                        gAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTg                                                                              >  2:2848340/1‑101 (MQ=255)
                        gAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTg                                                                              <  2:5033895/101‑1 (MQ=255)
                         aGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGTGCGTgg                                                                             >  2:2824840/1‑101 (MQ=255)
                           cggctgcggtcgCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACG.CGGCAACCGGCCCGGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                           <  1:694555/93‑1 (MQ=0)
                            acaTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGGTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                          <  2:3517850/101‑1 (MQ=255)
                            acaTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                          <  2:1152893/101‑1 (MQ=255)
                              ctgcggtcgCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACG.CGGCAACCGGCCCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                        <  2:184872/96‑1 (MQ=2)
                               tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                       >  1:5457366/1‑101 (MQ=255)
                                 cggtcgCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACG.CGGCAACCGGCCCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                     >  2:3055741/3‑101 (MQ=2)
                                  cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaaag                                                                    <  1:4242528/101‑3 (MQ=255)
                                  cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg                                                                    <  1:3605098/101‑1 (MQ=255)
                                  cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg                                                                    >  2:4598578/1‑101 (MQ=255)
                                    tcgCTCCCGCCGCGGGGGCGTGGATCGAAAAGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACgtg                                                                  <  2:4791062/101‑1 (MQ=255)
                                     cgCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACgtgt                                                                 >  2:816497/1‑101 (MQ=255)
                                     cgCGCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACgtgt                                                                 <  2:1860693/101‑1 (MQ=255)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTc                                                                >  2:5798820/1‑101 (MQ=255)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTc                                                                <  1:1649528/101‑1 (MQ=255)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCt                                                               >  1:1976989/1‑101 (MQ=255)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCt                                                               >  1:5073593/1‑101 (MQ=255)
                                        tCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTg                                                              >  2:2748986/1‑101 (MQ=255)
                                        tCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTg                                                              <  1:2434014/101‑1 (MQ=255)
                                        cccccccGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTg                                                              <  2:2699523/100‑1 (MQ=255)
                                          cccccGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGcc                                                            >  1:238926/1‑101 (MQ=255)
                                          cccccGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGcc                                                            >  1:4543791/1‑101 (MQ=255)
                                           ccccGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCa                                                           <  2:5267360/101‑1 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCac                                                          <  2:2580727/101‑1 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCac                                                          >  1:3212384/1‑101 (MQ=255)
                                              cgcgGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCacga                                                        >  2:3837009/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGc                                                     >  2:5450233/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGc                                                     >  2:1970983/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGc                                                     >  2:5720392/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGc                                                     >  2:1563953/1‑101 (MQ=255)
                                                  ggggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGcc                                                    >  2:3718165/1‑101 (MQ=255)
                                                   gggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGccc                                                   <  1:5392722/101‑1 (MQ=255)
                                                   gggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGccc                                                   <  1:4354509/101‑1 (MQ=255)
                                                   gggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGccc                                                   >  2:1025994/1‑101 (MQ=255)
                                                   gggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GGGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGccc                                                   <  1:4868778/101‑1 (MQ=255)
                                                    ggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCt                                                  <  2:5563948/101‑1 (MQ=255)
                                                    ggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCt                                                  <  2:5135846/101‑1 (MQ=255)
                                                    ggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCt                                                  >  2:4952151/1‑101 (MQ=255)
                                                    ggCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCt                                                  <  1:4295234/101‑1 (MQ=255)
                                                      cGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGc                                                >  2:3951898/1‑101 (MQ=255)
                                                      cGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGc                                                >  1:550717/1‑101 (MQ=255)
                                                      cGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGc                                                >  2:5262754/1‑101 (MQ=255)
                                                        tGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCAc                                              >  1:1417566/1‑101 (MQ=255)
                                                         ggATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCAcc                                             >  2:142519/1‑101 (MQ=255)
                                                          gATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCa                                            >  2:1919553/1‑101 (MQ=255)
                                                          gATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCa                                            >  1:3845603/1‑101 (MQ=255)
                                                          gATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCa                                            >  1:4909791/1‑101 (MQ=255)
                                                            tCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGa                                          <  1:2217545/101‑1 (MQ=255)
                                                             cGAAACGTCCAGCTCGCGCCACATGTCTTCGCCC.GTGCGCCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGAt                                         <  1:3945448/101‑1 (MQ=255)
                                                               aaaCGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATcc                                       <  2:3153903/101‑1 (MQ=255)
                                                               aaaCGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATcc                                       >  2:758897/1‑101 (MQ=255)
                                                                aaCGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCcg                                      <  2:1051989/101‑1 (MQ=255)
                                                                     ccAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCg                                 >  2:4619573/1‑101 (MQ=255)
                                                                     ccAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCg                                 <  1:2749147/101‑1 (MQ=255)
                                                                          tCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCt                            >  2:3450100/1‑101 (MQ=255)
                                                                           cgcgCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTc                           <  1:3725027/101‑1 (MQ=255)
                                                                            gcgcCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTcc                          <  2:182194/101‑1 (MQ=255)
                                                                                 acaTGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCcg                     <  1:126510/101‑1 (MQ=255)
                                                                                 acaTGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCcg                     <  1:4581826/101‑1 (MQ=255)
                                                                                   aTGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCcgcg                   >  2:4013935/1‑101 (MQ=255)
                                                                                   aTGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCcgcg                   >  2:4062058/1‑101 (MQ=255)
                                                                                   aTGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGCCGCTCCCCCcgcg                   >  1:7571/1‑101 (MQ=255)
                                                                                      tCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCgggg                <  1:2087202/101‑1 (MQ=255)
                                                                                      tCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCgggg                >  2:2270507/1‑101 (MQ=255)
                                                                                       cTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCggggg               >  1:1328427/1‑101 (MQ=255)
                                                                                       cTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCggggg               <  1:1026069/101‑1 (MQ=255)
                                                                                       cTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCggggg               >  1:1956111/1‑101 (MQ=255)
                                                                                         tCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGCCGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCg             >  1:3454148/1‑101 (MQ=255)
                                                                                         tCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCg             >  2:169523/1‑101 (MQ=255)
                                                                                         tCGACC.GCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTTCCCCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCg             <  1:809066/101‑1 (MQ=255)
                                                                                          cGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGt            <  1:3952142/101‑1 (MQ=255)
                                                                                           gACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTg           >  2:218013/1‑101 (MQ=255)
                                                                                            gcc.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTgg          <  1:2865246/100‑1 (MQ=255)
                                                                                              c.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGAt        <  2:4543519/101‑1 (MQ=255)
                                                                                                 tGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCg      >  2:542216/1‑101 (MQ=255)
                                                                                                  gcgcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGa     <  2:4272988/101‑1 (MQ=255)
                                                                                                    gcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   >  1:1842094/1‑101 (MQ=255)
                                                                                                    gcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   <  1:998230/101‑1 (MQ=255)
                                                                                                    gcCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa   <  1:5263044/101‑1 (MQ=255)
                                                                                                     ccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac  <  2:4116134/101‑1 (MQ=255)
                                                                                                     ccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac  >  1:3013265/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
CTGCCGTCTTGAGTGATGGGGCGCGAGCACATGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCCAGCTCGCGCCACATGTCTTCGACC.GTGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTGTCTGCCACGACGCCCTGCACCAGATCCGCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC  >  NC_014034/1397065‑1397265

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