Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,185 A→C 100% intergenic (‑113/+112) RCAP_rcc01564 ← / ← mmsA isobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,696,1850A→C100.0% 10.4 231intergenic (‑113/+112)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (99/123):  ref base (7/2):  total (106/125)
Fisher's exact test strand distribution p-value = 8.39e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.89e-03

Only 100 of 237 total aligned reads displayed.

AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCC  >  NC_014034/1696085‑1696283
                                                                                                     |                                                                                                  
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                                <  1:1126888/101‑2 (MQ=37)
   tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAccccg                                                                                               >  2:4312960/1‑99 (MQ=255)
    gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGGGGGGC.GC.TGCCCCTCCTGGGCCTCCGG.CAAAACCACCCagg                                                                                              >  2:3255501/1‑101 (MQ=255)
    gCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccccg                                                                                              >  2:1341969/1‑98 (MQ=37)
      cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCCCCTCGCCTTTCGGCCCATTCACCCCGGGa                                                                                             >  1:186739/1‑101 (MQ=25)
       tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGat                                                                                            <  2:651641/101‑2 (MQ=255)
         aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgt                                                                                          >  1:3327866/1‑98 (MQ=255)
         aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgt                                                                                          >  1:3727796/1‑98 (MQ=255)
         aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCACCCCGGGAGGt                                                                                         >  1:1269887/1‑101 (MQ=255)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                         <  1:3467403/101‑5 (MQ=255)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                         <  1:2224912/101‑5 (MQ=255)
            tGAACGTCTGTTTAATTTCGAATTCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTg                                                                                       >  1:3821670/1‑101 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                                      >  2:1923292/1‑101 (MQ=25)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                                      <  2:3179951/101‑1 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                                      <  2:4478891/101‑1 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                                      >  1:4372930/1‑101 (MQ=25)
                  tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaa                                                                                 <  1:1605482/101‑1 (MQ=255)
                   cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCCCCCCCGGTTGTTTGGGCCatg                                                                                >  2:4744232/1‑99 (MQ=19)
                    tGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaagg                                                                               <  2:1584926/101‑3 (MQ=19)
                    tGTTTAATTTCGAATCCTTGCACGAATTTTACAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaagg                                                                               >  2:3885032/1‑99 (MQ=255)
                     gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggt                                                                              <  1:928700/101‑4 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                             <  1:4212243/101‑5 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                             <  1:2777073/101‑5 (MQ=255)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                            <  1:554618/101‑6 (MQ=16)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                            >  1:3910152/1‑96 (MQ=255)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                            <  2:3316550/101‑6 (MQ=16)
                        tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                           <  1:2939372/101‑7 (MQ=255)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCCATTCACCCCGGGATGTTTGGGCCAAGgggaaa                                                                          >  1:4033543/1‑98 (MQ=255)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                          <  2:672659/101‑8 (MQ=255)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                          <  2:1171172/101‑8 (MQ=255)
                          aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaagc                                                                         <  1:2498042/101‑9 (MQ=16)
                          aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                       >  2:1888661/1‑101 (MQ=255)
                            ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                     <  2:3067531/101‑3 (MQ=255)
                             tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATACACCCCGGGACGTTTGGGCCAA..GGTGAAacacc                                                                    >  1:2672944/1‑98 (MQ=255)
                              cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTTGGTCAA..GGTGAAacacct                                                                   >  2:3486337/1‑97 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                  <  2:3841958/101‑6 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                  <  1:782048/101‑6 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                  <  1:4313207/101‑6 (MQ=255)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                 <  2:468101/101‑7 (MQ=12)
                                  tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                               <  1:3640741/101‑9 (MQ=255)
                                   ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                              >  2:3613016/1‑92 (MQ=255)
                                   ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                              <  2:2765040/101‑10 (MQ=12)
                                       gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                        <  1:2703290/101‑1 (MQ=255)
                                       gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                        <  2:2563312/101‑1 (MQ=255)
                                         aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                      <  1:1057439/101‑1 (MQ=9)
                                          cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                     <  1:653614/101‑1 (MQ=255)
                                          cGAATTTTGCAAGGAGGGGAG.AGGGTCGC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                     <  2:1153863/101‑1 (MQ=255)
                                          cGAATTTTGCAAGGAGGCGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                     <  2:4266530/101‑1 (MQ=255)
                                              ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                 >  2:2030324/1‑101 (MQ=9)
                                               tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                                >  1:3023257/1‑101 (MQ=255)
                                               tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                                >  1:2371884/1‑101 (MQ=255)
                                               tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                                >  1:4195384/1‑101 (MQ=255)
                                                ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                               <  1:1356035/101‑1 (MQ=255)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              >  1:4128518/1‑101 (MQ=9)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              <  1:2031155/101‑1 (MQ=255)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                              >  1:4670387/1‑101 (MQ=9)
                                                  gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                             <  2:1727680/101‑1 (MQ=9)
                                                   cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            <  1:1334362/101‑1 (MQ=255)
                                                   cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                            <  1:2808289/101‑1 (MQ=255)
                                                    aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           >  2:4569099/1‑101 (MQ=9)
                                                    aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           >  1:1508644/1‑101 (MQ=9)
                                                    aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                           >  1:712446/1‑101 (MQ=9)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                        >  2:4614446/1‑101 (MQ=9)
                                                        aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                       <  2:1293761/101‑1 (MQ=9)
                                                          gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                     <  2:575617/101‑1 (MQ=255)
                                                               agaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGc                                  <  1:4463129/91‑1 (MQ=255)
                                                                gaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                 <  2:1813504/92‑1 (MQ=9)
                                                                gaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                 <  2:835501/92‑1 (MQ=9)
                                                                gaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                 >  1:4235738/10‑101 (MQ=16)
                                                                 aggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                <  1:1119077/93‑1 (MQ=255)
                                                                 aggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                >  2:179491/9‑101 (MQ=16)
                                                                 aggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                <  2:441246/93‑1 (MQ=255)
                                                                 aggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                >  1:3169384/9‑101 (MQ=16)
                                                                 aggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                >  2:4129428/9‑101 (MQ=16)
                                                                 aggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                <  2:2097078/93‑1 (MQ=255)
                                                                   gggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGa                              >  1:1666108/7‑101 (MQ=16)
                                                                    ggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                             >  2:4194955/6‑101 (MQ=16)
                                                                     gcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                            <  1:1582882/97‑1 (MQ=255)
                                                                      cgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                           >  2:3676459/4‑101 (MQ=16)
                                                                        ctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                         <  1:3968355/100‑1 (MQ=19)
                                                                        ctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                         >  1:2624254/2‑101 (MQ=16)
                                                                         ttccgatcttGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                        <  1:3562819/92‑1 (MQ=16)
                                                                          gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                       <  2:2964375/101‑1 (MQ=255)
                                                                           ccccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                      >  2:2877813/1‑101 (MQ=16)
                                                                             ccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                    >  1:171891/1‑101 (MQ=19)
                                                                              cTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                   <  1:2242794/101‑1 (MQ=255)
                                                                                 ctGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCgcgc                >  1:575723/2‑101 (MQ=19)
                                                                                 ctGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCgcgc                >  2:2593146/2‑101 (MQ=255)
                                                                                   ggCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg              >  2:3005208/1‑101 (MQ=255)
                                                                                     ccTTTCGGCCAATTGACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa            <  2:1236462/101‑1 (MQ=19)
                                                                                     ccTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa            >  2:2395017/1‑101 (MQ=255)
                                                                                     ccTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa            <  1:1923404/101‑1 (MQ=19)
                                                                                     ccTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa            <  1:179516/101‑1 (MQ=19)
                                                                                      cTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg           >  2:2162204/1‑101 (MQ=255)
                                                                                         tcggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcg        >  2:2299640/6‑101 (MQ=255)
                                                                                         tcggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcg        >  2:2583068/6‑101 (MQ=255)
                                                                                         tcggccAATTCACCCCGGGATGTTTGGGCC..AGGGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcg        >  1:4205818/6‑101 (MQ=255)
                                                                                              cAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcc   >  2:4542153/1‑101 (MQ=255)
                                                                                               aaTTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATCCGCCTCGGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc  >  1:3862916/1‑101 (MQ=255)
                                                                                               aaTTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc  <  2:737674/101‑1 (MQ=25)
                                                                                                     |                                                                                                  
AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCC  >  NC_014034/1696085‑1696283

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG