Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,797,132 2 bp→CG 100% coding (156‑157/2238 nt) RCAP_rcc01660 → signal transduction histidine kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,797,1320G→C100.0% 636.1 172R52R (CGG→CGCRCAP_rcc01660signal transduction histidine kinase
Reads supporting (aligned to +/- strand):  new base (92/80):  ref base (0/0):  total (92/80)
*NC_0140341,797,1330C→G100.0% 654.2 172H53D (CAT→GAT) RCAP_rcc01660signal transduction histidine kinase
Reads supporting (aligned to +/- strand):  new base (92/80):  ref base (0/0):  total (92/80)

Only 100 of 174 total aligned reads displayed.

CGAGAAGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGATACGCGGGGCCGCCTGGGTGCAGTTCGAACCGGCCGGCATGCTGGAGG.ACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTGGCCGAGG  >  NC_014034/1797036‑1797220
                                                                                                ||                                                                                        
cGAGAAGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGc                                                                                      >  2:2507255/1‑101 (MQ=255)
cGAGAAGCTTCTGACCATCGTCGGCGCGCTGATCCGACATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGc                                                                                      <  1:1121970/101‑1 (MQ=255)
 gagaAGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCt                                                                                     >  1:916216/1‑101 (MQ=255)
  agaAGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTg                                                                                    >  1:3249386/1‑101 (MQ=255)
  agaAGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTg                                                                                    >  1:3136575/1‑101 (MQ=255)
   gaAGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTgg                                                                                   >  1:1008727/1‑101 (MQ=255)
   gaAGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTgg                                                                                   >  2:1709942/1‑101 (MQ=255)
    aaGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTgga                                                                                  <  1:2097632/101‑1 (MQ=255)
      gCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg                                                                                <  1:24466/101‑1 (MQ=255)
      gCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg                                                                                >  1:349809/1‑101 (MQ=255)
       cTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.a                                                                              >  1:2637022/1‑101 (MQ=255)
       cTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.a                                                                              <  1:3702501/101‑1 (MQ=255)
       cTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.a                                                                              >  2:13394/1‑101 (MQ=255)
       cTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.a                                                                              >  2:2092383/1‑101 (MQ=255)
       cTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGCGCAGTTGGAACGCGCCGCGATGCTggagg.a                                                                              >  1:2949396/1‑101 (MQ=255)
         tCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.aac                                                                            >  1:4021731/1‑99 (MQ=255)
         tCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.aac                                                                            >  1:341418/1‑99 (MQ=255)
          cTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.aaca                                                                           <  1:2692679/101‑4 (MQ=255)
           tGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.aacag                                                                          >  1:1906435/1‑97 (MQ=255)
            gACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTggagg.aacagg                                                                         <  1:2364819/101‑6 (MQ=255)
             gccATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGt                                                                         <  2:666023/100‑1 (MQ=255)
              ccATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTc                                                                        >  1:1081827/1‑101 (MQ=255)
              ccATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTc                                                                        >  2:2796011/1‑101 (MQ=255)
               cATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTcc                                                                       >  2:2461895/1‑101 (MQ=255)
               cATCGGCGGGGCGCTGATCCGCCATGCCGGAGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTcc                                                                       <  1:2018737/101‑1 (MQ=255)
                 tcgtcgGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCcgc                                                                     <  2:349758/101‑1 (MQ=255)
                      ggCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCg                                                                <  2:3136344/101‑1 (MQ=255)
                         gcgcTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCAc                                                             <  2:143013/101‑1 (MQ=255)
                         gcgcTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCAc                                                             <  1:12494/101‑1 (MQ=255)
                          cgcTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCAcc                                                            >  1:3561486/1‑101 (MQ=255)
                           gcTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCAccc                                                           <  1:501566/101‑1 (MQ=255)
                            cTGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCcg                                                          >  1:1413158/1‑101 (MQ=255)
                             tGATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCcgc                                                         >  1:2252441/1‑101 (MQ=255)
                              gATCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCcgcg                                                        >  1:1893335/1‑101 (MQ=255)
                                tCCGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGAt                                                      <  2:3787198/101‑1 (MQ=255)
                                 ccGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATc                                                     >  2:549887/1‑101 (MQ=255)
                                  cGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCt                                                    <  2:2754589/101‑1 (MQ=255)
                                  cGCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCt                                                    >  2:3877827/1‑101 (MQ=255)
                                   gCCATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTg                                                   <  1:1005657/101‑1 (MQ=255)
                                    ccATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTgg                                                  >  2:2570336/1‑101 (MQ=255)
                                     cATGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGACCTGGa                                                 >  1:1620335/1‑101 (MQ=255)
                                      aTGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGaa                                                <  2:1791379/101‑1 (MQ=255)
                                      aTGCCGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGaa                                                >  1:2732963/1‑101 (MQ=255)
                                        gCCGAGGAAACCCCAGAGGAT.CGCGGCGCCGACTGGGTGCAGTTCGAACGCGACGAGATACTGGAGGAACAGGTACGAGAACGCACCCGCGATCTGGAACg                                              >  1:3461424/1‑101 (MQ=255)
                                         ccGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGTACAGGTCCGCGACCGCACCCGCGGTCTGAGACgg                                             >  2:4082950/1‑101 (MQ=255)
                                         ccGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACgg                                             <  1:1995911/101‑1 (MQ=255)
                                          cGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACggg                                            <  2:108011/101‑1 (MQ=255)
                                          cGAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACggg                                            >  1:894657/1‑101 (MQ=255)
                                           gAGGAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGgc                                           >  2:1247639/1‑101 (MQ=255)
                                             ggAAACCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGgcgc                                         >  2:3431485/1‑101 (MQ=255)
                                               aaaCCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGGCCGCACCCGCCACCTGGACCGGGCGCtg                                       >  2:593737/1‑100 (MQ=255)
                                               aaaCCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTc                                       <  2:1893215/101‑1 (MQ=255)
                                               aaaCCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTc                                       >  1:926271/1‑101 (MQ=255)
                                               aaaCCCCCGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTc                                       <  2:4041681/101‑1 (MQ=255)
                                                  cccccGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGAt                                    >  1:2569153/1‑101 (MQ=255)
                                                   ccccGAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCCCACCCCCGATCTGGAACGGGCGCTCCATc                                   >  2:2952917/1‑101 (MQ=255)
                                                       gAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCt                               <  2:2873275/101‑1 (MQ=255)
                                                       gAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCt                               >  1:2793847/1‑101 (MQ=255)
                                                       gAGGAT.CGCGGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCt                               >  2:4189502/1‑101 (MQ=255)
                                                          ggatcgcgGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCa                             >  2:2378577/5‑101 (MQ=255)
                                                           gatcgcgGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCaa                            <  1:216219/98‑1 (MQ=255)
                                                            gaggcgGTGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAAt                           >  2:2740190/4‑101 (MQ=255)
                                                            atcgcgGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAAt                           <  2:1645977/99‑1 (MQ=255)
                                                            atcgcgGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAAt                           <  1:4108637/99‑1 (MQ=255)
                                                              cgcgGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGc                         >  1:818004/1‑101 (MQ=255)
                                                               gcgGCGCCGCCTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGcc                        <  1:2092521/101‑1 (MQ=255)
                                                                 ggcgccgccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTc                      >  2:4010312/4‑101 (MQ=255)
                                                                 ggcgccgccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTc                      >  1:2622601/4‑101 (MQ=255)
                                                                 ggcgccgccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTc                      <  2:3168426/98‑1 (MQ=255)
                                                                   cgccgccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGa                    <  2:2272760/100‑1 (MQ=255)
                                                                    gccgccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGaa                   <  1:3878086/101‑1 (MQ=255)
                                                                     ccgccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAAt                  <  2:2743151/101‑1 (MQ=255)
                                                                      cgccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATg                 <  1:3239678/101‑1 (MQ=255)
                                                                       gccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATgc                <  2:268965/101‑1 (MQ=255)
                                                                       gccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATgc                <  1:2507439/101‑1 (MQ=255)
                                                                       gccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATgc                <  1:547472/101‑1 (MQ=255)
                                                                       gccTGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGACTCGAATgc                >  1:338895/1‑101 (MQ=255)
                                                                          tGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGcgcg             >  2:112317/1‑101 (MQ=255)
                                                                          tGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGcgcg             >  2:408289/1‑101 (MQ=255)
                                                                          tGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGcgcg             >  1:1309272/1‑101 (MQ=255)
                                                                          tGGGTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGcgcg             >  2:802316/1‑101 (MQ=255)
                                                                           gggTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCgg            >  1:431943/1‑101 (MQ=255)
                                                                           gggTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCgg            <  2:3946798/101‑1 (MQ=255)
                                                                           gggTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCgg            >  2:346842/1‑101 (MQ=255)
                                                                           gggTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCgg            <  1:1061564/101‑1 (MQ=255)
                                                                            ggTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGc           <  1:2878599/101‑1 (MQ=255)
                                                                            ggTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGc           <  2:3907432/101‑1 (MQ=255)
                                                                             gTGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCt          <  1:13404/101‑1 (MQ=255)
                                                                              tGCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTg         <  2:30872/101‑1 (MQ=255)
                                                                               gCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGTGGCTgg        >  1:686928/1‑101 (MQ=255)
                                                                               gCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTgg        <  2:3249145/101‑1 (MQ=255)
                                                                               gCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTgg        <  1:4214054/101‑1 (MQ=255)
                                                                               gCAGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTgg        >  1:2856719/1‑101 (MQ=255)
                                                                                 aGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCTCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAAGGCGCGGCTGGcc      >  1:712532/1‑101 (MQ=255)
                                                                                 aGTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTGGcc      >  1:102282/1‑101 (MQ=255)
                                                                                  gTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTGGCCg     <  2:1667678/101‑1 (MQ=255)
                                                                                  gTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTGGCCg     >  1:3506051/1‑101 (MQ=255)
                                                                                  gTTCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTGGCCg     <  1:2740396/101‑1 (MQ=255)
                                                                                    tCGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTGGCCGAg   >  2:939094/1‑101 (MQ=255)
                                                                                     cGAACGCGCCGCGATGCTGGAGGAACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTGGCCGAgg  <  1:3867128/101‑1 (MQ=255)
                                                                                                ||                                                                                        
CGAGAAGCTTCTGACCATCGTCGGCGCGCTGATCCGCCATGCCGAGGAAACCCCCGAGGATACGCGGGGCCGCCTGGGTGCAGTTCGAACCGGCCGGCATGCTGGAGG.ACAGGTCCGCGACCGCACCCGCGATCTGGAACGGGCGCTCGATCTGCTCAATGCCTCGAATGCGCGGCTGGCCGAGG  >  NC_014034/1797036‑1797220

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG