Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,820,393 +C 2.3% intergenic (+84/‑102) RCAP_rcc01683 → / → RCAP_rcc01684 terminase‑like family protein/HK97 family phage portal protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,820,3932.→C 2.3% 2.3 ‑4.5 ‑0.1132intergenic (+84/‑102)RCAP_rcc01683/RCAP_rcc01684terminase‑like family protein/HK97 family phage portal protein
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (102/27):  total (102/30)
Fisher's exact test strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.56e-01

Only 100 of 163 total aligned reads displayed.

AGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACGG  >  NC_014034/1820293‑1820493
                                                                                                       |                                                                                                     
aGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGGGTGaaa                                                                                                        >  1:4077496/1‑101 (MQ=255)
 ggTGCGCGGGCTTTGAGCGCGCAAGGGCGTACTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..c                                                                                                     >  1:4130809/1‑101 (MQ=255)
  gTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..cc                                                                                                    >  2:638404/1‑101 (MQ=255)
  gTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..cc                                                                                                    >  1:4157130/1‑101 (MQ=255)
  gTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..cc                                                                                                    >  1:809859/1‑101 (MQ=255)
   tGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTTCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..ccc                                                                                                   <  2:3108460/101‑1 (MQ=255)
      gcgGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCg                                                                                                >  2:2613064/1‑101 (MQ=255)
       cgGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGc                                                                                               >  1:2251519/1‑101 (MQ=255)
       cgGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..AGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGc                                                                                               >  2:3260670/1‑101 (MQ=255)
        gggCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCt                                                                                              >  1:1385695/1‑101 (MQ=255)
        gggCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCt                                                                                              >  2:3182940/1‑101 (MQ=255)
         ggCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCtt                                                                                             <  2:3800063/101‑1 (MQ=255)
          gCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTc                                                                                            <  1:2097880/101‑1 (MQ=255)
           cTTTGAGCGCGCAAGGGCTTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTAAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCt                                                                                           <  2:2239545/101‑1 (MQ=255)
             ttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGc                                                                                         >  1:2811991/1‑101 (MQ=255)
             ttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGc                                                                                         >  2:1830340/1‑101 (MQ=255)
              tGAGCGCGCAAGGCCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGTGGGGGTTAAA..CCCCCGCTTCtctc                                                                                        >  1:2271538/1‑98 (MQ=25)
              ggggcgcgaAGGGGCTTTTTGGGACAAGAAAAAAGCTGCTTTTTTTTTTTTGGTAAAAATC..CGCCCCGGGGGTGCGGGGGTGTGAAC..CCCCCGCTTCTGcc                                                                                        <  1:465094/98‑1 (MQ=16)
                  cgcgcAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCgg                                                                                    >  2:1363600/1‑101 (MQ=255)
                   gcgcAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGa                                                                                   >  1:2344980/1‑101 (MQ=255)
                    cgcAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGaa                                                                                  >  2:2248411/1‑101 (MQ=255)
                    cgcAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCACCGCTTCTGCCGCGGaa                                                                                  >  2:4244487/1‑101 (MQ=255)
                      cAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAAt                                                                                >  1:3545823/1‑101 (MQ=255)
                          ggCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTgg                                                                            >  2:2364894/1‑101 (MQ=255)
                           gCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGt                                                                           >  2:2051022/1‑101 (MQ=255)
                              tATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGa                                                                        >  1:2293479/1‑101 (MQ=255)
                              tATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGa                                                                        >  2:4441939/1‑101 (MQ=255)
                               aTTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGaa                                                                       >  2:2820404/1‑101 (MQ=255)
                                tttGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGaaa                                                                      >  1:2995348/1‑101 (MQ=255)
                                 ttGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAg                                                                     >  2:534461/1‑101 (MQ=255)
                                  tGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAgg                                                                    >  2:749556/1‑101 (MQ=255)
                                  tGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAgg                                                                    >  1:1718424/1‑101 (MQ=255)
                                  tGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAgg                                                                    >  1:1320600/1‑101 (MQ=255)
                                   ggACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAggg                                                                   >  1:3124103/1‑101 (MQ=255)
                                   ggACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAggg                                                                   >  1:1556623/1‑101 (MQ=255)
                                    gACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGt                                                                  >  1:2979629/1‑101 (MQ=255)
                                     aCCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTc                                                                 >  2:3393794/1‑101 (MQ=255)
                                     aCCAGGAAAAAGCTGCGTGTTTTTTCTGGGTGAAAATA..CGCCAGGGGGGTGCGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTc                                                                 <  1:1376436/101‑1 (MQ=255)
                                       cAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCgg                                                               >  1:4099535/1‑101 (MQ=255)
                                       aaGGAAAAAGTGGGTTTTTTTTTATGGGAAAAAAGA..CGCCGGGGGGGGGCGGGGGGGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCgg                                                               <  1:1868429/100‑1 (MQ=255)
                                        aGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCggg                                                              >  2:2510090/1‑101 (MQ=255)
                                         ggAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCgggg                                                             >  2:3226615/1‑101 (MQ=255)
                                           aaaaaGGTGCGTGTTTTTTCTTGGTGAAAATA..CGCCGGGGGGTTGGGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGtt                                                           <  1:1455453/101‑1 (MQ=37)
                                            aaaaGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGttt                                                          >  1:1771610/1‑101 (MQ=255)
                                             aaaGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGtttt                                                         >  1:382563/1‑101 (MQ=255)
                                              aaacgttgttttttttTTTTGGGAAAAAA..CGCGCCGGGGGGGGGGGGGGGGGAAACCCCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGttt                                                          <  1:149212/91‑1 (MQ=255)
                                               atgttggttttttttTTTTGGGAAAAAA..CGGGCGGGGGGGGGGGGGGGGTGAAACCCCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGtttt                                                         <  2:4067135/92‑1 (MQ=255)
                                               aGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTcc                                                       >  1:3496111/1‑101 (MQ=255)
                                                gttGCGTTTTTTTTCTTGTGGAAAATA..CGCCAGGGGGGTGGGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTccc                                                      <  1:1927889/99‑1 (MQ=255)
                                                gctggttgtTTTTTCTTGGTGAAAACC..CGCCGGGGGGGGGGGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTccc                                                      <  1:2911593/95‑1 (MQ=255)
                                                 cTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCg                                                     >  2:3906365/1‑101 (MQ=255)
                                                 agttgtgtTTTTTTTTTGGGAAAAAAGCCCCCAGGGGGGGGGGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTcc                                                       <  1:444423/97‑1 (MQ=255)
                                                    cGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGcc                                                  >  1:434345/1‑101 (MQ=255)
                                                      tgtTTTTTCTTGTTGAAAACC..CGCGACGGGGGTGCGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGTAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCt                                                <  2:4577863/101‑1 (MQ=38)
                                                      tgtTTTTTCTTGGTGAAAATC..CGCCAGGGGGGTGCGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCt                                                <  2:2976543/101‑1 (MQ=255)
                                                       gtTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtt                                               >  2:3530491/1‑101 (MQ=255)
                                                       gtTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtt                                               >  2:848678/1‑101 (MQ=255)
                                                        ttttttTTTGGTAAAAATG..CGCCACGGGGGTGCGGGGGTTTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCttt                                              <  2:4524087/101‑1 (MQ=38)
                                                        ttttttCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCttt                                              >  1:2678759/1‑101 (MQ=255)
                                                         tttttCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtttt                                             >  1:1078508/1‑101 (MQ=255)
                                                         tttttCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtttt                                             >  2:790496/1‑101 (MQ=255)
                                                         tttttCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtttt                                             >  1:1202023/1‑101 (MQ=255)
                                                             ttttttggAAAAAA..CGCCCCGGGGGGGGGGGGGGTTGCAACCCCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTc                                           <  1:1078602/94‑1 (MQ=255)
                                                               ttcgtaaaaaTA..CGCCACGGGGGTGGGGGCGTGTAAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTg                                       <  2:773520/95‑1 (MQ=38)
                                                                tGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGCCGGGGTTTTCCCGGCCCTTTTTCTTTgg                                      >  2:487504/1‑101 (MQ=255)
                                                                   taaaaaTC..CGCCGGGGGGGTGGGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCg                                   <  2:2293317/99‑1 (MQ=255)
                                                                    gAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCgg                                  >  1:91121/1‑101 (MQ=255)
                                                                      aaaTA..CGCCACGGGGGTTCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGACgggg                                >  1:721412/1‑101 (MQ=255)
                                                                       aaTA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGTCCCAAAGGGCCGGGCTTTTCCGGGCCTTTTTTCTTTGG.CGGGGcc                              >  2:4575271/1‑101 (MQ=25)
                                                                       aaTA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGc                               <  1:2135860/101‑1 (MQ=255)
                                                                       aaTA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGc                               >  2:4633758/1‑101 (MQ=255)
                                                                         tA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCt                             >  2:1423183/1‑101 (MQ=255)
                                                                         tA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCt                             >  1:2458485/1‑101 (MQ=255)
                                                                         aaaacccccggggggggggggggggTGAAACCCCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCggg                                 <  1:230370/91‑1 (MQ=255)
                                                                          atgcGCCACGGGGGTCCNGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGcc                              <  2:2102234/98‑1 (MQ=255)
                                                                          a..CGCCACGGGGGTGGGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCtt                            <  2:1433183/101‑1 (MQ=255)
                                                                           tgcgcccggggggggggggggttTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCt                             <  1:638465/93‑1 (MQ=255)
                                                                           ggcggggggggggggggggggtgCCACC..CCCGCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCt                             <  1:3183183/93‑1 (MQ=255)
                                                                            gcGCCAGGGGGGTGCGGGGGTGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCtt                            <  2:809785/100‑1 (MQ=255)
                                                                              gcggggggggtggggggttgTGACC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgg                          <  2:2344818/96‑1 (MQ=255)
                                                                                caggggggtgcgggggtgTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggc                        <  2:2811775/98‑1 (MQ=255)
                                                                                cagcggggtgcgggggtgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggc                        <  2:3111590/97‑1 (MQ=255)
                                                                                 aCGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGGTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggcc                       >  2:104571/1‑101 (MQ=255)
                                                                                  ggggggtgcgggggtgGGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccg                      <  2:3636831/100‑1 (MQ=255)
                                                                                  cgggggtgcgggggtgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccg                      >  2:205629/1‑101 (MQ=255)
                                                                                   gggggtgcgggggtgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccgg                     >  2:1557650/1‑101 (MQ=255)
                                                                                   gggggtgcgggggtgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccgg                     >  1:1914719/1‑101 (MQ=255)
                                                                                         gcgggggtgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCcgcg               >  1:4195125/1‑101 (MQ=255)
                                                                                           gggggtgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGcgac             >  2:1912765/1‑101 (MQ=255)
                                                                                             gggtgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCTTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGcgacga           >  1:1982176/1‑101 (MQ=255)
                                                                                              ggtgTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgag          <  2:3076218/101‑1 (MQ=255)
                                                                                                tgtgAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgagga        >  1:3689129/1‑101 (MQ=255)
                                                                                                 gtgAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgaggag       >  2:4262917/1‑101 (MQ=255)
                                                                                                  tgAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGgaga      >  1:1264856/1‑101 (MQ=255)
                                                                                                   gAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAGAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGaa     >  1:4641763/1‑101 (MQ=255)
                                                                                                   gAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGaa     >  2:2941244/1‑101 (MQ=255)
                                                                                                   gAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGaa     >  2:4478137/1‑101 (MQ=255)
                                                                                                     aa..CCCCCGCTTCTGCCGCGGAAATCTGGTCGCAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACg   >  1:4176008/1‑101 (MQ=255)
                                                                                                     aa..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACg   >  2:4395052/1‑101 (MQ=255)
                                                                                                      a..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACgg  >  2:2231069/1‑101 (MQ=255)
                                                                                                       |                                                                                                     
AGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATA..CGCCACGGGGGTGCGGGGGTGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACGG  >  NC_014034/1820293‑1820493

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG