Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,293,808 +G 1.1% coding (191/1707 nt) ureC → urease subunit alpha

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,293,8081.→G 1.1% 3.9 ‑2.0 ‑0.7363coding (191/1707 nt)ureCurease subunit alpha
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (197/162):  total (201/162)
Fisher's exact test strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.95e-01

Only 100 of 398 total aligned reads displayed.

TCGAGGTTGAAAAAGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATCGCGGGG  >  NC_014034/1293709‑1293908
                                                                                                    |                                                                                                    
  gAGGTTGAAAAAGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGCC.ggg                                                                                                   >  2:269820/1‑101 (MQ=255)
     gTTGAAAAAGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.gggggg                                                                                                >  2:349406/1‑101 (MQ=255)
      ttGAAAAAGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.ggggggg                                                                                               >  2:1497413/1‑101 (MQ=255)
       tGAAAAAGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.gggggggg                                                                                              >  1:2150597/1‑101 (MQ=255)
         aaaaaGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCg                                                                                            <  2:6546935/101‑1 (MQ=255)
         aaaaaGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCg                                                                                            >  1:4663247/1‑101 (MQ=255)
             aGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATgg                                                                                        <  2:3303636/101‑1 (MQ=255)
                tctcACACCCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGaca                                                                                     <  2:5495767/101‑1 (MQ=255)
                  tcACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACAcc                                                                                   <  1:6410156/101‑1 (MQ=255)
                  tcACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACAcc                                                                                   <  1:4795376/101‑1 (MQ=255)
                  tcACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACAcc                                                                                   >  2:858462/1‑101 (MQ=255)
                   cacaACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCg                                                                                  >  1:2065843/1‑101 (MQ=255)
                    acaACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGt                                                                                 >  1:6158741/1‑101 (MQ=255)
                    acaACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGt                                                                                 >  1:368272/1‑101 (MQ=255)
                     caACCTACGGAGCAGAGGTGAAATTCGGCGGCGGCAAAGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtc                                                                                >  1:6321549/1‑101 (MQ=255)
                     caACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtc                                                                                >  2:2188066/1‑101 (MQ=255)
                     caACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtc                                                                                <  1:5747272/101‑1 (MQ=255)
                     caACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtc                                                                                <  1:5863674/101‑1 (MQ=255)
                     caACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtc                                                                                <  1:1296823/101‑1 (MQ=255)
                      aaCCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtca                                                                               <  2:2515648/101‑1 (MQ=255)
                         cTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtcatca                                                                            >  1:421535/1‑101 (MQ=255)
                         cTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtcatca                                                                            >  1:1078799/1‑101 (MQ=255)
                         cTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGtcatca                                                                            <  1:4506169/101‑1 (MQ=255)
                         cTACGGAGAAGAGGTGAAATCCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGGTGGACACCGgcatca                                                                            <  1:5630407/101‑1 (MQ=255)
                            cGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGGCGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCa                                                                         <  1:4121162/101‑1 (MQ=255)
                            cGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCa                                                                         <  2:3651737/101‑1 (MQ=255)
                            cGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCa                                                                         <  2:4827714/101‑1 (MQ=255)
                             ggAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GTGGGGGGCGATGGACACCGTCATCACCaa                                                                        <  2:2399141/101‑1 (MQ=255)
                             ggAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCaa                                                                        >  2:710155/1‑101 (MQ=255)
                              gagaAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAAt                                                                       >  1:2338109/1‑101 (MQ=255)
                               agaagaTGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATg                                                                      <  1:4048392/101‑1 (MQ=255)
                               agaagaGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATg                                                                      <  1:1944260/101‑1 (MQ=255)
                               agaagaGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATg                                                                      >  2:1100050/1‑101 (MQ=255)
                               agaagaGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATg                                                                      <  2:1333543/101‑1 (MQ=255)
                                gaagaGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATgc                                                                     >  2:3895724/1‑101 (MQ=255)
                                gaagaGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATgc                                                                     >  1:5438859/1‑101 (MQ=255)
                                   gagGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCt                                                                  >  2:173730/1‑101 (MQ=255)
                                   gagGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCt                                                                  >  2:6773925/1‑101 (MQ=255)
                                    agGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTg                                                                 >  2:6082681/1‑101 (MQ=255)
                                     ggTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGa                                                                >  2:4737427/1‑101 (MQ=255)
                                      gTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGAt                                                               >  1:3224686/1‑101 (MQ=255)
                                       tGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGAtc                                                              <  2:5594522/101‑1 (MQ=255)
                                        gAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGAtcg                                                             >  2:2341078/1‑101 (MQ=255)
                                        gAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGAtcg                                                             <  2:4625003/101‑1 (MQ=255)
                                           aTTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGAtcgtcg                                                          >  1:3083150/1‑101 (MQ=255)
                                            ttCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGa                                                         <  1:5290580/101‑1 (MQ=255)
                                            ttCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGa                                                         <  1:2642295/101‑1 (MQ=255)
                                              cggcggcggCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATc                                                       <  2:4167480/101‑1 (MQ=255)
                                               ggcggcggcAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATca                                                      >  2:2831454/1‑101 (MQ=255)
                                                  ggcggcAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCacac                                                   <  1:1734075/101‑1 (MQ=255)
                                                  ggcggcAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCacac                                                   <  1:3539395/101‑1 (MQ=255)
                                                    cggcAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACgg                                                 <  1:5349278/101‑1 (MQ=255)
                                                     ggcAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTAATCACCAATGCGCTGATCGCCGATCAAACggc                                                >  2:5251524/1‑100 (MQ=255)
                                                       cAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCa                                              >  1:760919/1‑101 (MQ=255)
                                                       cAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCa                                              >  2:4014543/1‑101 (MQ=255)
                                                          ggTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCt                                           <  2:6321224/101‑1 (MQ=255)
                                                           gTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTa                                          <  2:295648/101‑1 (MQ=255)
                                                            tGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTAc                                         <  1:4285872/101‑1 (MQ=255)
                                                             gATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACa                                        >  2:2077643/1‑101 (MQ=255)
                                                                ccGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAgg                                     <  2:3748879/101‑1 (MQ=255)
                                                                ccGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAgg                                     <  2:4249195/101‑1 (MQ=255)
                                                                  gcgACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGcc                                   <  2:6206620/101‑1 (MQ=255)
                                                                    gACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGa                                 >  2:3140981/1‑101 (MQ=255)
                                                                    gACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGa                                 >  1:3462130/1‑101 (MQ=255)
                                                                    gACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGa                                 <  2:193782/101‑1 (MQ=255)
                                                                       ggCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATgg                              >  1:1041395/1‑100 (MQ=255)
                                                                       ggCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGt                              >  2:1682575/1‑101 (MQ=255)
                                                                        gCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGCTCACACGGGCATCTAGAAGGCCGATGTc                             <  1:5063548/101‑1 (MQ=255)
                                                                        gCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTc                             >  1:2453351/1‑101 (MQ=255)
                                                                        gCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTc                             <  1:269824/101‑1 (MQ=255)
                                                                        gCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTc                             >  1:2792081/1‑101 (MQ=255)
                                                                        gCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTc                             >  1:6818971/1‑101 (MQ=255)
                                                                         cATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCg                            >  2:695968/1‑101 (MQ=255)
                                                                          aTGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCgg                           >  1:1532397/1‑101 (MQ=255)
                                                                           tGGGGCAATCGCAAATCCCGCGGTCGGGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCgg                           >  2:44233/1‑101 (MQ=255)
                                                                           tGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGt                          <  2:4063973/101‑1 (MQ=255)
                                                                           tGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGt                          >  1:4925032/1‑101 (MQ=255)
                                                                           tGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGt                          <  1:4591340/101‑1 (MQ=255)
                                                                             gggCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCt                        >  1:1370691/1‑101 (MQ=255)
                                                                              ggCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTg                       <  2:6472907/101‑1 (MQ=255)
                                                                                 aaTCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTgcgc                    >  1:2864315/1‑101 (MQ=255)
                                                                                  aTCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGcgcg                   <  2:3025341/101‑1 (MQ=255)
                                                                                    cGCAAAGCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGAt                 <  1:6077190/101‑1 (MQ=255)
                                                                                      cAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATgg               <  2:120510/101‑1 (MQ=255)
                                                                                       aaaTCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGa              >  1:871114/1‑101 (MQ=255)
                                                                                       aaaTCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGACGTCGGTCTGCGCGATGGa              >  1:6982038/1‑101 (MQ=255)
                                                                                       aaaTCCCGCGGTC.GGGGGGGGCGATGGACACCGGCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGGCTGCGCGATGGa              >  2:6675019/1‑101 (MQ=255)
                                                                                           cccGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGa          >  1:6596886/1‑101 (MQ=255)
                                                                                            ccGCGGTCGGGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGa          >  1:2572875/1‑101 (MQ=255)
                                                                                             cgcgGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATc        <  1:970464/101‑1 (MQ=255)
                                                                                              gcgGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATcg       >  1:2940798/1‑101 (MQ=255)
                                                                                              gcgGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATcg       <  1:5933977/101‑1 (MQ=255)
                                                                                               gggTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATcgc      <  2:1892885/100‑1 (MQ=255)
                                                                                               cgGTC.GGGGGGGGCGATGGACACCGTCATCCCCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATcgc      <  1:6670958/101‑1 (MQ=255)
                                                                                               cgGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATcgc      >  1:1173591/1‑101 (MQ=255)
                                                                                                ggTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACGAGGCCGATGTCGGTCTGCGCGATGGACGGATcgcg     >  1:3374686/1‑101 (MQ=255)
                                                                                                ggTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATcgcg     >  1:3835278/1‑101 (MQ=255)
                                                                                                ggTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATcgcg     >  2:3598901/1‑101 (MQ=255)
                                                                                                  tC.GGGGGGGGCGATCGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATCGCggg   <  2:5885183/101‑1 (MQ=255)
                                                                                                   c.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATCGCgggg  >  1:6217904/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
TCGAGGTTGAAAAAGATCTCACAACCTACGGAGAAGAGGTGAAATTCGGCGGCGGCAAGGTGATCCGCGACGGCATGGGGCAATCGCAAATCCCGCGGTC.GGGGGGGGCGATGGACACCGTCATCACCAATGCGCTGATCGTCGATCACACGGGCATCTACAAGGCCGATGTCGGTCTGCGCGATGGACGGATCGCGGGG  >  NC_014034/1293709‑1293908

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG