Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,078 C→G 3.2% T193T (ACG→ACC)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0780C→G 3.2% 23.5 ‑0.7 0.0522T193T (ACG→ACC) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (7/11):  ref base (243/261):  total (250/272)
Fisher's exact test strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 862 total aligned reads displayed.

GGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGC  >  NC_014034/221974‑222178
                                                                                                        |                                                                                                    
    cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGGCCGTGGGATCGgtcgggtcc                                                                                                      >  1:1125334/1‑101 (MQ=11)
    cggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGGTCGGTCGGGTCCGGGGGATCGgtcgggtcc                                                                                                      >  2:1082918/1‑101 (MQ=11)
     ggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccg                                                                                                     >  1:502489/1‑101 (MQ=11)
     ggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccg                                                                                                     >  1:1505642/1‑101 (MQ=11)
      gtcgggtccgtggggtccgtcgtATCCGTGGGACCGGTCGGATCTGTGGGATCGGTCGGATCCTTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                                                                                                    <  1:40657/101‑1 (MQ=11)
      gtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                                                                                                    >  1:3381230/1‑101 (MQ=11)
      gtcggggccgtggggtccggcggATCCGGGGAATCGGTGGGATCCGTGGGATGGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                                                                                                    <  1:4257244/101‑1 (MQ=11)
      gtcggatccgtgggttccgtcggACCCGTGGGATCGCTCGGATTCTGGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgt                                                                                                    <  2:582978/101‑1 (MQ=11)
       tcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtg                                                                                                   >  1:76083/1‑100 (MQ=11)
       tcgggtccgtcggatcggtgggATCCGTGGGATCGGTAGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccggg                                                                                                   >  2:3125710/1‑99 (MQ=11)
       tcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtc                                                                                                   >  1:4927141/1‑101 (MQ=255)
       tcggatccgtgggaccggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGTATCCGTGGGATCGGTCGGATCCGTAGGATCGgtcgggtccgta                                                                                                   <  2:4745859/101‑2 (MQ=11)
        cgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgg                                                                                                  <  1:3730649/101‑3 (MQ=11)
         ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGGTCAGTGGGATCGgcggggtccgtggg                                                                                                 >  2:3360319/4‑101 (MQ=11)
          ggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGGCCGTGGGATCGgtcgggtccgtggga                                                                                                >  2:4548153/1‑100 (MQ=11)
          gctcggtgggatcggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgga                                                                                                <  2:4042373/96‑2 (MQ=11)
          gatgcgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                                                                >  1:3498152/5‑101 (MQ=11)
           gtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATACGTGGGATCGGTCGGGGCTGTGGGATCGgtcgggtacgcgggat                                                                                               >  2:1650653/1‑95 (MQ=11)
            tccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtc                                                                                              <  1:1760653/101‑1 (MQ=11)
            tccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcggatc                                                                                              <  2:1010228/101‑1 (MQ=11)
             ccgtggggtcggtcggGTTCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtcc                                                                                             <  1:1593860/101‑1 (MQ=11)
             ccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGATCGgtcggatccgggggatcg                                                                                             >  2:4198239/1‑93 (MQ=11)
             ccgtgggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                                             <  1:100380/101‑1 (MQ=14)
             ccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTCGGATCCgcggggtccgtcgggtcc                                                                                             >  2:5237618/1‑101 (MQ=14)
              cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGtcgggtccgtggggtccg                                                                                            >  1:5161733/1‑101 (MQ=14)
              cgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCGGtggggtccggcgggtccg                                                                                            >  2:118550/1‑101 (MQ=11)
               gtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcggggccgcggggtccgt                                                                                           >  2:3482133/1‑101 (MQ=11)
                tggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtc                                                                                          >  2:2361374/1‑100 (MQ=14)
                tggggtccgtcggATCCGGGGGATCGGCCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGCGGGGTCGgccgggaccgtggggacggtc                                                                                          >  2:2474915/1‑95 (MQ=2)
                tgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTg                                                                                          >  1:954902/1‑101 (MQ=11)
                 gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCCGTgg                                                                                         >  2:5534468/1‑101 (MQ=11)
                 gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgg                                                                                         <  1:3227381/101‑1 (MQ=11)
                 cggttcggttggGTTCGTGGGATCGGTCTGGTCCGTCGGCTCGGGGGGATCCGTGGGATCTGTCGGGTCAGTGGGATCGgtcgggtccgtcggatcgggcg                                                                                         <  1:3304972/97‑6 (MQ=2)
                 aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTgg                                                                                         <  1:4361177/100‑1 (MQ=14)
                 aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTgg                                                                                         >  1:2777783/2‑101 (MQ=14)
                  gggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTggg                                                                                        <  2:5180967/92‑1 (MQ=14)
                  ggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTggg                                                                                        <  2:2924732/101‑1 (MQ=14)
                     tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtc                                                                                     >  1:4562948/1‑101 (MQ=14)
                     tcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGtc                                                                                     <  1:4500774/101‑1 (MQ=11)
                      ccgtgggATCGGTCGGATCCGTGGGGGCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTAGGAGCGGTCGGGGGCGTCGCGTCCGTGGGAtcg                                                                                    >  2:1537444/6‑101 (MQ=9)
                       ggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCGGTCGGATCCGTgggggccg                                                                                   >  2:389203/1‑96 (MQ=11)
                        gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCGGTGGGATCGGTGGGATCCGTGGGGTGCGTGGGAtcggt                                                                                  >  1:901419/1‑101 (MQ=1)
                        gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTGGGATCCGTGGGGGCCGTGGGAtcggt                                                                                  >  2:2076077/1‑101 (MQ=1)
                        gtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgt                                                                                  <  2:1857166/101‑1 (MQ=11)
                          cggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcgggcg                                                                                >  2:2139626/1‑101 (MQ=14)
                          cgcaTCCGTATGAGCGGGAGGGTCCGTCGGATCGGTCGGGTCCGTCGGATCTGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGtccgtcg                                                                                <  2:844931/98‑1 (MQ=2)
                           gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgg                                                                               >  1:260756/4‑101 (MQ=11)
                           gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                                               >  1:5404887/4‑101 (MQ=14)
                           ggATCCGTGGGGTCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtccgtggg                                                                               <  1:2978759/101‑1 (MQ=11)
                           cggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgg                                                                               <  1:2887801/98‑1 (MQ=11)
                            ggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTCGGGTCCGTCGGGTCCGTGGGGtcggtcggg                                                                              >  1:3972460/3‑100 (MQ=14)
                            gATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGa                                                                              <  1:3113207/101‑1 (MQ=11)
                             aTCCGTGGGATCGGTCGTATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                                                                             <  2:1273984/101‑1 (MQ=11)
                               ccggaggATCGGTCGGGTCCGGCGGATCGGTCGGATCCGTGGGCTCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                           <  2:5124492/96‑1 (MQ=14)
                                cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCg                                                                          >  2:2554194/1‑101 (MQ=14)
                                cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAGCGGTCGGGTCCg                                                                          >  2:3583581/1‑101 (MQ=14)
                                 gtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGt                                                                         <  1:5481442/98‑1 (MQ=18)
                                  taggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGtc                                                                        <  2:3753911/99‑2 (MQ=14)
                                   cggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                                       >  2:285289/2‑101 (MQ=255)
                                   aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGGTCGGTCGGGTCCGtcg                                                                       >  1:3876583/2‑99 (MQ=14)
                                    ggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTGGGGTCCGTGGGGTCGGTCGGGTCCGTcgg                                                                      >  2:3213754/1‑101 (MQ=14)
                                      gtcggtcggGTCGGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                    <  1:4668293/100‑1 (MQ=17)
                                      ctcggtcggGTCCGTTGGAGCGGTGGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGGGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                                                    <  1:130119/100‑1 (MQ=14)
                                      aTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                    >  2:4341076/1‑101 (MQ=14)
                                       tcgggcggGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCCGATcctacggacc                                                                   <  2:3972216/101‑9 (MQ=11)
                                        cggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                                  <  2:2230143/101‑1 (MQ=14)
                                          gtcggGGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                <  1:2066826/101‑1 (MQ=17)
                                           tcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtc                                                               <  2:4951314/101‑1 (MQ=14)
                                           tcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtc                                                               <  1:1914200/101‑1 (MQ=14)
                                             gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGCGGGAtcggtcgg                                                             >  1:3578755/1‑101 (MQ=17)
                                             ggatCCGTAGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg                                                             >  2:3065465/4‑101 (MQ=14)
                                              gttgcGTGGGATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCTGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCggg                                                            <  2:776428/97‑1 (MQ=14)
                                                 ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                         >  1:5651628/1‑101 (MQ=255)
                                                  cGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                        <  2:5392439/101‑1 (MQ=11)
                                                       gATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggg                                                   >  2:928554/1‑101 (MQ=255)
                                                        gtcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                                                  <  1:2709028/100‑3 (MQ=14)
                                                          cggtcggGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTcgggtcc                                                >  1:4199801/1‑98 (MQ=11)
                                                          cggtcggATCCGTGGGGTCGGTCGGGTCCGTCTGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcg                                                <  1:2423181/101‑5 (MQ=11)
                                                            gtcggGTCCGTGGGATCGGTCGGGTCCGGAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggt                                              <  2:1324330/101‑1 (MQ=255)
                                                             tcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtc                                             <  1:1680954/101‑8 (MQ=11)
                                                             tcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtg                                             >  1:1355056/1‑101 (MQ=255)
                                                             gcggATCCGTGGGTTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGTTCCGTGGGATCGGTCGGGTCCGTaggatcggtc                                             <  2:4469939/100‑8 (MQ=14)
                                                              cggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGACCGGTCGGGTCCGTGGGGCtggtgg                                            >  2:4260396/1‑101 (MQ=255)
                                                               gggtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTggg                                           >  2:697725/4‑101 (MQ=255)
                                                                gATCCGTGGGATCGGTGGTGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcggg                                          <  1:2285848/101‑11 (MQ=11)
                                                                gATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcggg                                          <  1:2144912/101‑11 (MQ=14)
                                                                gATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcggg                                          >  1:1242001/1‑91 (MQ=11)
                                                                  tCCGTCGGAACGGTCGGGTCGGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc                                        <  1:536096/101‑1 (MQ=255)
                                                                     ctGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggt                                     <  2:1770680/100‑1 (MQ=255)
                                                                      tGGGATCGGTCGGGTCGGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtc                                    <  1:198782/101‑1 (MQ=255)
                                                                      tGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGAACCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcgggc                                    >  1:5755836/1‑99 (MQ=255)
                                                                          aTCGGTCGGGTCCGTAGGATCGGTCGGGTCCATCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAc                                >  2:1938078/1‑101 (MQ=255)
                                                                             ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGc                             >  1:3061479/1‑101 (MQ=255)
                                                                              gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCt                            >  1:622402/1‑101 (MQ=255)
                                                                              gtcggGTCCGTAGGATCGGGCGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGGCCGCt                            <  2:5607446/101‑1 (MQ=255)
                                                                                cggGTCCGTGGGATTGGTCGGATCCGTGGGGTCCGTGGGATCGGTGGGGCCCGTGGGATCGGTCGGCCCCGTAGGGTCGGTGGGGCCGGTCGGGccagggg                          >  2:1456474/1‑96 (MQ=2)
                                                                                      cGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGGGGGATCGGCCGGACCGCTGTcgcccc                    >  2:5708314/1‑99 (MQ=255)
                                                                                            aTCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCgg              <  2:2681634/101‑1 (MQ=255)
                                                                                               ggtcggGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACa           <  2:5325871/101‑1 (MQ=255)
                                                                                                        cgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAggc  >  1:4833489/1‑101 (MQ=255)
                                                                                                        |                                                                                                    
GGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGC  >  NC_014034/221974‑222178

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG