Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,561 C→G 14.6% P148R (CCG→CGG)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,5610C→G14.6% 107.0 ‑2.0 ‑0.3302P148R (CCG→CGG) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (18/26):  ref base (139/119):  total (157/145)
Fisher's exact test strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.59e-01

Only 100 of 348 total aligned reads displayed.

GATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCGT  >  NC_014034/1396461‑1396660
                                                                                                     |                                                                                                   
gATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGGACGTCGCTCCCCCCGCGGGGGCACGGATCGAAACTGCTTTC.GTGGGACTGGCCAgcc                                                                                                     >  1:4944646/1‑101 (MQ=255)
gATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGACGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgcc                                                                                                     <  2:5235006/101‑1 (MQ=255)
 aTCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAgccg                                                                                                    <  2:3624056/101‑1 (MQ=255)
    gAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGccgccc                                                                                                 >  1:4870810/1‑101 (MQ=255)
     aaaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGccgcccg                                                                                                >  1:177192/1‑101 (MQ=255)
     aaaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGccgcccg                                                                                                <  1:2238461/101‑1 (MQ=255)
      aaCGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGccgcccgc                                                                                               >  2:5770246/1‑101 (MQ=255)
        cGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCgcgc                                                                                             >  1:2536103/1‑101 (MQ=255)
          gACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCa                                                                                           >  2:2531954/1‑101 (MQ=255)
            ccTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAg                                                                                         <  1:3803614/101‑1 (MQ=255)
             cTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGt                                                                                        <  1:523379/101‑1 (MQ=255)
               gAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCg                                                                                      >  1:1333043/1‑101 (MQ=255)
                aaaTTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGc                                                                                     >  2:3227789/1‑101 (MQ=255)
                    tCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTccc                                                                                 >  2:2963900/1‑101 (MQ=255)
                     ccGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTcccc                                                                                <  1:4190277/101‑1 (MQ=255)
                       gcgGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTcccccc                                                                              <  1:425523/101‑1 (MQ=255)
                       gcgGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTcccccc                                                                              <  1:2861442/101‑1 (MQ=255)
                       gcgGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTcccccc                                                                              >  1:1037161/1‑101 (MQ=255)
                          ggAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCcgcg                                                                           <  2:3996520/101‑1 (MQ=255)
                            aaaCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCggg                                                                         >  1:1107208/1‑101 (MQ=255)
                            aaaCGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCggg                                                                         >  1:3953004/1‑101 (MQ=255)
                            aaaCGCGGCGCCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGCAGCTCCACCCGCggg                                                                         >  1:4997425/1‑101 (MQ=37)
                               cgcgGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGc                                                                      <  1:4424071/101‑1 (MQ=255)
                               cgcgGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGc                                                                      <  1:708033/101‑1 (MQ=255)
                                gcggcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCg                                                                     >  1:3034521/1‑101 (MQ=255)
                                gcggcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCg                                                                     >  2:713690/1‑101 (MQ=255)
                                  ggcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCggg                                                                   >  1:1968921/1‑99 (MQ=255)
                                  ggcggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTg                                                                   <  1:5487474/101‑1 (MQ=255)
                                    cggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                 >  2:5578806/1‑101 (MQ=255)
                                     ggcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                <  2:509823/101‑1 (MQ=255)
                                      gcgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                               <  1:4603931/101‑1 (MQ=255)
                                       ggacggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATc                                                               <  1:1435881/94‑1 (MQ=2)
                                       cgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATca                                                              <  2:4765783/101‑2 (MQ=255)
                                       cgAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                              >  1:84896/1‑101 (MQ=255)
                                        gAGCACGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                             <  2:3560033/101‑1 (MQ=255)
                                         aGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                            >  2:1850235/1‑101 (MQ=255)
                                          cggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGaa                                                            >  2:1159602/5‑101 (MQ=2)
                                           ggaggcGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGaaa                                                           <  2:1755427/96‑1 (MQ=2)
                                           ggagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGaaa                                                           >  2:1358472/4‑101 (MQ=2)
                                           caggTCACTCCCGCCGCGGGGGGGGGAATCGAAACTGATTAC.GTGGGCCTGGCCACCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          <  2:4862953/98‑1 (MQ=21)
                                           cACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                          <  2:679546/101‑1 (MQ=255)
                                            gagTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGGTCGAAAc                                                          >  2:8490/3‑101 (MQ=2)
                                            aCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACg                                                         <  2:3471147/101‑1 (MQ=255)
                                            aCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCCCGGGGGCGTGGATCGAAACg                                                         >  2:2563958/1‑101 (MQ=255)
                                             cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGa                                                        >  1:5312097/1‑101 (MQ=255)
                                             cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGa                                                        <  2:2069779/101‑1 (MQ=255)
                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAt                                                       >  1:1929139/1‑101 (MQ=255)
                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAtt                                                      >  1:5073883/1‑101 (MQ=255)
                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgt                                                       <  1:2439447/101‑4 (MQ=2)
                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                                     <  1:2964067/101‑1 (MQ=255)
                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                                     <  1:642726/101‑1 (MQ=255)
                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGAttt                                                     <  1:1675793/101‑1 (MQ=255)
                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtg                                                      >  1:2177471/1‑97 (MQ=2)
                                                  cTCCCCCGGCGGGGGCGTGGATCGAAACTTCTTTC.GTGGGACTGGCCACCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcg                                                   <  1:2932799/101‑1 (MQ=255)
                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcg                                                   >  2:1537911/1‑101 (MQ=255)
                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcg                                                   <  2:5647171/101‑1 (MQ=255)
                                                   tCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccg                                                   >  2:1284337/1‑94 (MQ=2)
                                                   tCCCCCCGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccg                                                   >  2:4046364/1‑94 (MQ=2)
                                                    ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgtc                                                 <  1:1830719/101‑1 (MQ=255)
                                                    ccccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTtcgtc                                                 <  2:3325059/101‑1 (MQ=255)
                                                    ccccccGCGGGGGCGTGGATCGAAACAGCTTTCAATGGGGGCGGGCGCCGGGCCAACCTATGTCGCTCACCCCGCGGGGGCGTGGATCGAAActgtgccgg                                                  <  1:3640680/101‑9 (MQ=2)
                                                      ccccGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGc                                               >  2:3069464/1‑101 (MQ=255)
                                                         cgcgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTa                                            <  2:4937605/101‑1 (MQ=255)
                                                          gcgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACGATTTCGTCGCTTac                                           >  2:4226252/1‑101 (MQ=255)
                                                          gcgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACGATTTCGTCGGTTac                                           >  1:784584/1‑101 (MQ=255)
                                                          gcgGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTac                                           <  1:5614687/101‑1 (MQ=255)
                                                            gggggCGGGGAGCGAAGCTCCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTacga                                         <  2:4508750/101‑1 (MQ=38)
                                                                 cGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCt                                    >  2:429627/1‑101 (MQ=255)
                                                                  gTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTa                                   >  1:2418782/1‑101 (MQ=255)
                                                                   tGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTAt                                  >  1:5133408/1‑101 (MQ=255)
                                                                   tGGATCGAAACTGCTTGC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTAt                                  <  1:2035679/101‑1 (MQ=255)
                                                                    ggATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATc                                 >  1:541826/1‑101 (MQ=255)
                                                                    ggATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATc                                 <  2:2434462/101‑1 (MQ=255)
                                                                    ggATCGAAACTGCTTTC.GTGGGACTGACCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATc                                 <  2:3675600/101‑1 (MQ=255)
                                                                     gATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCt                                <  1:35182/101‑1 (MQ=255)
                                                                     gATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCt                                >  2:4588188/1‑101 (MQ=255)
                                                                      aTCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTg                               >  2:5795011/1‑101 (MQ=255)
                                                                      aTCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTg                               <  1:3485836/101‑1 (MQ=255)
                                                                         gAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAtt                            >  1:983601/1‑101 (MQ=255)
                                                                          aaaCTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAttt                           <  1:1518244/101‑1 (MQ=255)
                                                                          aaaCTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGAttt                           >  2:2733531/1‑101 (MQ=255)
                                                                             cTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCgg                        >  2:291720/1‑101 (MQ=255)
                                                                             cTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCgg                        >  2:2877766/1‑101 (MQ=255)
                                                                              tGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCggg                       <  1:2893667/101‑1 (MQ=255)
                                                                              tGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCggg                       >  2:2343081/1‑101 (MQ=255)
                                                                                cTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTc                     <  1:2124313/101‑1 (MQ=255)
                                                                                 tttC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCg                    >  2:3767707/1‑101 (MQ=255)
                                                                                    c.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTc                 <  2:3360183/101‑1 (MQ=255)
                                                                                        gggACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTcccc              <  2:4473153/101‑1 (MQ=255)
                                                                                        gggACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTcccc              >  2:542128/1‑101 (MQ=255)
                                                                                        gggACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGCCGCTTACGACGCCTATCTGATTTCGGGTCGCTcccc              >  2:4707502/1‑101 (MQ=255)
                                                                                         ggACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTccccc             <  2:21982/101‑1 (MQ=255)
                                                                                             tGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCcgcg         >  2:3889131/1‑101 (MQ=255)
                                                                                             tGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCcgcg         >  2:4558728/1‑101 (MQ=255)
                                                                                             tGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCcgcg         >  2:782303/1‑101 (MQ=255)
                                                                                              ggCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCgg        >  2:1025359/1‑101 (MQ=255)
                                                                                              ggCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCgg        >  1:3158435/1‑101 (MQ=255)
                                                                                               gCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCggg       <  2:5014946/101‑1 (MQ=255)
                                                                                               gCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCggg       <  1:2653785/101‑1 (MQ=255)
                                                                                                 cAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCggggg     >  1:1465314/1‑101 (MQ=255)
                                                                                                     |                                                                                                   
GATCGAAACGGACCTGAAATTCCGCGGGAAACGCGGCGGCGAGCACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCTTTC.GTGGGACTGGCCAGCCGCCCGCGCCAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGATTTCGTCGCTTACGACGCCTATCTGATTTCGGGTCGCTCCCCCCGCGGGGGCGT  >  NC_014034/1396461‑1396660

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG