Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,214 Δ1 bp 100% intergenic (‑142/+83) 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,2140C→.100.0% 24.5 138intergenic (‑142/+83)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (67/61):  ref base (0/10):  total (67/71)
Fisher's exact test strand distribution p-value = 1.44e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.14e-07

Only 100 of 177 total aligned reads displayed.

ATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCG.G.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCTTGA  >  NC_014034/1696111‑1696313
                                                                                                         |                                                                                                   
aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                                                     <  2:891629/101‑1 (MQ=16)
aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                                                     >  1:1397737/1‑101 (MQ=255)
  ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                                                   >  1:4219327/1‑99 (MQ=12)
   tCGAATCCTTGCACTAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                                                  <  2:3180335/101‑4 (MQ=12)
   tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                                                  >  1:167125/1‑98 (MQ=255)
   tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                                                  >  1:1724768/1‑98 (MQ=255)
   tCGAATCCTTGCACGAATTTTACAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                                                  >  1:3439253/1‑98 (MQ=255)
     gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                                                <  2:3944851/101‑6 (MQ=255)
     gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                                                <  2:1704277/101‑6 (MQ=255)
     gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                                                <  1:301854/101‑6 (MQ=255)
      aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                                               <  1:67973/101‑7 (MQ=12)
        tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                                                             <  2:2846411/101‑9 (MQ=255)
         ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                                                            <  2:373607/101‑10 (MQ=9)
         ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                                                            <  1:939666/101‑10 (MQ=9)
         ccTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgctt                                                                                            <  1:4195133/101‑10 (MQ=9)
          cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                                                           >  1:3527130/1‑91 (MQ=255)
          cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                                                           <  1:3125239/101‑11 (MQ=255)
             gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                                                      >  2:3648841/1‑101 (MQ=255)
              cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                                                     >  1:2849499/1‑101 (MQ=9)
              cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                                                     <  1:1346254/101‑1 (MQ=255)
               aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCCATTGACCCCGGGATGTTTGGGCCAA..GGTGAAA.TACCCT.GCTTTTTCtt                                                                                    >  2:393602/1‑101 (MQ=255)
               aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                                    >  1:516253/1‑101 (MQ=255)
               aCGAATTGTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                                    <  2:2022773/101‑1 (MQ=2)
                cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                                   <  1:4094066/101‑1 (MQ=255)
                  aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCGAC.CAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                                                                 <  2:1667877/101‑1 (MQ=255)
                  aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                                                                 >  1:3861251/1‑101 (MQ=255)
                   aTTTTGCAAGGAGGGGAG.AGGGGAGC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCCA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                                                >  2:1797521/1‑101 (MQ=255)
                    ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                                               >  1:1889287/1‑101 (MQ=9)
                    ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTAGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                                               >  2:3575452/1‑101 (MQ=9)
                     tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCC.AAGGG.GAAA..AACCTTGCTTTTTCTTGGCGCa                                                                              >  1:3720776/1‑101 (MQ=255)
                      ttGCAAGGAGGGGATAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                                                             <  2:3871389/101‑1 (MQ=255)
                      ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCCATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                                                             >  1:4403641/1‑101 (MQ=255)
                       tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..NGNNNAA.CA.CCTTGNTTTTTCTTGGGGCaaa                                                                            >  1:4905/1‑101 (MQ=2)
                        gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaaa                                                                           >  2:3339231/1‑100 (MQ=255)
                        gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                                           <  2:3814493/101‑1 (MQ=9)
                        gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                                           >  2:4629246/1‑101 (MQ=255)
                          aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                         <  2:2614232/101‑1 (MQ=255)
                          aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                         >  2:1518390/1‑101 (MQ=9)
                          aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                         <  2:101686/101‑1 (MQ=255)
                          aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                         <  1:4078852/101‑1 (MQ=255)
                          aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                         >  2:4253141/1‑101 (MQ=9)
                            ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                                                       <  2:4429655/101‑1 (MQ=255)
                            ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                                                       <  1:1076429/101‑1 (MQ=255)
                             gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                                                      <  1:1213397/101‑1 (MQ=255)
                              aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                                                     <  1:4253429/101‑1 (MQ=9)
                               ggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctc                                                                    <  1:1548295/101‑1 (MQ=255)
                               ggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctc                                                                    >  1:2678334/1‑101 (MQ=12)
                                gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                                   <  1:4293058/101‑1 (MQ=255)
                                gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                                   <  2:4413375/101‑1 (MQ=255)
                                     agaggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                                                <  1:3326053/91‑1 (MQ=255)
                                      gaggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                                               <  2:3143671/92‑1 (MQ=9)
                                      gaggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                                               <  2:551620/92‑1 (MQ=9)
                                      gaggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                                               <  1:972764/92‑1 (MQ=12)
                                      gaggggcgctGCACCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                                               >  1:2519476/10‑101 (MQ=12)
                                       aggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                                              <  2:3044577/93‑1 (MQ=255)
                                        ggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAAC..TCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                                             <  2:4364600/94‑1 (MQ=255)
                                        ggggcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                                             <  1:1620916/94‑1 (MQ=255)
                                           gcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                                                          <  1:109706/97‑1 (MQ=255)
                                           gcgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                                                          <  1:2745883/97‑1 (MQ=255)
                                            tgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                                                         <  2:1442065/98‑1 (MQ=16)
                                            cgctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAgg                                                         <  2:3699720/98‑1 (MQ=16)
                                             gctGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAggg                                                        >  2:306874/3‑101 (MQ=16)
                                               tGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                                      >  2:4050151/1‑101 (MQ=19)
                                               tGCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                                      <  1:3943468/101‑1 (MQ=16)
                                                gCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                                                     <  2:2849283/101‑1 (MQ=255)
                                                gCCCCTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                                                     >  2:1205929/1‑101 (MQ=19)
                                                 ccccTCCTGGCCTTTCG.GCCAATTCTCCCCTGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                                    <  2:703134/101‑1 (MQ=16)
                                                 ccccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                                    >  1:2346823/1‑101 (MQ=16)
                                                   ccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                                                  <  1:1478686/101‑1 (MQ=255)
                                                   ccTCCTGGCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                                                  >  1:760081/1‑101 (MQ=255)
                                                          gCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa                                           >  1:3848836/1‑101 (MQ=255)
                                                          gCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa                                           <  1:2979462/101‑1 (MQ=19)
                                                          gCCTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa                                           >  2:2462267/1‑101 (MQ=255)
                                                            cTTTCG.GCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGTTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                                         <  1:1087991/101‑1 (MQ=19)
                                                            cTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGTTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                                         <  2:1369991/101‑1 (MQ=19)
                                                            cTTTCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                                         >  2:1419420/1‑101 (MQ=255)
                                                             tttCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg                                        >  1:3130591/1‑101 (MQ=255)
                                                             tttCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg                                        >  2:2228689/1‑101 (MQ=19)
                                                              ttCG.GCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgc                                       >  2:3247382/1‑101 (MQ=255)
                                                                tcggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcg                                      <  2:760006/96‑1 (MQ=19)
                                                                   gccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCgg                                   >  1:3763943/3‑101 (MQ=255)
                                                                      aaTTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc                                >  1:415995/1‑101 (MQ=25)
                                                                      aaTTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccc                                >  2:2488738/1‑101 (MQ=255)
                                                                        ttcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccccc                              >  2:4625903/6‑101 (MQ=255)
                                                                        ttCCCCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGccccc                              <  2:336457/101‑1 (MQ=25)
                                                                          ccccccGGGATTGTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCt                            <  1:838480/98‑1 (MQ=255)
                                                                            ccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGa                          <  1:1367221/100‑1 (MQ=37)
                                                                               ccgggatgTTTGGCC.AAGGGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATc                        <  2:1089553/94‑1 (MQ=255)
                                                                               cGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCg                       >  1:3162341/1‑101 (MQ=255)
                                                                                gggATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGAGGCCCCCCTGATCGt                      >  1:865735/1‑101 (MQ=255)
                                                                                  gatgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGc                    >  2:717685/6‑101 (MQ=255)
                                                                                    tgtttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCt                  <  1:4627096/98‑1 (MQ=37)
                                                                                     gtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTc                 >  1:649042/3‑101 (MQ=255)
                                                                                     gtttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATCCCCCTCTGCCCGGAAGGCCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTc                 <  2:3848551/99‑1 (MQ=255)
                                                                                      tttGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCa                <  1:3486613/100‑1 (MQ=37)
                                                                                      tttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCa                >  1:219409/2‑101 (MQ=255)
                                                                                            ggccaaggTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCg            <  2:2678149/95‑1 (MQ=38)
                                                                                                 aggTGAAAC..CCCTTGTTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAgcgc       <  2:1745006/100‑1 (MQ=255)
                                                                                                  ggTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACCCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGcc      <  1:4065566/101‑1 (MQ=255)
                                                                                                       aaCCA.CCTTGCTTTTTCTGGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCTTGa  <  2:3538558/101‑1 (MQ=255)
                                                                                                         |                                                                                                   
ATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTTCG.G.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCTTGA  >  NC_014034/1696111‑1696313

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG