Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,696,217 Δ1 bp 100% intergenic (‑145/+80) 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,2170C→.100.0% 77.0 159intergenic (‑145/+80)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (68/68):  ref base (8/15):  total (76/83)
Fisher's exact test strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.55e-04

Only 100 of 167 total aligned reads displayed.

GCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTT.CGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCTTG  >  NC_014034/1696124‑1696312
                                                                                               |                                                                                               
gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                                                     <  1:4274413/101‑1 (MQ=255)
gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                                                     <  1:1407018/101‑1 (MQ=255)
gCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTc                                                                                     <  2:3239260/101‑1 (MQ=255)
 cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                                                    <  1:4762993/101‑1 (MQ=255)
 cACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCt                                                                                    >  2:2475489/1‑101 (MQ=9)
  aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                                   <  1:4993133/101‑1 (MQ=9)
  aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                                   >  2:1713609/1‑101 (MQ=255)
  aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                                   <  2:3719521/101‑1 (MQ=9)
  aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                                   <  2:33411/101‑1 (MQ=9)
  aCGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCtt                                                                                   >  2:2193312/1‑101 (MQ=255)
   cGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTg                                                                                  <  2:595561/101‑1 (MQ=255)
    gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGGTGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                                                 >  1:3797807/1‑101 (MQ=9)
    gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                                                 <  1:3926370/101‑1 (MQ=255)
    gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                                                 >  2:1956478/1‑101 (MQ=9)
     aaTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgc                                                                                <  2:3878203/101‑1 (MQ=9)
      aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                                               <  2:2278915/101‑1 (MQ=255)
      aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                                               <  1:2295652/101‑1 (MQ=255)
       ttttGCCAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                                              <  2:1449166/101‑1 (MQ=255)
       ttttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcgc                                                                              >  1:3214769/1‑101 (MQ=9)
        tttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.TGGCCAATTCACCCCGGGATGGTTGGGCC.AAGGG.GAAA.CA.CCTTGCTTTTTCTTGGCGCa                                                                             >  1:2121468/1‑101 (MQ=255)
         ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                                                            >  2:4594879/1‑101 (MQ=255)
          tGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaaa                                                                           >  2:1874234/1‑101 (MQ=9)
           gCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAAt                                                                          <  1:492109/101‑1 (MQ=9)
            cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCCATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                                         >  2:2731217/1‑101 (MQ=255)
            cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                                         >  2:4529348/1‑101 (MQ=255)
            cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                                         <  2:4151953/101‑1 (MQ=255)
            cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                                         <  2:5037115/101‑1 (MQ=255)
            cAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                                         >  2:355748/1‑101 (MQ=255)
             aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                        >  2:2800421/1‑101 (MQ=9)
             aaGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATAc                                                                        >  2:2371384/1‑101 (MQ=9)
              aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                                                       <  2:4568490/101‑1 (MQ=9)
              aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                                                       <  1:2371661/101‑1 (MQ=9)
              aggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACg                                                                       <  2:3797522/101‑1 (MQ=9)
               ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTTCCGC.CAATTCACCCCGGGATGTTTGGGCC.AAGGGGGAA..CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                                                      >  1:2698196/1‑101 (MQ=255)
               ggaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGc                                                                      >  2:4879442/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:1872402/101‑1 (MQ=255)
                gaggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGcc                                                                     <  1:1932128/101‑1 (MQ=255)
                 aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGGAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                                                    >  1:4375659/1‑101 (MQ=255)
                  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)
                   gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                                  <  2:1681825/101‑1 (MQ=255)
                          aggggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                                             <  1:427612/93‑1 (MQ=255)
                           ggggccctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                                            <  1:3455533/94‑1 (MQ=255)
                            gggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGa                                                           <  2:1238408/95‑1 (MQ=12)
                              gcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAg                                                         <  1:2184985/97‑1 (MQ=255)
                                 ctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                                                      >  1:4032118/2‑101 (MQ=16)
                                  tGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                                     <  2:2411898/101‑1 (MQ=16)
                                  tGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                                     <  1:3666867/101‑1 (MQ=16)
                                   gCCCCTCCTGGCCTTT.CGGCCGATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                                                    >  1:2871045/1‑101 (MQ=16)
                                   gCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCaa                                                    <  1:992890/101‑1 (MQ=255)
                                    ccccTCCTGGGCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                                   <  2:2804917/101‑1 (MQ=255)
                                    ccccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                                   <  2:1576724/101‑1 (MQ=255)
                                    ccccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                                   >  2:4247008/1‑101 (MQ=19)
                                     cccgcctGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGc                                                  <  2:2946894/95‑1 (MQ=16)
                                      ccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                                                 <  2:3210640/101‑1 (MQ=255)
                                       cTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcg                                                >  2:3907141/1‑101 (MQ=255)
                                        tcctGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgc                                               >  1:370736/4‑101 (MQ=255)
                                        tcctGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgc                                               <  2:2073793/98‑1 (MQ=19)
                                           tGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc                                            <  2:1633382/101‑1 (MQ=19)
                                            ggCCTTT.CGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                                           <  1:34164/101‑1 (MQ=25)
                                            ggCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                                           >  2:3534139/1‑101 (MQ=255)
                                              ccTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa                                         >  2:4836454/1‑101 (MQ=255)
                                              ccTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa                                         >  1:1110783/1‑101 (MQ=255)
                                               cTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAg                                        >  1:3159239/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)
                                                    cggccAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc                                    <  1:2552609/97‑1 (MQ=25)
                                                    cggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc                                    <  1:4716805/97‑1 (MQ=25)
                                                     ggccAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg                                   <  1:2362892/98‑1 (MQ=25)
                                                     ggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg                                   >  2:692925/4‑101 (MQ=255)
                                                      gccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCgg                                  >  2:3045488/3‑101 (MQ=255)
                                                          attcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccc                              <  2:735126/95‑1 (MQ=255)
                                                            tcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGcccccc                            >  1:4717429/5‑101 (MQ=255)
                                                            tcaccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGAcccccc                            >  1:2767977/5‑101 (MQ=25)
                                                             caccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCt                           >  2:3173522/4‑101 (MQ=255)
                                                             caccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCt                           >  2:3310289/4‑101 (MQ=255)
                                                              cccccGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTg                          <  1:610035/99‑1 (MQ=25)
                                                              accccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTg                          >  2:382133/3‑101 (MQ=255)
                                                                 ccGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATc                       <  2:3916789/101‑1 (MQ=37)
                                                                 ccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATc                       >  1:3839864/1‑101 (MQ=255)
                                                                 ccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATc                       <  1:1139751/101‑1 (MQ=37)
                                                                 ccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTAGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGACCCCCTGATc                       >  2:4759846/1‑101 (MQ=25)
                                                                  cGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCg                      >  1:1602145/1‑101 (MQ=255)
                                                                     gatgtttGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGc                   <  1:4620110/96‑1 (MQ=37)
                                                                     gatgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGc                   >  2:2944188/6‑101 (MQ=255)
                                                                       tgtttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCt                 >  1:4136189/4‑101 (MQ=255)
                                                                        gtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTc                >  2:3484727/3‑101 (MQ=255)
                                                                        gtttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTc                <  2:1996281/99‑1 (MQ=37)
                                                                         tttGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCa               <  1:3025811/100‑1 (MQ=37)
                                                                         tttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCa               >  1:2185681/2‑101 (MQ=255)
                                                                          ttGGGCCAA..GGTGAAAC..CCCTTGTTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAg              <  1:1900184/101‑1 (MQ=255)
                                                                          ttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAg              >  1:2935178/1‑101 (MQ=255)
                                                                             tgggccaaggTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGt             <  1:4443689/93‑1 (MQ=38)
                                                                             tgggccaaggTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGt             >  1:602219/9‑101 (MQ=255)
                                                                              gggccaaggTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTc            <  1:439506/94‑1 (MQ=255)
                                                                                 ccaaggTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAg         <  2:2775164/97‑1 (MQ=255)
                                                                                         tgaacccccTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCtt   <  2:4598118/95‑1 (MQ=255)
                                                                                         tgaaaccccTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCtt   <  1:3146183/95‑1 (MQ=255)
                                                                                          gaaaccccTTGCTTTTTCTTGGCCAAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCTTg  <  2:1744138/96‑1 (MQ=255)
                                                                                               |                                                                                               
GCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTT.CGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGCCTTG  >  NC_014034/1696124‑1696312

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG