Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,075 G→A 21.5% D194D (GAC→GAT RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0750G→A21.5% 226.3 ‑0.3 0.0410D194D (GAC→GATRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (43/45):  ref base (164/156):  total (208/202)
Fisher's exact test strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 703 total aligned reads displayed.

GGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCA  >  NC_014034/221973‑222175
                                                                                                      |                                                                                                    
   atcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTAGGGTCGGTGGGATCGgtggggg                                                                                                     >  2:3010551/1‑100 (MQ=14)
     cggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtcc                                                                                                   <  2:3053704/101‑1 (MQ=255)
      ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTGGGCTCGgtcgggtcgg                                                                                                  >  2:908137/1‑99 (MQ=11)
      ggtcgggtccgtgggagcggccggGGCGGTGGGATCGGGCGGGTCCGTGGGATGGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccg                                                                                                  <  1:2597964/101‑1 (MQ=11)
      ggtcgggtccggcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccg                                                                                                  <  1:2686170/101‑1 (MQ=11)
      ggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCGgtggggtacg                                                                                                  >  2:2821684/1‑101 (MQ=11)
       gtcgggtccgtggggtccgtcggATCCGTGGGATCGGCCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                                                                                                 <  1:3936528/101‑1 (MQ=14)
       gtcgactacgtggaataggtcggAGGCGTGGGATGGGTCGGGTCCGTGGGATCGGTCGGATCCGTGAGATCCGTCGTAACCGTGGGATCGgtcggatccgt                                                                                                 <  2:617351/93‑1 (MQ=0)
        tcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggggccgtg                                                                                                >  1:3873746/1‑100 (MQ=11)
        tcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTACGTGGGAGCGgtcggggccgtg                                                                                                >  2:1735947/1‑100 (MQ=11)
        tcgggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatccgtg                                                                                                <  2:1176909/101‑2 (MQ=11)
         gggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcg                                                                                               <  2:4435885/97‑1 (MQ=11)
         cgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgggg                                                                                               >  2:3468085/1‑98 (MQ=11)
         cgggtccgtgggatcgggcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcg                                                                                               <  1:3283602/101‑1 (MQ=11)
         cgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGggcgggtccgtcg                                                                                               >  2:2539762/1‑101 (MQ=11)
         cggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcg                                                                                               >  1:689657/1‑101 (MQ=255)
         cggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcg                                                                                               <  1:3065660/101‑1 (MQ=11)
          gggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggttcggagg                                                                                              >  2:236226/1‑97 (MQ=11)
          ggatccgtggggtccgtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGggggggtccgtcgg                                                                                              >  2:586071/4‑101 (MQ=11)
          ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                                              >  1:3760059/4‑101 (MQ=255)
          ggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGgtcgggtccgtcgg                                                                                              >  2:2634046/4‑101 (MQ=11)
          ggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGgtcgggtccgtcgg                                                                                              >  1:2346802/4‑101 (MQ=11)
           ggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGGCCGTGGGGTCCgtcgggtcccggggg                                                                                             >  1:4397289/1‑95 (MQ=11)
           gatccgtggtatcgggcggATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGCTCCGTAGGATCGGTCGGGTCCGTCGGATCGgtcggatccgtggga                                                                                             <  2:3132974/99‑2 (MQ=11)
           gatccgtaggatgggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                                                             <  1:1003866/99‑1 (MQ=11)
            atccgtaggtgcggtcggGTCCGTAGTATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggt                                                                                            <  2:1466833/100‑1 (MQ=11)
              ccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggatcg                                                                                          <  2:1253427/101‑2 (MQ=255)
              ccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                                          >  1:2753295/1‑101 (MQ=14)
                gtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGATCGGTcgggtccgtggggtccgt                                                                                        >  1:2216112/1‑101 (MQ=11)
                gtgggatcggtcggATCCGTGGTATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGCCGGGACCGTAGGATCGgtcgggtccgtgggatcggc                                                                                        >  1:3133218/1‑98 (MQ=11)
                gtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtggggtcggt                                                                                        >  1:1820643/1‑101 (MQ=14)
                gtggcgtcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgt                                                                                        <  2:53305/95‑1 (MQ=11)
                 tggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTg                                                                                       >  2:2242310/1‑101 (MQ=11)
                 tgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGGTCGGGCGGAGCCGTTGGGTCCG.CGGGGCCGGgg                                                                                      >  2:1975116/1‑101 (MQ=11)
                 tgggagcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggtc                                                                                       <  2:1354084/101‑2 (MQ=11)
                 taggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTg                                                                                       <  2:2703186/99‑1 (MQ=14)
                  ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtcg                                                                                      >  1:120762/1‑99 (MQ=11)
                  aggatcggtcggGTCCGTAGGATCGGTCAGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGCCCGGgg                                                                                      >  2:2889398/2‑101 (MQ=11)
                    ggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGa                                                                                    <  2:3180634/96‑1 (MQ=11)
                    gatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGa                                                                                    <  1:4239285/101‑1 (MQ=14)
                     gtcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAt                                                                                   <  1:1492107/100‑1 (MQ=11)
                     gtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATAGGTCGGGTCCGTCGGGCTCGTGGGAt                                                                                   >  2:722339/8‑101 (MQ=11)
                     gtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAt                                                                                   >  2:3327019/5‑101 (MQ=11)
                        cgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCAGGATCGGTCGGGCCCGTCGGGGGCGTGGCAccgg                                                                                >  2:3613420/5‑101 (MQ=11)
                         gtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggt                                                                               <  1:1731346/98‑1 (MQ=14)
                          tgggATCTGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtc                                                                              <  2:4160421/99‑1 (MQ=17)
                          tcggATCCGTGGGATCGGTCGGATCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGAGCGGTCGGGGCCCGGGGGGCCGTCGGAtccgtg                                                                              >  2:2924872/1‑100 (MQ=11)
                           cggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcg                                                                             <  2:3925250/101‑1 (MQ=11)
                            ggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                                            <  2:4319772/101‑1 (MQ=17)
                            ggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTGGGGtcggccgg                                                                            >  2:4001580/1‑101 (MQ=14)
                            ggATCGGTCGGATCCGTGGGGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                                            <  1:4666147/101‑1 (MQ=17)
                             ggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcggg                                                                           >  1:553197/3‑100 (MQ=11)
                              gtcggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGAt                                                                          <  2:4651814/94‑1 (MQ=17)
                               gccGTGGGGTCCGGCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATc                                                                         <  2:4489470/100‑1 (MQ=11)
                                 ggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCg                                                                       <  1:354585/97‑1 (MQ=18)
                                  gtcggATCCGTGGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGt                                                                      <  1:4501247/98‑1 (MQ=17)
                                  gTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGGTCGGGCGGGTGCGTGGGATCGGGGGGGTGGGTGGGacaggg                                                                      >  2:1583326/1‑96 (MQ=0)
                                    gggATGGGTCGAGCCGGCGGGATCGGTTGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                    <  1:2834746/101‑1 (MQ=14)
                                    gggATCGGTCGGGTCGGTGGGATCGGTCGGATCGGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                    <  1:908217/101‑1 (MQ=14)
                                    gggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                    >  1:3414375/1‑101 (MQ=14)
                                      gATCGGTCGGGTCCGTGGGATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                                                  <  2:3886071/101‑1 (MQ=17)
                                       aTCCGTGGGGTTCGTGGGATCGGTCGGGTCCGTGGGATCCGTCGGGTCCGTAGGATCGGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                                                 <  2:2472939/101‑1 (MQ=18)
                                        tccgtggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGGTCGGTCGGATCCGTCGGAtg                                                                >  1:241886/7‑100 (MQ=17)
                                        gcggtcggGTCCGGCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtc                                                                <  2:1910660/100‑1 (MQ=14)
                                         cggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATTGGTGGGATCGGTGGGTTCGGTGGGACCGCTCGGGGCCGGGGGAccg                                                               >  2:4208307/1‑101 (MQ=11)
                                         cggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTGGGGTCCGTCGGAtcg                                                               >  2:1040878/1‑101 (MQ=14)
                                          ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                              <  1:61080/101‑1 (MQ=17)
                                            tggggTCCGCGGGGTGGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAGCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtc                                                            <  1:1058716/99‑1 (MQ=17)
                                             cggGTCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                           <  1:4342012/101‑1 (MQ=17)
                                               ggTCCGTGGGATCCGTCGGATCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCggg                                                         <  2:4081237/101‑1 (MQ=17)
                                                  ccGTGGGATCGGGCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                      <  1:197149/101‑1 (MQ=11)
                                                    gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCgg                                                    >  1:3216504/1‑100 (MQ=255)
                                                    gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                                                    >  1:357352/1‑101 (MQ=255)
                                                    gTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                    >  2:2156011/1‑101 (MQ=17)
                                                     tcggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTg                                                   >  1:2128067/3‑101 (MQ=17)
                                                      gggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                                                  <  1:73777/101‑1 (MQ=255)
                                                       ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggg                                                 <  2:2764530/101‑1 (MQ=17)
                                                        gATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTcggg                                                <  1:3659409/101‑1 (MQ=11)
                                                           cggtcggATCCTTGCGATCGGTGGGGTCCGTGGCATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcg                                             <  1:4023107/101‑5 (MQ=14)
                                                            gatctgATCCGGGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcgg                                            <  1:4101320/96‑6 (MQ=11)
                                                              tcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtc                                          <  1:4587335/101‑8 (MQ=11)
                                                              gcggATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtc                                          <  2:2034848/100‑8 (MQ=14)
                                                               cggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGGggggacggggg                                         >  2:3730726/1‑94 (MQ=14)
                                                                  gtCCCTGGGGTCCGTCGGGTCGGTAGGATTGGTCGGGTCCGTGGGGTCGGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCCtgctggcga                                      <  2:3259017/100‑4 (MQ=17)
                                                                     cGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcgg                                   >  1:2017267/1‑101 (MQ=255)
                                                                      gTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggt                                  >  1:2861731/1‑101 (MQ=255)
                                                                       ttggACCGGTGGGGCCGGTGGCATCGGTCGAGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtctgtc                                 <  2:49212/99‑1 (MQ=18)
                                                                       tGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTTGTGGGAtcggtc                                 <  1:2522619/101‑1 (MQ=255)
                                                                        gggGTCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcg                                <  1:3010780/101‑1 (MQ=255)
                                                                           aTCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAc                             >  1:1266687/1‑101 (MQ=255)
                                                                              ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGc                          >  1:2950663/1‑101 (MQ=255)
                                                                                  gggTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTGGGGCTGGTGGGGTCGGTCGGACAGCTggc                      >  2:848169/1‑99 (MQ=255)
                                                                                     tCCGTAGGGTCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgcc                   <  1:1391020/101‑1 (MQ=255)
                                                                                     tCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgcc                   <  2:1673343/101‑1 (MQ=255)
                                                                                     tCCGTAGGATCGGTCGGGTCCGTCGGGTCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgcc                   <  1:4162808/101‑1 (MQ=255)
                                                                                       cttgggctgtgtcgggtccgtcgggtccGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgc                 <  2:3414121/91‑1 (MQ=255)
                                                                                        gTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgcc                >  1:1649984/1‑101 (MQ=255)
                                                                                                  tcgggtcggtggggtccgTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGa      <  1:859293/101‑1 (MQ=255)
                                                                                                    gggtccgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGAcc    >  1:1558062/1‑101 (MQ=255)
                                                                                                    gggtccgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGAcc    >  2:2406909/1‑101 (MQ=255)
                                                                                                      |                                                                                                    
GGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCA  >  NC_014034/221973‑222175

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