Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,849 C→G 14.9% intergenic (‑273/‑130) RCAP_rcc01601 ← / → tgt FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,732,8490C→G14.9% 8.5 ‑1.7 0.053intergenic (‑273/‑130)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (13/31):  total (14/39)
Fisher's exact test strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....GT...C...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCG  >  NC_014034/1732749‑1732946
                                                                                                    |                                                                                                                   
aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGc                                                                                                                     >  1:1639666/1‑101 (MQ=255)
 gCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCg                                                                                                                    >  2:908278/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgg                                                                                                                   >  2:225347/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgg                                                                                                                   >  1:2076891/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgg                                                                                                                   >  2:5793311/1‑101 (MQ=255)
   gggACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTGGGCAAAATACACCCTTCTTTCTTTTCGGGGGGGCggt                                                                                                                  >  1:5384219/1‑100 (MQ=255)
   gggACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggg                                                                                                                  >  2:2636154/1‑101 (MQ=255)
   gggACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggg                                                                                                                  >  1:3012101/1‑101 (MQ=255)
     gACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggg                                                                                                                <  1:514764/101‑1 (MQ=255)
     gACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggg                                                                                                                <  1:5296376/101‑1 (MQ=255)
     gACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggg                                                                                                                <  2:4310789/101‑1 (MQ=255)
      aCCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgggggg                                                                                                               >  1:4831564/1‑101 (MQ=255)
        cAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgggggggg                                                                                                             >  1:1427810/1‑99 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggg                                                                                                            >  1:3230074/1‑98 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggg                                                                                                            >  1:3487358/1‑98 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggg                                                                                                            >  2:4763931/1‑98 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGTGGTGtcg                                                                                                            <  2:1432655/101‑3 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGa                                                                                                            <  2:3533168/101‑1 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGaa                                                                                                           <  2:2660358/101‑1 (MQ=255)
            ccccccTGCAGCCGCATGCCGTCCGCGTCAACCGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGAGGc                                                                                                         >  1:2718732/1‑101 (MQ=255)
             cccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgggggggcgaggg                                                                                                        >  1:5521743/1‑94 (MQ=255)
              ccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGccc                                                                                                       >  1:1361725/1‑101 (MQ=255)
              ccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGccc                                                                                                       >  2:508623/1‑101 (MQ=255)
               cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGcccc                                                                                                      <  1:4138062/101‑1 (MQ=255)
               cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGcccc                                                                                                      <  1:2040799/101‑1 (MQ=255)
                ccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggggggcccc                                                                                                     >  1:3358738/1‑91 (MQ=255)
                ccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGgccc...c                                                                                                  >  1:1349763/1‑101 (MQ=255)
                ccTGCAGCCGCATCCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGcccc...c                                                                                                  <  1:4972910/101‑1 (MQ=255)
                 cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGcccc...c.c                                                                                                <  2:3586581/101‑1 (MQ=255)
                 cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGcgcc...c.c                                                                                                >  1:689170/1‑101 (MQ=255)
                  ggCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....g                                                                                           <  1:2636303/100‑1 (MQ=255)
                    cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....gcg                                                                                         <  2:2782920/101‑3 (MQ=255)
                    cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGGGcccc...c.cccg                                                                                             >  2:3538509/1‑98 (MQ=255)
                     aGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGgccc...c.cccgg                                                                                            >  2:4642426/1‑97 (MQ=255)
                     aGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGcccc...c.cccgg                                                                                            >  1:5503465/1‑97 (MQ=255)
                      gCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGgccc...c.ccccgg                                                                                           >  1:5588846/1‑96 (MQ=255)
                         gCATGCCGTCCGCGTCAACGGGCAAGGCCCCCTCTTTTCCTAGGCAAAATACCCCCTTTTTTCTTTGCGGGGGGGCGGGGGGTGAAGCCCC...C.C....gcgggtcg                                                                                    <  1:1370606/101‑8 (MQ=25)
                         gCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGgccc...c.cccggggcg                                                                                        >  1:3769685/1‑93 (MQ=255)
                          cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGTGGGGGGTGAAGCCCC...C.C....gcggggcgc                                                                                   >  1:4863383/1‑93 (MQ=255)
                                         aaCAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGa                                                             <  1:3845554/101‑1 (MQ=21)
                                            aGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGc                                                          <  2:991574/101‑1 (MQ=21)
                                                      ccggtttttCCTAGGAAAAATACGCCTTTCTTTCTTTTCGGGGGGGCGGGGGGTGAACCCCC...C.C....GC...GGGTCG....CCTCGGGCTTGCCCGAGGCGAAAGCTTcc                                                <  2:70871/96‑1 (MQ=17)
                                                       ccGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....GCGGGT...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTccc                                               <  1:3610843/101‑1 (MQ=21)
                                                          cTTTTCCTAGGCAAAATCCGCCCTTCTTTCTTTGCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTcccccc                                            <  2:3336542/101‑1 (MQ=18)
                                                          cTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTcccccc                                            <  2:5781576/101‑1 (MQ=21)
                                                             tcccTGGGCAAAATCCGCCCTTCTTTCTTTGCGGGGGGGCGGGGGGTGAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                         <  2:3229874/99‑1 (MQ=17)
                                                               ccgaGGAAAAATTCGCCCTTTTTGCTTTGGGGGGGGGCGGGGGGGCAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaa                                       <  1:599072/98‑1 (MQ=2)
                                                               ccaaGGAAAAATCCGCCTTTTTTTCTTTTGGGGGGGGCGGGGGGGAAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaa                                       <  2:335494/98‑1 (MQ=9)
                                                                cTAGGCAAAATCCCCCCTTTTTTCTTTGGGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaaa                                      <  1:2496240/101‑1 (MQ=14)
                                                                  agtcAAAATCCTCCTTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAAc                                    <  2:2537918/98‑1 (MQ=17)
                                                                    gCAAAGTACCCCCTTTTTTTTTTTCGGGGGGGGGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCt                                  <  2:5385751/101‑1 (MQ=17)
                                                                    gCAAAATACGCCCTTCTTTCTTGTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCt                                  <  2:1652602/101‑1 (MQ=18)
                                                                     taaaaTCCGCCCTCGTTTCTTTTGGGGGGGGCGGGAGGTGAAGCCCCTGGC.G....GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCtt                                 <  2:2685915/100‑1 (MQ=12)
                                                                     cAAAATCCCCCTTTTTTTTTTTGGGGGGGGGGGGGGGGTGACCCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCtt                                 <  1:2407497/101‑1 (MQ=2)
                                                                     cAAAATCCCCCTTTCTTTTTTTGGGGGGGGGGGGGGGGTGAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCtt                                 <  2:5068540/101‑1 (MQ=9)
                                                                     aaaattccGCCCTTCTTTCTTTGGGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCtt                                 <  2:1954077/94‑1 (MQ=14)
                                                                     aaaaactcccccTTTTTTTTTTGGGGGGGGGGGGGGGGGGAACCCCC...C.CCGCGGG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCt                                  <  2:3011927/92‑1 (MQ=2)
                                                                        aaTACGCCCTTCTTTCTTTTGGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGcc                              <  2:1669152/101‑1 (MQ=18)
                                                                         aTACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCt                             <  2:868409/101‑1 (MQ=21)
                                                                            cGCTCTTCTTTCTTTGCGGGGGGGCGGGGGGGGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcg                          <  2:805220/101‑1 (MQ=17)
                                                                               ccTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcgcgc                       <  2:5079872/101‑1 (MQ=21)
                                                                                    tttctttTCGGGGGGGCGGGGGGGGGAGCCCC...C.CCCG.GG...G...GG....CCCCCGGCGTGCCCCAGGCGGAAACGTCCCCCCCCGGAAACCCTGCCCCCGCGCGCCGc                  >  2:4423967/1‑101 (MQ=255)
                                                                                       tttttGGGGGGGGGGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt               <  2:5312555/100‑1 (MQ=18)
                                                                                        ttttCGGGGGGGGGGGGGGTGAAACCCC...C.CCGCGGG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt               <  2:1349675/101‑1 (MQ=21)
                                                                                        ttttCGGGGGGGCGGGGGGGGAAGCCCC...C.C....CC...C...GGGGGCCCCCGGGCTTGCCCCGGGGGAAAGCTTCCCCCCCGGAAACCCTGCCCGCcgcgcgcgccttc               >  2:2409782/1‑94 (MQ=255)
                                                                                        tttgggggggggGGGGGGGTGACCCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTc              <  1:4642702/96‑1 (MQ=18)
                                                                                            ggggggggCGGGGGGTGACCCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCtt          <  2:1427717/100‑1 (MQ=18)
                                                                                             gggggggGGGGGGGGGAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCttt         <  2:5614901/101‑1 (MQ=18)
                                                                                              ggggggCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTc        <  2:5463259/101‑1 (MQ=21)
                                                                                              ggggggCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTc        <  1:3473759/101‑1 (MQ=21)
                                                                                              ggggggCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTc        <  2:3769442/101‑1 (MQ=21)
                                                                                               gggggGGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCg       <  2:4416929/101‑1 (MQ=21)
                                                                                                ggggCGGGGGGGGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCgg      <  1:3502339/101‑1 (MQ=21)
                                                                                                  ggCGGGGGGTCAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGcc    <  2:2413746/101‑1 (MQ=21)
                                                                                                    cGGGGGGTGAAGCCCG...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTTCCTGCGCGCGTCGCTTTCCCTTTTGGCCCg  >  1:4881740/1‑101 (MQ=18)
                                                                                                    cGGGGGGTGAAGCCCCT..CGC....GTG..T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCg  <  1:3548126/101‑1 (MQ=25)
                                                                                                    cGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCg  >  2:2458531/1‑101 (MQ=21)
                                                                                                    |                                                                                                                   
AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....GT...C...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCG  >  NC_014034/1732749‑1732946

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG