Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,735 T→G 3.1% intergenic (+39/+1859) RCAP_rcc01307 → / ← RCAP_rcc01309 hypothetical protein/IS66 family transposase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,7350T→G 3.1% 5.3 ‑0.4 ‑0.8160intergenic (+39/+1859)RCAP_rcc01307/RCAP_rcc01309hypothetical protein/IS66 family transposase
Reads supporting (aligned to +/- strand):  new base (2/3):  ref base (86/69):  total (88/72)
Fisher's exact test strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.64e-01

Only 100 of 163 total aligned reads displayed.

GCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGCGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGTGTCCAG  >  NC_014034/1397636‑1397806
                                                                                                     |                                                                       
gCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTc                                                                        >  2:3518898/1‑101 (MQ=255)
gCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTc                                                                        <  1:1906829/101‑1 (MQ=255)
 cTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCt                                                                       <  2:45752/101‑1 (MQ=255)
 cTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCt                                                                       <  2:2361011/101‑1 (MQ=255)
 cTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCt                                                                       <  1:3678165/101‑1 (MQ=255)
  tcgcgccGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTg                                                                      <  1:297372/96‑1 (MQ=255)
  tCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTg                                                                      <  1:3312305/101‑1 (MQ=255)
  tCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTGGTTCTGTTC.CGGTAACAGGGCGCACCCCCCGCGGGGGCGTGGATCGAAACTTGTCTg                                                                      >  1:3990477/1‑101 (MQ=38)
    cccccGCGGGGGCGTGGATCGAAACTACAGAGTGC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGAc                                                                    <  1:4435188/101‑1 (MQ=255)
     ccccGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGAcg                                                                   <  1:157802/101‑1 (MQ=255)
      cccGCGGGGGCGTGGATCGAAACGGCGTAGATCTCGGTCTGGGTGCGCTCGCGGTTGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCtatggc                                                                  <  1:3366579/101‑6 (MQ=0)
       ccGCGGGGGCGTGGATCGAAACGGCGTAGATCTCGGTCTGGGTGCGCTCGCGGTTGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCtatggcg                                                                 <  1:341240/101‑7 (MQ=0)
       ccGCGGGGGCGTGGATCGAAACGGCGTAGATCTCGGTCTGGGTGCGCTCGCGGTTGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCtatggcg                                                                 <  2:1821775/101‑7 (MQ=0)
       ccGCGGGGGCGTGGATCGAAACGGCGTAGATCTCGGTCTGGGTGCGCTCGCGGTTGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCtatggcg                                                                 >  1:1082998/1‑95 (MQ=0)
        cgcgGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACgcgc                                                                <  1:230528/101‑1 (MQ=255)
          cgGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtc                                                              >  2:2650603/1‑101 (MQ=255)
           gggggCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtca                                                             >  1:2763970/1‑101 (MQ=255)
           gggggCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtca                                                             <  2:3823229/101‑1 (MQ=255)
           gggggCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtca                                                             <  1:3924987/101‑1 (MQ=255)
            gtggCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtcat                                                            <  1:1517329/99‑1 (MQ=255)
            ggggCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCtcat                                                            <  1:730189/101‑1 (MQ=255)
               gCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCCGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAg                                                         >  1:3252673/1‑101 (MQ=255)
                gggggttcaaaaCTCCGGAGTGC.GGTTCTTTTTCTGTTC.CGGTAAAAGAT.GCCCCCCCCCGGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCa                                                       <  1:4589841/92‑1 (MQ=255)
                cGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGc                                                        >  2:3864592/1‑101 (MQ=255)
                 gTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCa                                                       >  2:2166038/1‑101 (MQ=255)
                  tGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCAc                                                      >  1:948296/1‑101 (MQ=255)
                  tGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCAc                                                      >  1:4309724/1‑101 (MQ=255)
                  tGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCAc                                                      >  1:1503730/1‑101 (MQ=255)
                  tGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCAc                                                      >  2:4104400/1‑101 (MQ=255)
                  tGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCAc                                                      >  2:4232192/1‑101 (MQ=255)
                  gggATCGCAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCAc                                                      <  2:4556224/100‑1 (MQ=255)
                   ggATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACg                                                     <  2:853383/101‑1 (MQ=255)
                   ggATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACg                                                     >  2:1443938/1‑101 (MQ=255)
                       cGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCaa                                                 >  2:689129/1‑101 (MQ=255)
                        gAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAg                                                <  1:938455/101‑1 (MQ=255)
                        gAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGATCCCCCCTCGGGGGCTTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAg                                                >  1:4636984/1‑101 (MQ=255)
                         aaaCTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGa                                               >  1:4477491/1‑101 (MQ=255)
                          aaCTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGAt                                              <  1:3666844/101‑1 (MQ=255)
                          aaCTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGAt                                              >  1:1868443/1‑101 (MQ=255)
                          aaCTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGAt                                              >  1:2536817/1‑101 (MQ=255)
                          aaCTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGAt                                              <  2:3362166/101‑1 (MQ=255)
                          aaCTACAAAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCCCGTCAAGAt                                              >  1:3149339/1‑101 (MQ=255)
                           aCTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGAtt                                             <  1:1045572/101‑1 (MQ=255)
                           aCTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGAtt                                             >  2:3112680/1‑101 (MQ=255)
                            cTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTc                                            >  2:3007739/1‑101 (MQ=255)
                             tACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCa                                           >  1:3517479/1‑101 (MQ=255)
                               cAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGt                                         >  1:2054227/1‑101 (MQ=255)
                                agagTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGGCAAGATTCAGTc                                        >  1:1045386/1‑101 (MQ=255)
                                 gagTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCg                                       <  2:3944542/101‑1 (MQ=255)
                                 gagTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCg                                       <  2:2448390/101‑1 (MQ=255)
                                  agTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGc                                      >  2:1029787/1‑101 (MQ=255)
                                    ttC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGGGGATCGACACTTGTCTGACGCGCTCATCAGAAGGTCAAGATTCAGTCGCTc                                    >  1:4416581/1‑101 (MQ=38)
                                     tC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcc                                   <  1:2631034/101‑1 (MQ=255)
                                         ggTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTccccc                                >  2:1211778/1‑101 (MQ=255)
                                          gTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcccccc                               <  1:2559895/101‑1 (MQ=255)
                                          gTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcccccc                               <  2:2633324/101‑1 (MQ=255)
                                          gTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcccccc                               <  1:1776451/101‑1 (MQ=255)
                                          gTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcccccc                               <  2:1487871/101‑1 (MQ=255)
                                          gTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTcccccc                               <  2:1632664/101‑1 (MQ=255)
                                            cTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCcgc                             >  2:2282421/1‑101 (MQ=255)
                                            cTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCcgc                             >  2:1871977/1‑101 (MQ=255)
                                              tgttctgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCgg                           >  2:3209398/1‑101 (MQ=255)
                                               gttctgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCggg                          >  1:2760447/1‑101 (MQ=255)
                                                ttctgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCgggg                         >  2:731241/1‑101 (MQ=255)
                                                ttctgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCgggg                         >  1:882432/1‑101 (MQ=255)
                                                 tctgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCggggg                        >  1:395013/1‑101 (MQ=255)
                                                 tctcttc.CGGTAACAGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCggggg                        <  1:1939773/101‑1 (MQ=255)
                                                  ctgttc.CGGTAACGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGc                       <  2:666927/101‑1 (MQ=255)
                                                  ctgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGc                       <  2:1805646/101‑1 (MQ=255)
                                                  ctgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGc                       <  2:132736/101‑1 (MQ=255)
                                                  ctgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGc                       <  2:82905/101‑1 (MQ=255)
                                                   tgttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCg                      >  2:3368287/1‑101 (MQ=255)
                                                    gttc.CGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGt                     >  1:1683540/1‑101 (MQ=255)
                                                       c.cGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGa                  <  1:2105088/101‑1 (MQ=255)
                                                       c.cGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGa                  >  1:1076447/1‑101 (MQ=255)
                                                         cGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATGGAAACTTGTCTGACGCGCTCATCAGCAAGTCAAGATTAAGTCGCTCCCCCCCCGGGGGCGTGGAt                 >  1:3831785/1‑101 (MQ=38)
                                                         cGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGAt                 <  2:3149095/101‑1 (MQ=255)
                                                         cGGTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGAt                 >  2:779625/1‑101 (MQ=255)
                                                          ggTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATc                <  1:3915071/101‑1 (MQ=255)
                                                           gTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCg               >  2:2800329/1‑101 (MQ=255)
                                                           gTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCg               <  1:2074519/101‑1 (MQ=255)
                                                           gTAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCg               >  1:208338/1‑101 (MQ=255)
                                                            tAACAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa              >  1:667181/1‑101 (MQ=255)
                                                             aaCAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa             <  2:1982229/101‑1 (MQ=255)
                                                               cAGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc           <  2:948216/101‑1 (MQ=255)
                                                                aGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct          >  2:3614685/1‑100 (MQ=255)
                                                                aGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct          >  2:3990610/1‑100 (MQ=255)
                                                                aGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct          >  1:3405207/1‑100 (MQ=255)
                                                                aGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct          >  2:1519604/1‑100 (MQ=255)
                                                                 ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc         >  2:1112222/1‑99 (MQ=255)
                                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcg        >  2:1557959/1‑98 (MQ=255)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTTGCGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCCCGGGGGCGTGGATCGaaagtcggc      >  2:1023900/1‑95 (MQ=37)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggca     >  2:2705131/1‑95 (MQ=255)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggca     >  2:72617/1‑95 (MQ=255)
                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCT.CCCCCGCGGGGGCGTGGATCGAAAcccggcaag   >  1:1449185/1‑93 (MQ=37)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcaag   <  1:667059/101‑9 (MQ=255)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcaag   >  2:1179491/1‑93 (MQ=255)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcaag   <  2:655568/101‑9 (MQ=255)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcaag   >  2:3696373/1‑93 (MQ=255)
                                                                        ccccccGCGGGGGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcaagg  <  1:3112924/101‑10 (MQ=38)
                                                                                                     |                                                                       
GCTCCCCCCGCGGGGGCGTGGATCGAAACTACAGAGTTC.CGGTCTTGTTCTGTTC.CGGTAACAGGTCGCTCCCCCCGCGGGCGCGTGGATCGAAACTTGTCTGACGCGCTCATCAGCACGTCAAGATTCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGTGTCCAG  >  NC_014034/1397636‑1397806

Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG