Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,177 +C 95.1% intergenic (‑105/+120) 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,1771.→C95.1% 25.2 ‑4.4 ‑1.5223intergenic (‑105/+120)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (89/123):  ref base (9/2):  total (98/125)
Fisher's exact test strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.22e-02

Only 100 of 242 total aligned reads displayed.

TCCTCCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGC  >  NC_014034/1696077‑1696276
                                                                                                     |                                                                                                   
tcctccGAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.cc                                                                                                   <  1:759101/101‑2 (MQ=255)
tcctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.cc                                                                                                   <  2:627595/101‑2 (MQ=255)
tcctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.cc                                                                                                   <  2:2327514/101‑2 (MQ=255)
tcctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.cc                                                                                                   <  1:3873887/101‑2 (MQ=255)
 cctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.cca                                                                                                  >  1:368125/1‑99 (MQ=255)
  ctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.ccaa                                                                                                 >  2:2971665/1‑98 (MQ=255)
  ctccCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.ccaa                                                                                                 <  1:1775401/101‑4 (MQ=255)
    cccAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.ccaatt                                                                                               >  2:3492291/1‑96 (MQ=255)
     ccAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTc                                                                                               >  1:1727785/1‑101 (MQ=255)
       aaaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTcac                                                                                             >  1:2955219/1‑99 (MQ=255)
          gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                         <  1:1126888/101‑2 (MQ=37)
          gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCAcccc                                                                                         >  2:2406485/1‑100 (MQ=255)
           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.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGa                                                                                      <  2:4284249/101‑1 (MQ=255)
              cTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCCCCTCGCCTTTCGGCCCATTCACCCCGGGa                                                                                      >  1:186739/1‑101 (MQ=25)
               tGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGat                                                                                     >  1:2765204/1‑100 (MQ=37)
                gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatg                                                                                    <  1:1519087/101‑3 (MQ=25)
                 aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgt                                                                                   >  1:3327866/1‑98 (MQ=255)
                 aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgt                                                                                   <  2:652248/101‑4 (MQ=255)
                 aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGG.CCAATTCACCCCGGGAGGt                                                                                  >  1:1269887/1‑101 (MQ=255)
                 aaTTGAACGTCTGTTTAATTTAGAATCCTTGAACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgt                                                                                   >  1:2699879/1‑98 (MQ=255)
                  aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                  <  1:2224912/101‑5 (MQ=255)
                  aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGatgtt                                                                                  <  1:3467403/101‑5 (MQ=255)
                    tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTg                                                                                <  1:710718/101‑1 (MQ=255)
                     gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                               >  1:4372930/1‑101 (MQ=25)
                     gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTgg                                                                               >  2:2681950/1‑101 (MQ=25)
                      aaCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTggg                                                                              >  1:3949002/1‑101 (MQ=255)
                         gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCa                                                                           >  2:3485688/1‑101 (MQ=255)
                          tCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaa                                                                          <  1:1605482/101‑1 (MQ=255)
                           cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCCCCCCCGGTTGTTTGGGCCatg                                                                         >  2:4744232/1‑99 (MQ=19)
                           cTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAG.AGGGGAGC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCAaag                                                                         >  2:4590206/1‑100 (MQ=255)
                            tGTTTAATTTCGAATCCTTGCACGAATTTTACAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaagg                                                                        >  2:3885032/1‑99 (MQ=255)
                              tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                      <  1:4212243/101‑5 (MQ=255)
                              tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                      >  2:1231102/1‑97 (MQ=19)
                               ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGTGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                     <  1:1434500/101‑6 (MQ=16)
                               ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                     <  1:554618/101‑6 (MQ=16)
                               ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                     <  2:3316550/101‑6 (MQ=16)
                               ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                     <  1:2432864/101‑6 (MQ=16)
                                 aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCCATTCACCCCGGGATGTTTGGGCCAAGgggaaa                                                                   >  1:4033543/1‑98 (MQ=255)
                                 aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                   >  2:278642/1‑94 (MQ=16)
                                  aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                >  2:1888661/1‑101 (MQ=255)
                                    ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                              <  2:3067531/101‑3 (MQ=255)
                                       gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                           <  1:4313207/101‑6 (MQ=255)
                                        aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                          <  1:3987760/101‑7 (MQ=12)
                                        aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                          <  2:468101/101‑7 (MQ=12)
                                          tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                        <  1:3640741/101‑9 (MQ=255)
                                           ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                       <  1:4234604/101‑10 (MQ=9)
                                           ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                       >  1:2591201/1‑92 (MQ=255)
                                               gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                 <  2:2563312/101‑1 (MQ=255)
                                                cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                <  2:1157285/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)
                                                  cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                              <  2:4117171/101‑1 (MQ=255)
                                                  cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                              <  2:4507805/101‑1 (MQ=255)
                                                  cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                              <  1:3842979/101‑1 (MQ=255)
                                                   gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                             <  2:509463/101‑1 (MQ=255)
                                                    aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                            >  1:4688075/1‑101 (MQ=255)
                                                     aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                           <  2:1926297/101‑1 (MQ=255)
                                                      ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCCATTCACCCCGGGGTGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTTTTGgcgc                                          >  1:877508/1‑101 (MQ=255)
                                                      ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                          >  1:2926806/1‑101 (MQ=9)
                                                      ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                          >  2:2030324/1‑101 (MQ=9)
                                                       tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                         >  1:2371884/1‑101 (MQ=255)
                                                       tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCa                                         >  1:1940386/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)
                                                        ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                        <  2:2955043/101‑1 (MQ=255)
                                                          gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                      <  2:1727680/101‑1 (MQ=9)
                                                           cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                     <  1:1334362/101‑1 (MQ=255)
                                                            aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                    >  2:4569099/1‑101 (MQ=9)
                                                            aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                    <  2:3106469/101‑1 (MQ=255)
                                                             aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                   <  1:4175932/101‑1 (MQ=9)
                                                             aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                   <  2:2591049/101‑1 (MQ=9)
                                                             aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                   <  2:4439867/101‑1 (MQ=9)
                                                                aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                <  1:3767219/101‑1 (MQ=9)
                                                                aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                >  2:1253216/1‑101 (MQ=12)
                                                                       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)
                                                                         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:2097078/93‑1 (MQ=255)
                                                                         aggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                         >  1:3169384/9‑101 (MQ=16)
                                                                          ggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                        <  2:718598/94‑1 (MQ=255)
                                                                           gggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGa                       >  1:1666108/7‑101 (MQ=16)
                                                                            ggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                      <  1:2806646/96‑1 (MQ=16)
                                                                            ggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                      >  2:4194955/6‑101 (MQ=16)
                                                                                ctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                  <  2:3821435/100‑1 (MQ=255)
                                                                                 ttccgatcttGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                 <  1:3562819/92‑1 (MQ=16)
                                                                                   ccccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg               <  1:2999592/101‑1 (MQ=255)
                                                                                   ccccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg               <  2:362688/101‑1 (MQ=19)
                                                                                     ccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc             <  2:368081/101‑1 (MQ=255)
                                                                                      cTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg            <  1:2242794/101‑1 (MQ=255)
                                                                                         ctGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCgcgc         <  2:142872/100‑1 (MQ=19)
                                                                                           ggCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg       >  2:3005208/1‑101 (MQ=255)
                                                                                             ccTTTCGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa     <  1:958628/101‑1 (MQ=25)
                                                                                             ccTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa     <  1:1923404/101‑1 (MQ=19)
                                                                                              cTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg    >  2:2162204/1‑101 (MQ=255)
                                                                                                ttCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgc  >  1:265520/1‑101 (MQ=255)
                                                                                                     |                                                                                                   
TCCTCCCAAAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGC  >  NC_014034/1696077‑1696276

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