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% 23.2 159intergenic (‑142/+83)RCAP_rcc01564/mmsAisobutyryl‑CoA dehydrogenase/methylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  new base (76/73):  ref base (0/10):  total (76/83)
Fisher's exact test strand distribution p-value = 1.59e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.25e-08

Only 100 of 199 total aligned reads displayed.

ATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTT.CGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGC  >  NC_014034/1696111‑1696308
                                                                                                         |                                                                                              
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                                                                                               <  1:345317/101‑2 (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.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacct                                                                                            <  1:1251321/101‑5 (MQ=255)
     gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                                           <  2:3276096/101‑6 (MQ=255)
     gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                                           <  2:4523873/101‑6 (MQ=255)
     gAATCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctt                                                                                           >  2:3666589/1‑96 (MQ=12)
      aaTCCTTGCACGAATTTTGCAGGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                                          >  1:2864404/1‑95 (MQ=255)
      aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                                          >  1:3217323/1‑95 (MQ=255)
      aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttg                                                                                          <  2:5018585/101‑7 (MQ=12)
      aaTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacacctgg                                                                                          >  2:1390540/1‑95 (MQ=255)
       aTCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgc                                                                                         <  1:573757/101‑8 (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:1574110/101‑9 (MQ=255)
        tCCTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgct                                                                                        <  2:3637404/101‑9 (MQ=255)
          cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                                                      >  2:4579389/1‑91 (MQ=255)
          cTTGCACGAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAacaccttgcttt                                                                                      <  2:3552800/101‑11 (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)
               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)
                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)
                 gAATTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTgg                                                                             >  2:1956478/1‑101 (MQ=9)
                   aTTTTGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGgcg                                                                           <  2:2278915/101‑1 (MQ=255)
                     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)
                      ttGCAAGGAGGGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCaa                                                                        >  1:3350144/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.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATa                                                                     <  2:4151953/101‑1 (MQ=255)
                         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)
                           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)
                             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)
                              aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGGAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                                                >  1:4375659/1‑101 (MQ=255)
                              aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                                                >  2:2769706/1‑101 (MQ=12)
                              aggGGAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCct                                                                >  2:3835162/1‑101 (MQ=12)
                                gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                              <  2:2341899/101‑1 (MQ=255)
                                gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                              <  2:1681825/101‑1 (MQ=255)
                                gggAGAGGGGC.GC.TGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCctct                                                              >  1:1566168/1‑101 (MQ=9)
                                     agaggggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGc                                                           >  2:512860/11‑101 (MQ=9)
                                      gaggggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGcc                                                          >  2:3899256/10‑101 (MQ=16)
                                       aggggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                                         <  1:427612/93‑1 (MQ=255)
                                       aggggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGccc                                                         >  2:3501248/9‑101 (MQ=16)
                                        ggggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCg                                                        <  2:2121320/94‑1 (MQ=255)
                                          ggcgctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGaa                                                      >  2:496754/6‑101 (MQ=16)
                                              ctGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGc                                                  >  2:309115/2‑101 (MQ=16)
                                               tGCCCCTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCa                                                 <  1:3666867/101‑1 (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:3216431/101‑1 (MQ=255)
                                                 ccccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAg                                               >  2:4247008/1‑101 (MQ=19)
                                                   ccTCCTGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGcc                                             <  2:2373273/101‑1 (MQ=255)
                                                     tcctGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgc                                           >  1:370736/4‑101 (MQ=255)
                                                      cctGGCCTTT.CGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTTTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCcgcg                                          <  1:3513663/99‑1 (MQ=255)
                                                      cctGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGTCCCGAAGGGCAAgatcgg                                          >  1:4090752/3‑96 (MQ=255)
                                                        tGGCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGcc                                        <  2:1633382/101‑1 (MQ=19)
                                                         ggCCTTT.CGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                                       <  1:2164448/101‑1 (MQ=25)
                                                         ggCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCg                                       >  1:4154517/1‑101 (MQ=255)
                                                          gCCTTT.CGGCCAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGa                                      <  1:2332952/101‑1 (MQ=19)
                                                           ccTTT.CGGCCAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGaa                                     <  1:1898918/101‑1 (MQ=25)
                                                           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                                  >  1:3835587/1‑101 (MQ=255)
                                                                 cggccAATTCCCCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGgcgc                                <  1:2552609/97‑1 (MQ=25)
                                                                  ggccAATTCACCCCGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGcgcg                               >  2:692925/4‑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:2064566/4‑101 (MQ=255)
                                                                           cccccGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTg                      <  1:610035/99‑1 (MQ=25)
                                                                            ccccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGa                     >  1:4702833/2‑101 (MQ=255)
                                                                              ccGGGATGTTTGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATc                   <  2:3916789/101‑1 (MQ=37)
                                                                              ccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATc                   <  1:1139751/101‑1 (MQ=37)
                                                                              ccGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTAGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGACCCCCTGATc                   >  2:4759846/1‑101 (MQ=25)
                                                                               cGGGATGTTTGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCg                  >  2:2953442/1‑101 (MQ=255)
                                                                                  gatgtttGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGc               <  1:4620110/96‑1 (MQ=37)
                                                                                  gatgtttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGc               >  2:2944188/6‑101 (MQ=255)
                                                                                     gtttGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTc            <  1:2265746/99‑1 (MQ=37)
                                                                                     gtttGGGCCAA..GGTGAAA.CA.CCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTc            <  2:2546887/99‑1 (MQ=37)
                                                                                      tttGGGCCAA..GGTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCa           <  1:3025811/100‑1 (MQ=37)
                                                                                      tttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCa           >  1:2185681/2‑101 (MQ=255)
                                                                                       ttGGGCCAA..GGTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAg          >  1:2935178/1‑101 (MQ=255)
                                                                                           gggccaaggTGAAAC..CCCTTGCTTTTTCTTGGCGCAAATCCGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTc        <  1:439506/94‑1 (MQ=255)
                                                                                                 aggTGAAA.CACCTT.GCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAgcgc  >  2:3457461/2‑101 (MQ=255)
                                                                                                         |                                                                                              
ATTTCGAATCCTTGCACGAATTTTGCAAGGAGGGGAGAAGGGCAGCTTGCCCCTCATGGCCTTT.CGG.CAATTCCACCCGGGAAGGTTGGGCCAAAGGGTGAAACCACCCTGGCTTTTTCTTGGCGCAAATACGCCTCTGGCCCGAAGGGCAAGCCGCGCCGAAGGCGCGGCCCCCCTGATCGTGCCTCAGTCGAGCGC  >  NC_014034/1696111‑1696308

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