Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,949 C→A 5.4% E236D (GAG→GAT RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9490C→A 5.4% 18.1 ‑3.1 ‑1.8187E236D (GAG→GATRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/1):  ref base (95/81):  total (104/83)
Fisher's exact test strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.75e-02

Only 100 of 223 total aligned reads displayed.

AGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221849‑222049
                                                                                                    |                                                                                                    
 gCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCt                                                                                                     <  2:4043345/101‑1 (MQ=255)
 gCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCt                                                                                                     >  2:3001175/1‑101 (MQ=255)
  ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                                    >  2:2452152/1‑101 (MQ=255)
  ccGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtc                                                                                                    >  2:3553496/1‑101 (MQ=255)
    ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcgg                                                                                                  <  1:1611983/101‑1 (MQ=255)
    ggTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGGGGGCtcgg                                                                                                  >  1:2332632/1‑101 (MQ=255)
     gTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggt                                                                                                 <  1:324169/101‑1 (MQ=255)
      tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                                >  2:1758597/1‑101 (MQ=255)
      tttCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCtcggtc                                                                                                <  1:922756/101‑1 (MQ=255)
        tCGCCGGTGCCCGACACCGTCGACGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTAGGCTGCGTTGGGtcggtcgg                                                                                              >  2:1114265/1‑101 (MQ=255)
           ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCTCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                           <  2:2118579/101‑1 (MQ=255)
           ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                           <  1:182344/101‑1 (MQ=255)
           ccGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatc                                                                                           <  1:1489538/101‑1 (MQ=255)
            cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                          >  1:730281/1‑101 (MQ=255)
            cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                          >  2:3366552/1‑101 (MQ=255)
            cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatca                                                                                          >  1:2114522/1‑100 (MQ=255)
             ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                         >  2:3130278/1‑101 (MQ=255)
             ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                         >  2:3527668/1‑101 (MQ=255)
             ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                         <  1:3522368/101‑1 (MQ=255)
               tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                       <  2:2686388/101‑1 (MQ=255)
               tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                       <  2:3770876/101‑1 (MQ=255)
                 cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgg                                                                                     >  2:1773162/1‑101 (MQ=255)
                  ccGACACCGTCGCCGAGGCTTCGGTGTTGGGGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcggg                                                                                    >  1:3763635/1‑101 (MQ=255)
                   cGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggt                                                                                   >  1:2224625/1‑101 (MQ=255)
                     gcacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  1:3292170/100‑1 (MQ=255)
                     acacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGGGTGGGCTCGGTCggatcggtcgggtcc                                                                                 <  2:3351481/101‑1 (MQ=255)
                      cacCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccg                                                                                >  1:1969809/1‑101 (MQ=255)
                         cGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgg                                                                             <  2:1068697/101‑1 (MQ=255)
                          ggcGCCGGGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                            <  2:2766727/99‑1 (MQ=255)
                          gTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggg                                                                            >  1:3437845/1‑101 (MQ=255)
                           tCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtggga                                                                           <  1:3735014/101‑1 (MQ=255)
                                gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggt                                                                      >  1:102121/1‑101 (MQ=255)
                                gAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggagcggg                                                                      >  2:1960184/1‑100 (MQ=255)
                                gAGGCTTCGGTGTTGGAGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgaggccgtgggctcggc                                                                      >  2:2583725/1‑100 (MQ=255)
                                 aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                     >  1:2641239/1‑101 (MQ=255)
                                  ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcg                                                                    <  1:4116663/101‑1 (MQ=255)
                                     ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGTCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                 >  1:1326319/1‑101 (MQ=255)
                                     ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                 <  2:3080498/101‑1 (MQ=255)
                                     ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                 >  2:773034/1‑101 (MQ=255)
                                     ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                 >  1:2775774/1‑101 (MQ=255)
                                        ggTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccg                                                              >  2:1388500/1‑101 (MQ=255)
                                           gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                           >  2:3047793/1‑101 (MQ=255)
                                             tggtggTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAtcggtcgggtccgtgggatcggtcgggtccgtggga                                                         >  1:327101/1‑101 (MQ=255)
                                               gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                       <  1:2754564/101‑1 (MQ=255)
                                                tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGCTCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                      >  1:3182088/1‑101 (MQ=255)
                                                tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGgtcgggtccgtgggatcggtcggggccgtgggatcg                                                      >  1:2708375/1‑101 (MQ=255)
                                                   tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcggggccgtgggatcggtc                                                   >  1:999414/1‑101 (MQ=255)
                                                      gCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGa                                                >  1:3108231/1‑101 (MQ=255)
                                                           gCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCgg                                           >  1:1130117/1‑100 (MQ=255)
                                                           gCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt                                           >  2:1317439/1‑101 (MQ=255)
                                                           gCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCCGCTCGGTCGGGGCCGTGGGACCGGGCGGGTCCGGGGGGTCGGTCGGATCCGt                                           >  1:1079973/1‑101 (MQ=255)
                                                            ccGGATCGCGCTCATCGACGGTGGGCGTCGGATACGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                          <  1:3628306/101‑1 (MQ=255)
                                                               gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTgggg                                       >  2:1546483/1‑100 (MQ=255)
                                                               gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                       <  1:347727/101‑1 (MQ=255)
                                                                 tCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGTCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                     <  2:1326219/101‑1 (MQ=255)
                                                                      cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                >  2:1564479/1‑101 (MQ=255)
                                                                       tCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                               >  2:3687661/1‑101 (MQ=255)
                                                                        cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg                              >  1:1876909/1‑101 (MQ=255)
                                                                          tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCggggc                            >  2:3735724/1‑99 (MQ=255)
                                                                          tCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                            <  1:1507999/101‑1 (MQ=255)
                                                                           cGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                           >  1:252169/1‑101 (MQ=255)
                                                                              cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                        <  2:615902/101‑1 (MQ=255)
                                                                               ggTGGGCGTCGGATGCGTGGGGTGGGTCGGTTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAACCGTGGGATCGGTCGGGTCCGTgg                       <  2:44776/101‑1 (MQ=18)
                                                                               ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                       >  1:1307806/1‑101 (MQ=255)
                                                                               ggTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                       >  2:3863140/1‑101 (MQ=255)
                                                                                gTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                      >  1:2962715/1‑101 (MQ=255)
                                                                                   ggggTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                   >  1:2413805/4‑101 (MQ=255)
                                                                                   ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCGGTGGGGtc                   >  2:1901930/1‑101 (MQ=255)
                                                                                   ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                   >  1:4196785/1‑101 (MQ=255)
                                                                                   ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                   <  1:2052160/101‑1 (MQ=255)
                                                                                     ggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccg                 >  2:2213199/2‑99 (MQ=11)
                                                                                     cGTCGGATCTGTGGGATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                 <  1:1884945/101‑1 (MQ=11)
                                                                                      gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                >  1:1434621/1‑101 (MQ=255)
                                                                                      gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                <  1:2042472/101‑1 (MQ=255)
                                                                                      gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                >  2:742564/1‑101 (MQ=255)
                                                                                       tggcaTGCGTAGGGTCGCGCGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  2:3590610/97‑1 (MQ=14)
                                                                                       tCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGtccgtg               <  2:2440845/101‑2 (MQ=11)
                                                                                       gcGGATGCGTGGGCTCGGTCGGTTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc               <  2:714899/100‑1 (MQ=255)
                                                                                        cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg              <  1:1925897/101‑1 (MQ=255)
                                                                                        cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg              >  2:3471668/1‑101 (MQ=255)
                                                                                          gATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGTTCGGa            <  2:2622138/101‑1 (MQ=255)
                                                                                            tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc          <  2:2113937/101‑1 (MQ=255)
                                                                                             tcGTGGGCGCGGTCGGATGTGGGGGGTGCGTGGTATCGGTCGGGTCCGTGGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc         <  1:1813002/100‑1 (MQ=14)
                                                                                             gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATcc         >  2:2075183/1‑101 (MQ=255)
                                                                                             gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc         <  2:3257619/101‑1 (MQ=255)
                                                                                             ccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc         <  1:2990604/100‑1 (MQ=1)
                                                                                              cGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg        <  2:330468/101‑1 (MQ=11)
                                                                                               gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt       <  1:3466737/101‑1 (MQ=255)
                                                                                               gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGATCCGTGGGATCGGACGGGTCGGTGGGATGGGTCGGATCCGt       >  2:1242439/1‑101 (MQ=17)
                                                                                               gTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGt       >  2:442325/1‑101 (MQ=1)
                                                                                                tGGGCTCGGTCGGATCGGTCGGGTCGCTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      <  1:3553749/101‑1 (MQ=255)
                                                                                                tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg      >  2:434983/1‑101 (MQ=255)
                                                                                                tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcggtcgggggcgtc      >  1:969407/1‑94 (MQ=11)
                                                                                                tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCggc      >  2:227711/1‑99 (MQ=11)
                                                                                                tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtc      <  2:3926693/101‑2 (MQ=11)
                                                                                                tGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtc      >  1:2909741/1‑100 (MQ=11)
                                                                                                tGGGATCGGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGtc      <  2:3157492/101‑2 (MQ=11)
                                                                                                  ggCTCGGCAGGATGGGTCGGGCCCGTGGTAAGGGTCGGGTACGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg    <  1:52978/101‑1 (MQ=14)
                                                                                                   gCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa   >  1:395282/1‑101 (MQ=255)
                                                                                                    cTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt  >  1:1755928/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
AGCCGGTTTCGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAT  >  NC_014034/221849‑222049

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG