Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,450,8720A→C23.8% 11.342562yrfDpredicted pilus assembly protein
Reads supporting (aligned to +/- strand):  new base (0/10):  ref base (30/2):  total (30/12)

CGCCAGTTTGCGCCCCCAGCGATAGCGTGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATT  >  REL606/3450837‑3450907
                                   |                                   
cGCCAGTTTGCGCCCCCAGCGATAGCGTGTCGCCca                                     >  2:2317575/1‑36 (MQ=255)
cGCCAGCTTGCGCCCCCAGCGATAGCGTGTCGCCca                                     >  2:363341/1‑36 (MQ=255)
  ccAGTTTGCGCCCCCAGCGATAGCGTGTCGCCcaca                                   <  1:2534026/36‑1 (MQ=255)
  ccAGTTTGCGCCCCCAGCGATAGCGTGTCGCCcaca                                   <  2:1768843/36‑1 (MQ=255)
  ccAGTTTGCGCCCCCAGCGATAGCGTGTCGCCcaca                                   <  1:3725662/36‑1 (MQ=255)
    aGTTTGGGCCCCCAGCGATAGCGTGTCGCCCACAGc                                 <  1:1765420/36‑1 (MQ=255)
    aGTTTGCGCCCCCAGCGATAGCGTGTCGCCCACAGc                                 >  1:933875/1‑36 (MQ=255)
     gTTTGTGCCCCCAGCGATAGCGTGTCGCTCACAGcc                                >  2:2963235/1‑36 (MQ=37)
      tttGCGCCCCCAGCGATAGCGTGTCGCCCACAGCCa                               <  2:2788947/36‑1 (MQ=255)
      tttGCGCCCCCACCGACAGCGCGTCGCCCACAACCa                               >  1:3274012/1‑36 (MQ=19)
        tGCGCCCCCAGCGATAGCGTGTCGCCCACAGCCACt                             >  2:1811151/1‑36 (MQ=255)
         gcgcCCCCAGCGATAGCGTGTCGCCCACAGCCACTg                            >  1:3518304/1‑36 (MQ=255)
          cccccgcAGCAATAGCGTCCCCCCCCCAGCCACTGt                           <  2:749610/34‑1 (MQ=2)
           gcCCCCAGCGATAGCGTGTCGCCCACAGCCACTGtt                          >  2:3269324/1‑36 (MQ=255)
             ccccAGCAATAGCTTCTCCCCCCCAGCCACTGTTCg                        <  1:1151208/36‑1 (MQ=12)
             ccacAGAGATAGCGTGCCCCCCCCAGCCACTGTTCg                        <  1:945458/33‑1 (MQ=16)
              cccAGCGAAGGCGTCTCCCCCCCAGCCACTGTTCGt                       <  1:3047353/36‑1 (MQ=16)
                cAGCTATAGCGTCCCCCCCCCAGCCACTGTTCGTTa                     <  1:1295169/36‑1 (MQ=16)
                cAGCGATAGCGTGTCGCCCACAGCCACTGTTCGCTa                     >  1:3112227/1‑36 (MQ=255)
                cAGAGATAGCGTGTCGCACACAGCCACTGTTCGTTa                     >  2:76883/1‑36 (MQ=38)
                 aGCGATAGCGTGTCGCCCACAGCCACTGTTCGTTAt                    >  2:2550618/1‑36 (MQ=255)
                  ataaTAGCGTCCCCCCCCCAGCCACTGTTCGTTATc                   <  1:1238161/33‑1 (MQ=16)
                    gATAGCGTGTCGCCCACAGCCACTGTTCGTTTTctc                 >  2:3408630/1‑34 (MQ=38)
                     aTAGCGTGTCGCCCACAGCCACTGTTCGTTATCACg                >  2:2791084/1‑36 (MQ=255)
                      tAGCGTGTCGCCCACAGCCACTGTTCGTTATCACgc               >  1:3023876/1‑36 (MQ=255)
                       aGCGTGTCGCCCACAGCCACTGTTCGTTATCACgcc              >  2:2503410/1‑36 (MQ=255)
                       aGCGTGTCGCCCACACCCACTGTTCGTTATCACgcc              >  2:3142592/1‑36 (MQ=255)
                        gCGTGTCGCCCACAGCCATTGTTCGTTATTACgcca             >  2:803468/1‑36 (MQ=37)
                         cGTGTCGCCCACAGCCACTGTTCGTTATCACgccag            >  2:629426/1‑36 (MQ=255)
                         cGTGTCGCCCACAGCCACTGTTCGTTATCACgccag            >  2:1345905/1‑36 (MQ=255)
                          gtgtCGCCCACAGCCACTGTTCGTTATCACaccagg           >  2:447140/1‑36 (MQ=255)
                           tgtCGCCCACAGCCACTGTTCGTTATCACgccaggc          >  1:2288408/1‑36 (MQ=255)
                           tgtCGCCCACAGCCACTGTTCGTTATCAAgccaggc          >  1:2515561/1‑36 (MQ=255)
                           tgtCGCCCACAGCCACTGTTCCTTATCACgccaggc          >  1:3431967/1‑36 (MQ=255)
                             ctccccCCAGCCACTGTTCGTTATCACgccaggcca        <  1:634140/33‑1 (MQ=37)
                              cGCCCACAGCCACTGTTCGTTATCACgcccggccag       >  2:3351599/1‑36 (MQ=255)
                              cGCCCACAGCCACTGTTCGTTATCACgccaggccag       >  2:1671321/1‑36 (MQ=255)
                              cGCCCACAGCCACTGTTCGTTATCACgccaggccag       >  2:227368/1‑36 (MQ=255)
                               gCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGa      >  2:1970992/1‑36 (MQ=255)
                               gCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGa      >  2:370817/1‑36 (MQ=255)
                               gCCCACAGCCACAGTTCGCTCTCACGCCACCCCAGa      >  2:1684570/1‑36 (MQ=16)
                               ccccCCAGCCACTGTTCGTTATCACGCCAGGCCAGa      <  1:367894/35‑1 (MQ=255)
                                cccCCAGCCACTGTTCGTTATCACGCCAGGCCAGAc     <  1:3333163/36‑1 (MQ=255)
                                cccCCAGCCACTGTTCGTTATCACGCCAGGCCAGAc     <  2:257521/36‑1 (MQ=255)
                                cccACAGCCACTGTTCGTTATCACGCCAGGCCAGAc     >  2:2659598/1‑36 (MQ=255)
                                 ccACAGCCACTGTTCGTTATCACGCCAGGCCAGACa    >  2:707958/1‑36 (MQ=255)
                                 ccACAGCCACTGTTCGTTATCACGCCAGGCCAGACa    >  1:1938221/1‑36 (MQ=255)
                                  cacaGCCACTGTTCGCTATCACGCCAGGCCAGACAt   >  1:77698/1‑36 (MQ=255)
                                   acaGCCACTGTTCGTTATCACGCCAGGCCAGACAtt  >  2:2523351/1‑36 (MQ=255)
                                   |                                   
CGCCAGTTTGCGCCCCCAGCGATAGCGTGTCGCCCACAGCCACTGTTCGTTATCACGCCAGGCCAGACATT  >  REL606/3450837‑3450907

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG