Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,184 2 bp→AC 100% intergenic (‑112/+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,1840C→A100.0% 16.7 247intergenic (‑112/+113)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (113/112):  ref base (7/15):  total (120/127)
Fisher's exact test strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.62e-05
*NC_0140341,696,1850A→C100.0% 14.8 248intergenic (‑113/+112)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (114/124):  ref base (6/4):  total (120/128)
Fisher's exact test strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.79e-04

Only 100 of 252 total aligned reads displayed.

AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTT.CGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGG  >  NC_014034/1696085‑1696280
                                                                                                     ||                                                                                               
aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGG.CCAATTCAcc                                                                                               >  1:1319899/1‑101 (MQ=255)
aaGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGG.CCAATTCAcc                                                                                               <  1:4135828/101‑1 (MQ=255)
 aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGG.CCAATTCAccc                                                                                              >  2:1972227/1‑101 (MQ=255)
 aGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGG.CCAATTCAccc                                                                                              >  1:2361277/1‑101 (MQ=255)
  gTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGG.CCAATTCAcccc                                                                                             <  2:4148993/101‑2 (MQ=37)
   tGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGG.CCAATTCAccccg                                                                                            <  2:2159604/101‑3 (MQ=255)
     ccTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCggg                                                                                           >  2:2247346/1‑101 (MQ=255)
     ccTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGG.CCACTTCACCCcggg                                                                                          >  2:4219513/1‑100 (MQ=255)
       tGAATTGAACGTGTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGat                                                                                         >  1:4598440/1‑100 (MQ=25)
        gAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGatg                                                                                        <  2:3407625/101‑3 (MQ=255)
         aaTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGatgt                                                                                       >  2:439441/1‑98 (MQ=255)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGatgtt                                                                                      >  1:4152270/1‑97 (MQ=25)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGatgtt                                                                                      >  1:1744268/1‑97 (MQ=25)
          aTTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGatgtt                                                                                      >  2:1347690/1‑97 (MQ=25)
            tGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTg                                                                                    >  2:3722601/1‑101 (MQ=255)
             gAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTgg                                                                                   >  2:427547/1‑101 (MQ=25)
              aaCGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTggg                                                                                  >  2:1409730/1‑101 (MQ=255)
                cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGcc                                                                                >  1:2279070/1‑101 (MQ=25)
                cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGcc                                                                                <  1:2167604/101‑1 (MQ=255)
                cGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGcc                                                                                <  2:3691060/101‑1 (MQ=255)
                 gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCa                                                                               <  1:2652348/101‑1 (MQ=25)
                 gTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCa                                                                               >  1:3152117/1‑101 (MQ=255)
                     gTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggt                                                                           <  1:4044466/101‑4 (MQ=255)
                     gTTTAATTTCGAATACTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGGCaaggt                                                                           >  1:729460/1‑98 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                          <  2:2598396/101‑5 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                          >  2:3025573/1‑97 (MQ=19)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                          <  1:1972573/101‑5 (MQ=255)
                      tttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtg                                                                          >  2:1725329/1‑97 (MQ=19)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                         >  1:523665/1‑96 (MQ=255)
                       ttAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGACGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtga                                                                         <  2:4998717/101‑6 (MQ=16)
                        tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                        <  2:1541355/101‑7 (MQ=255)
                        tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                        >  1:376747/1‑95 (MQ=255)
                        tAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaa                                                                        >  2:2659307/1‑95 (MQ=255)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCaaggtgaaa                                                                       >  2:3209495/1‑94 (MQ=16)
                         aaTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAAGgggaaa                                                                       >  1:3690187/1‑98 (MQ=16)
                          aTTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAac                                                                    >  1:3210891/1‑101 (MQ=255)
                           tttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAaca                                                                   <  1:4219830/101‑2 (MQ=255)
                           tttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAaca                                                                   >  2:3228170/1‑100 (MQ=255)
                            ttCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacac                                                                  <  1:3846889/101‑3 (MQ=255)
                             tCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacc                                                                 >  2:2439141/1‑98 (MQ=255)
                              cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.TGGCCCATTCCCCCCCGGAGGTTTGGGCCAAGGGGAAAAacct                                                                  >  1:932033/1‑100 (MQ=255)
                              cGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACACCGGGATGTTTGGGCCAA..GGTGAAacacct                                                                >  1:323223/1‑97 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                               <  2:4523873/101‑6 (MQ=255)
                               gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                               <  1:2161499/101‑6 (MQ=255)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                              <  1:3565055/101‑7 (MQ=12)
                                aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                              <  2:5018585/101‑7 (MQ=12)
                                 aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                             >  2:2332788/1‑94 (MQ=255)
                                 aTCCTTGCACGAATTTTGCAAGGAGGGGAG.AGGGGAGC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                             >  2:1518565/1‑94 (MQ=255)
                                  tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                            <  2:1218947/101‑9 (MQ=255)
                                  tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                            <  2:3637404/101‑9 (MQ=255)
                                    cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                          <  2:3552800/101‑11 (MQ=255)
                                    cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                          >  2:4579389/1‑91 (MQ=255)
                                       gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                     <  1:1407018/101‑1 (MQ=255)
                                       gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                     <  1:4274413/101‑1 (MQ=255)
                                        cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                    <  1:4762993/101‑1 (MQ=255)
                                         aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                   <  1:4993133/101‑1 (MQ=9)
                                         aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                   <  2:33411/101‑1 (MQ=9)
                                          cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                  <  2:595561/101‑1 (MQ=255)
                                           gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                 <  1:3926370/101‑1 (MQ=255)
                                             aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                               <  2:2278915/101‑1 (MQ=255)
                                                ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                            >  1:3350144/1‑101 (MQ=255)
                                                 tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                           >  2:1033090/1‑101 (MQ=9)
                                                  gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                          <  2:1904539/101‑1 (MQ=9)
                                                   cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                         <  1:1637800/101‑1 (MQ=255)
                                                    aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                        >  2:2800421/1‑101 (MQ=9)
                                                    aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTATTTCTTGGCGCAAATAc                                        >  2:4545305/1‑101 (MQ=9)
                                                     aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                       <  1:2371661/101‑1 (MQ=9)
                                                      ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCCGC.CAATTCACCCCGGGATGTTTGGGCC.AAGGGGGAA..CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                      >  1:2698196/1‑101 (MQ=255)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGcc                                     <  1:2439315/101‑1 (MQ=255)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                     <  1:1932128/101‑1 (MQ=255)
                                                       gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                     <  1:1872402/101‑1 (MQ=255)
                                                        aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGGAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                    >  1:4375659/1‑101 (MQ=255)
                                                        aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                    >  2:3835162/1‑101 (MQ=12)
                                                         ggggAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctc                                   <  1:492532/101‑1 (MQ=255)
                                                          gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                  >  1:1566168/1‑101 (MQ=9)
                                                                gaggggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                              >  2:3899256/10‑101 (MQ=16)
                                                                  ggggccctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                            <  1:3455533/94‑1 (MQ=255)
                                                                    ggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                          >  2:4719998/6‑101 (MQ=16)
                                                                        ctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                      >  2:309115/2‑101 (MQ=16)
                                                                         tGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                     <  2:780163/101‑1 (MQ=16)
                                                                         ccgatcTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGTCCCGAAGGGCa                     <  2:4090449/96‑1 (MQ=16)
                                                                          gCCCCTCCTGGCCTTT.CGGCCGATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                    >  1:2871045/1‑101 (MQ=16)
                                                                           ccccTCCTGGGCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                   <  2:2804917/101‑1 (MQ=255)
                                                                           ccccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                   <  2:1688056/101‑1 (MQ=19)
                                                                           ccccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                   <  2:1576724/101‑1 (MQ=255)
                                                                           ccccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                   >  1:1192504/1‑101 (MQ=16)
                                                                             ccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                 <  2:2373273/101‑1 (MQ=255)
                                                                                 ctGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCgcgc             >  1:3054877/2‑101 (MQ=19)
                                                                                  tGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc            <  2:1633382/101‑1 (MQ=19)
                                                                                    gCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa          <  1:2332952/101‑1 (MQ=19)
                                                                                    gCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa          >  1:2804275/1‑101 (MQ=255)
                                                                                     ccTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa         >  1:1110783/1‑101 (MQ=255)
                                                                                     ccTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa         >  2:4836454/1‑101 (MQ=255)
                                                                                       ttt.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAgg       <  2:1266650/101‑1 (MQ=19)
                                                                                        tt.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgc      <  2:523586/101‑1 (MQ=19)
                                                                                        tt.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgc      >  1:3835587/1‑101 (MQ=255)
                                                                                           cggccAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc    <  1:2552609/97‑1 (MQ=25)
                                                                                           cggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc    <  1:4716805/97‑1 (MQ=25)
                                                                                            ggccAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg   <  2:3272317/98‑1 (MQ=25)
                                                                                             gccAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCgg  <  1:4739985/99‑1 (MQ=25)
                                                                                                     ||                                                                                               
AAGTGCCTGAATTGAACGTCTGTTTAATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTT.CGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGG  >  NC_014034/1696085‑1696280

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG