Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,820,393 +C 2.0% 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.0% 2.3 ‑4.2 ‑0.1152intergenic (+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 (114/35):  total (114/38)
Fisher's exact test strand distribution p-value = 1.47e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.23e-01

Only 100 of 184 total aligned reads displayed.

AGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACGG  >  NC_014034/1820293‑1820493
                                                                                                      |                                                                                                     
aGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGaaa                                                                                                        >  2:1327494/1‑101 (MQ=255)
aGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGaaa                                                                                                        >  2:4303478/1‑101 (MQ=255)
aGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGaaa                                                                                                        >  2:411574/1‑101 (MQ=255)
 gttGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..c                                                                                                     >  1:1260981/3‑101 (MQ=255)
  gTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..cc                                                                                                    >  1:3551443/1‑101 (MQ=255)
   tGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..ccc                                                                                                   >  1:3359082/1‑101 (MQ=255)
      gcgGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCg                                                                                                >  2:630678/1‑101 (MQ=255)
      gcgGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCg                                                                                                >  2:3599923/1‑101 (MQ=255)
       cgGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGc                                                                                               >  2:360361/1‑101 (MQ=255)
        gggCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCt                                                                                              >  1:2125009/1‑101 (MQ=255)
        gggCTTGGAGCGCGAAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTTTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTAAAA..CCTTCGCt                                                                                              <  2:911143/101‑1 (MQ=25)
         ggCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAC..CCCCCGCtt                                                                                             <  2:2478841/101‑1 (MQ=255)
         ggCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCtt                                                                                             >  1:119923/1‑101 (MQ=255)
          gCTTTGAGCGCGCAAGGGCTTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTAAAAATACGCCAGGGGGGTGCGGGGG.TGTGAAC..CCCCCGCTTc                                                                                            <  2:3512566/101‑1 (MQ=25)
          gCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTc                                                                                            >  1:1699460/1‑101 (MQ=255)
            tttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTg                                                                                          >  1:3458268/1‑101 (MQ=255)
            tttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTg                                                                                          >  1:934770/1‑101 (MQ=255)
            tttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTg                                                                                          >  2:2417924/1‑101 (MQ=255)
            tttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTg                                                                                          >  2:3458263/1‑101 (MQ=255)
             ttGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGc                                                                                         >  2:2336166/1‑101 (MQ=255)
                 gcgcgcAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAC..CCCCCGCTTCTGCcgcg                                                                                     <  2:765484/101‑1 (MQ=255)
                   gcgcAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGa                                                                                   >  1:3897336/1‑101 (MQ=255)
                     gcAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGaaa                                                                                 >  2:2067941/1‑101 (MQ=255)
                     gcAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGaaa                                                                                 >  1:1006490/1‑101 (MQ=255)
                       aaGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATc                                                                               >  2:1876504/1‑101 (MQ=255)
                         gggCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTg                                                                             >  2:1666875/1‑101 (MQ=255)
                           gCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGt                                                                           >  2:598638/1‑101 (MQ=255)
                           gCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGt                                                                           >  2:1224387/1‑101 (MQ=255)
                            cGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTc                                                                          >  1:3919258/1‑101 (MQ=255)
                             gTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCg                                                                         >  2:3649147/1‑101 (MQ=255)
                             gTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCg                                                                         >  1:2848314/1‑101 (MQ=255)
                             gTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCg                                                                         >  1:963863/1‑101 (MQ=255)
                               aTTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGaa                                                                       >  1:3607004/1‑101 (MQ=255)
                               aTTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGaa                                                                       >  2:1518603/1‑101 (MQ=255)
                               aTTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGaa                                                                       >  1:669095/1‑101 (MQ=255)
                                 ttGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATCCGCCGCGGGGTGGGGGGGG.TGTAAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAg                                                                     <  2:2993713/101‑1 (MQ=37)
                                  tGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAgg                                                                    <  2:770746/101‑1 (MQ=255)
                                      ccAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCGGGGGGGTGCGGGGG.TGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCg                                                                <  2:1440871/101‑1 (MQ=255)
                                      ccAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCg                                                                >  2:1055522/1‑101 (MQ=255)
                                       cAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CGCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCgg                                                               >  1:3104456/1‑101 (MQ=255)
                                       cAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCgg                                                               >  1:211442/1‑101 (MQ=255)
                                        aGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCggg                                                              >  2:3407278/1‑101 (MQ=255)
                                          gAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGt                                                            >  2:21843/1‑101 (MQ=255)
                                           aaaaaGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGtt                                                           >  1:1812228/1‑101 (MQ=255)
                                           aaaaaGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTATGCCGCGGAAATCTGGTCGAAAGGGTCgggggt                                                           >  1:3802051/1‑99 (MQ=255)
                                            aaaaGCTGCGTGTTTTTTCTGGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGttt                                                          <  2:4413325/101‑1 (MQ=255)
                                             aaaGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGtttt                                                         >  1:4203781/1‑101 (MQ=255)
                                             aaaGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGtttt                                                         >  2:3437464/1‑101 (MQ=255)
                                             aaaGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGGGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGtttt                                                         >  1:627640/1‑101 (MQ=255)
                                              aaGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTc                                                        >  1:2263708/1‑101 (MQ=255)
                                              aaGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGGGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTc                                                        >  1:4388122/1‑101 (MQ=255)
                                               aGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTcc                                                       >  1:4245277/1‑101 (MQ=255)
                                               aGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTcc                                                       >  1:1517507/1‑101 (MQ=255)
                                               aGCGGCGTTTTTTTTCTGGGTGAAAATACCCCGCGGGGGTGCGGGGG.TTTAAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTcc                                                       <  2:2600597/101‑1 (MQ=255)
                                                gCTGGGTTTTTTTTTTTGGAAAAAATACGCCGGGGGGGGGGGGGGG.TGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGCCGGGGTTTTccc                                                      <  1:2416243/101‑1 (MQ=255)
                                                gCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTccc                                                      >  2:1892335/1‑101 (MQ=255)
                                                agtttgtgtTTTTTTTTTGGGAAAAAAGGCCCGGGGGGGGGGGGGG.GTTGAAACCCCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTc                                                        <  1:809904/96‑1 (MQ=255)
                                                 cTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCg                                                     >  2:1726855/1‑101 (MQ=255)
                                                  ggtgttttttttCTGGGAAAAAATACGCCGGGGGGGGGGGGGGG.TGTAAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCgg                                                    <  1:470489/95‑1 (MQ=255)
                                                   gCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGc                                                   >  2:3217983/1‑101 (MQ=255)
                                                    cGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGCCGGGGTTTTCCCGGcc                                                  >  2:1795203/1‑101 (MQ=255)
                                                         tttttCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtttt                                             >  1:2567794/1‑101 (MQ=255)
                                                         tttttCTTGCTAAAAATCCCCCCAGGGGGGGCGGGGG.TGAAACA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtttt                                             <  2:2424570/101‑1 (MQ=255)
                                                         tttttCTGGGTAAAAACCCGCGGGGGGGGCGGGGGGG.TGGGACC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCtttt                                             <  1:630746/101‑1 (MQ=255)
                                                           tttCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTc                                           <  1:3633224/101‑1 (MQ=255)
                                                             tCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCtt                                         >  2:14089/1‑101 (MQ=255)
                                                               ttGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTg                                       >  2:2144309/1‑101 (MQ=255)
                                                                tGGTGAAAATCCGCCAGGGGGGTGCGGGGG.TGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTgg                                      <  2:3685711/101‑1 (MQ=255)
                                                                tGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCACGGCCCTTTTTCTTTgg                                      >  1:4097103/1‑101 (MQ=255)
                                                                tGGTGAAAAACCCCCGCGGGGGGGGGGGGGTTGAAAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTg                                       <  2:1701039/101‑1 (MQ=255)
                                                                 ggagAAAAGCCGCCCGGGGGGTGCGGGGT.TGGGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGc                                     <  1:1327591/98‑1 (MQ=255)
                                                                  gtaaaaaTACGCCCCGGGGGTGCGGGGG.TGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGcc                                    <  2:2684176/98‑1 (MQ=255)
                                                                   tGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCg                                   >  1:1892716/1‑101 (MQ=255)
                                                                    gAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCgg                                  >  2:2260712/1‑101 (MQ=255)
                                                                     aaaaTACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCggg                                 >  1:4421475/1‑101 (MQ=255)
                                                                      aaaTACGCCACGGGGGTGTGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCgggg                                <  2:2960070/101‑1 (MQ=255)
                                                                      aaaTACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCgggg                                >  2:3108731/1‑101 (MQ=255)
                                                                      aaaTACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCgggg                                >  1:2075215/1‑101 (MQ=255)
                                                                       aaTACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGc                               >  1:2780776/1‑101 (MQ=255)
                                                                        accggcgggggggggggggggg.tgGAAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGcc                              <  1:1072371/92‑1 (MQ=255)
                                                                         cgcgcgggggggggggggggg.tgGAACC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCt                             <  2:2561527/93‑1 (MQ=255)
                                                                          aCGCCGGGGGGGGGGGGGGG.TGTGAAC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCtt                            <  1:4224776/101‑1 (MQ=255)
                                                                            gCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgg                          >  1:1012469/1‑101 (MQ=255)
                                                                            gCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgg                          >  2:4025605/1‑101 (MQ=255)
                                                                            cgcccggggctgcggggg.tgTTAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgg                          <  2:4270374/97‑1 (MQ=255)
                                                                             ccACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTggg                         >  2:3578226/1‑101 (MQ=255)
                                                                              cACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggc                        >  1:3635131/1‑101 (MQ=255)
                                                                                 gggggtgtggggg.tgTACAC..CCCCTGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccgg                     <  2:4387850/101‑1 (MQ=255)
                                                                                 gggggtgcggggg.tgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccgg                     >  2:1850708/1‑101 (MQ=255)
                                                                                   ggggggggggg.tgCGACC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccgggc                   <  2:669027/95‑1 (MQ=255)
                                                                                    ggtggggggg.tgCAACC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTgggccgggcc                  <  2:3358845/101‑1 (MQ=255)
                                                                                     gtgcggggg.tgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGggccgggccg                 >  1:1539178/1‑101 (MQ=255)
                                                                                       ggggggg.tgTAAAC..CCCCTGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCcgcg               <  2:2440088/99‑1 (MQ=255)
                                                                                         ggggg.tgTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGcgac             >  2:2455398/1‑101 (MQ=255)
                                                                                         ggggg.agTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGcgac             >  2:1987177/1‑101 (MQ=255)
                                                                                          gggg.tgTAACC..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGcgacg            <  2:445243/101‑1 (MQ=255)
                                                                                             g.tgtGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgagg         >  2:558827/1‑101 (MQ=255)
                                                                                                  gAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACgaggagca     >  2:3283319/1‑99 (MQ=255)
                                                                                                  gAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGaa     >  1:2191974/1‑101 (MQ=255)
                                                                                                     a..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACgg  >  2:3917793/1‑101 (MQ=255)
                                                                                                      |                                                                                                     
AGGTGCGCGGGCTTTGAGCGCGCAAGGGCGTATTTGGACCAGGAAAAAGCTGCGTGTTTTTTCTTGGTGAAAATACGCCACGGGGGTGCGGGGG.TGTGAAA..CCCCCGCTTCTGCCGCGGAAATCTGGTCGAAAGGGTCGGGGTTTTCCCGGCCCTTTTTCTTTGGCCGGGGCCTTGGGCCGGGCCGCGACGAGGAGAACGG  >  NC_014034/1820293‑1820493

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG