Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,206 Δ1 bp 98.4% intergenic (‑134/+91) 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,2060G→.98.4% 2.2 ‑2.3 ‑1.8192intergenic (‑134/+91)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (87/102):  ref base (3/0):  total (90/102)
Fisher's exact test strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.53e-02

Only 100 of 194 total aligned reads displayed.

AATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAG  >  NC_014034/1696110‑1696300
                                                                                                 |                                                                                              
  tttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAaca                                                                                       <  2:4692694/101‑2 (MQ=255)
   ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                                      <  2:3067531/101‑3 (MQ=255)
    tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATACACCCCGGGACGTTTGGGCCAA..GGTGAAacacc                                                                                     >  1:2672944/1‑98 (MQ=255)
      gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGGCAA..GGTGAAacacctt                                                                                   >  2:1759935/1‑96 (MQ=12)
      gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                                   <  1:782048/101‑6 (MQ=255)
      gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaactt                                                                                   >  2:4582449/1‑99 (MQ=12)
       aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                                  >  2:4522496/1‑95 (MQ=255)
          ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                                               >  2:3613016/1‑92 (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)
              gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                                         >  2:1797634/1‑101 (MQ=255)
               cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                                        <  2:1157285/101‑1 (MQ=255)
                aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                       <  2:710272/101‑1 (MQ=9)
                aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                       <  1:1057439/101‑1 (MQ=9)
                 cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                      <  1:1365413/101‑1 (MQ=255)
                 cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                      >  2:3849871/1‑101 (MQ=255)
                 cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                      <  1:3842979/101‑1 (MQ=255)
                 cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                      <  2:4117171/101‑1 (MQ=255)
                 cGAATTTTGCAAGGAGGCGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                      <  2:4266530/101‑1 (MQ=255)
                  gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAAC..GCCTTGCTTTTTCTTgg                                                                     <  1:3614949/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)
                   aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                                                    <  2:3948760/101‑1 (MQ=9)
                    aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                                   <  2:1926297/101‑1 (MQ=255)
                     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)
                        tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                               >  1:4128518/1‑101 (MQ=9)
                        tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                               >  1:2061442/1‑101 (MQ=9)
                         gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                              >  1:4199260/1‑101 (MQ=255)
                         gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTG.TTTCA.CCTTGCTTTTTCTTGGCGCAAAt                                                              <  2:3028697/101‑1 (MQ=2)
                          cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                             <  1:2808289/101‑1 (MQ=255)
                           aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATAc                                                            <  1:4300157/101‑1 (MQ=255)
                           aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                            <  2:3106469/101‑1 (MQ=255)
                           aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                            >  1:712446/1‑101 (MQ=9)
                            aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                                           <  2:4439867/101‑1 (MQ=9)
                               aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                                        <  2:1293761/101‑1 (MQ=9)
                                 gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                      <  2:4372683/101‑1 (MQ=255)
                                 gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                      <  1:4361545/101‑1 (MQ=255)
                                 gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                      >  1:1118131/1‑101 (MQ=9)
                                       gaggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                                  >  2:4559569/10‑101 (MQ=16)
                                        aggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                                 >  1:3169384/9‑101 (MQ=16)
                                         ggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                                <  1:1123222/94‑1 (MQ=255)
                                         ggggcgctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCC.AAGGG.GAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                                >  2:2333317/8‑101 (MQ=12)
                                          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)
                                               ctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                                          <  1:3661463/100‑1 (MQ=255)
                                               ctGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                                          <  2:3821435/100‑1 (MQ=255)
                                                ttccgatcttGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                         <  1:3562819/92‑1 (MQ=16)
                                                tGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                         <  1:1342061/101‑1 (MQ=19)
                                                tGCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                         >  1:4184023/1‑101 (MQ=16)
                                                 gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                                        <  2:2964375/101‑1 (MQ=255)
                                                 gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                                        <  1:421007/101‑1 (MQ=255)
                                                 gCCCCTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCC.AAGGG.GAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                                        >  1:3388995/1‑101 (MQ=19)
                                                  ccccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                       <  1:2999592/101‑1 (MQ=255)
                                                    ccTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                                     >  1:171891/1‑101 (MQ=19)
                                                     cTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                                    <  1:4029361/101‑1 (MQ=19)
                                                     cTCCTGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                                    <  1:2242794/101‑1 (MQ=255)
                                                        ctGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCgcgc                                 <  2:142872/100‑1 (MQ=19)
                                                        ctGGCCTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCgcgc                                 >  2:2593146/2‑101 (MQ=255)
                                                          ggCCTTTCGGCCAATTCGCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                               <  1:2325558/101‑1 (MQ=19)
                                                          ggCCTTTCGGCCAATTCACCCCNGGANNNTNGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                               <  1:2708938/101‑1 (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:2139435/101‑1 (MQ=19)
                                                             cTTTCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                            >  2:2162204/1‑101 (MQ=255)
                                                              tttCGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg                           >  2:1241919/1‑101 (MQ=255)
                                                                tcggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcg                         >  2:2583068/6‑101 (MQ=255)
                                                                  ggccAATTCACCCCGGGATGTTTGGGCCAA..GGCGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg                       >  1:1643640/4‑101 (MQ=19)
                                                                     gcaattcgcccGGGGATGTTTGGGCAAA.GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGc                     <  2:1836203/93‑1 (MQ=255)
                                                                      aaTTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc                   >  1:649195/1‑101 (MQ=25)
                                                                      aaTTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc                   >  1:4413879/1‑101 (MQ=255)
                                                                      aaTTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc                   <  2:737674/101‑1 (MQ=25)
                                                                       attcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc                  >  1:3264214/7‑101 (MQ=255)
                                                                        ttcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CATTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGCAGGCCCCGccccc                 >  2:2468519/6‑101 (MQ=19)
                                                                          caccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCt               >  2:349396/4‑101 (MQ=255)
                                                                           cccccGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGTTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTg              <  2:1875194/99‑1 (MQ=25)
                                                                           cccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTg              <  2:1344501/99‑1 (MQ=37)
                                                                            ccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGa             >  2:270900/2‑101 (MQ=255)
                                                                            ccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGa             <  2:4670131/100‑1 (MQ=37)
                                                                             cccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGAt            >  1:1102114/1‑101 (MQ=255)
                                                                              ccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATc           >  1:880304/1‑101 (MQ=255)
                                                                                 ggatgtttGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTg        <  2:1102028/95‑1 (MQ=37)
                                                                                 ggatgtttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAACTACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGAGCGGCCCCCCTGATCGTg        >  2:3230932/7‑101 (MQ=255)
                                                                                  gggtgTTTGGGCAAA.GGTGAAACC.CCTTGG.TTTTTCTGGGCGCAAATGCGCTTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTg        <  1:4614699/97‑1 (MQ=255)
                                                                                  ggaGGTTTGGGCAAA.GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTg        <  1:1231192/99‑1 (MQ=255)
                                                                                  gatgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGc       >  2:2770782/6‑101 (MQ=255)
                                                                                   atgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGcc      >  2:4144068/5‑101 (MQ=255)
                                                                                   atgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGcc      >  2:3404605/5‑101 (MQ=255)
                                                                                    tgtttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCt     >  2:4222012/4‑101 (MQ=37)
                                                                                      tttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTGGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCa   <  2:2898267/100‑1 (MQ=255)
                                                                                       ttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAg  >  1:4398642/1‑101 (MQ=255)
                                                                                       ttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAg  >  2:169454/1‑101 (MQ=255)
                                                                                       ttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAg  >  1:4135601/1‑101 (MQ=255)
                                                                                       ttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAg  <  2:1161192/101‑1 (MQ=255)
                                                                                                 |                                                                                              
AATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAG  >  NC_014034/1696110‑1696300

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG