Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606398,4450T→G16.7% 11.548V173G (GTG→GGG) secFprotein export protein SecF
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (7/33):  total (15/33)

AGGTTTCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCACGACGTTATCATTACGC  >  REL606/398410‑398480
                                   |                                   
gggTTTTCGCTTTGAGTGGCGACTGGCGGCAGGggt                                     <  1:686551/35‑1 (MQ=255)
  gTTTCCGCTTTGAGTGGCGACTGGCGGGAGGggggg                                   >  1:5379292/1‑36 (MQ=25)
  gTTTCCGCTTTGAGTGGCGACTGGCGGCTGGggggg                                   >  1:663628/1‑36 (MQ=25)
  gTTTCCGCTTTGAGTGGCGACTGGCGGCAGGggtgg                                   <  1:3505041/36‑1 (MQ=255)
  gTTTCCGCTTTGAGTGGCGACTGGCGGCAGGggtgg                                   <  1:2428956/36‑1 (MQ=255)
    ttCCGCTTTGAGTGGCGACTGGCGGGCGGggggggt                                 >  1:5755095/1‑34 (MQ=18)
    ttCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGtt                                 <  1:1739083/36‑1 (MQ=255)
    ttCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGtt                                 <  1:6785160/36‑1 (MQ=255)
     tCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGTTa                                <  1:7370978/36‑1 (MQ=255)
      ccGCTTTGAGTGGCGACTGGCGGCCGGGGGGGGTAt                               >  1:4011822/1‑36 (MQ=21)
         cTTTGAGTGGCGACTGGCGGCAGGGGTGGTTATTgc                            <  1:7152430/36‑1 (MQ=255)
          tttGAGTGGCGACTGGCGGCAGGGGTGGTTATTgcg                           <  1:3724984/36‑1 (MQ=255)
          tttGAGTGGCGACTGGCGGCAGGGGGGGGTTTTggg                           >  1:5540415/1‑34 (MQ=21)
           ttGAGTGGCGACTGGCGGCAGGGGTGGTTATTgcgc                          <  1:910973/36‑1 (MQ=255)
           ttGAGTGGCGACTGGCGGCAGGGGGGGGTTTTgagt                          >  1:3130999/1‑33 (MQ=18)
            tGAGTGGCGACTGGCGGGAGGGGGGGGGATTGCGCt                         >  1:7113416/1‑36 (MQ=18)
            tGAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCt                         <  1:2418772/36‑1 (MQ=255)
            tGAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCt                         <  1:7472132/36‑1 (MQ=255)
            tGAGTGGCGACTGGCGGCAGGGGGGGTTAATGGGCt                         >  1:6328960/1‑36 (MQ=21)
            ggAGGGGCGACTGGCGGCAGGGGTGGTTATTGCGCt                         <  1:664555/35‑1 (MQ=255)
             gAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTg                        <  1:3311228/36‑1 (MQ=255)
              aGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTgg                       <  1:1834616/36‑1 (MQ=255)
              aGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTgg                       <  1:808025/36‑1 (MQ=255)
               gTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTGgc                      <  1:3737726/36‑1 (MQ=255)
               gTGGCGACTGGCGGCAGGGGGGGGTATTGGGCTggg                      >  1:6387482/1‑35 (MQ=21)
                tGGCGACTGGCGGCAGGGGGGGTTATTCACCTGgcg                     >  1:5378044/1‑36 (MQ=18)
                  gcaaCTGGCGGCAGGGGTGGTTATTGCGCTGGCGCa                   <  1:2423151/33‑1 (MQ=255)
                  gCGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCa                   <  1:2964712/36‑1 (MQ=255)
                   caaCTGGCGGCAGGGGTGGTTATTGCGCTGGCGCac                  <  1:1709791/34‑1 (MQ=255)
                   cGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCac                  <  1:4862824/36‑1 (MQ=255)
                   cGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCac                  <  1:4220432/36‑1 (MQ=255)
                     aCTGGCGGCAGGGGTGGTTATTGCGCTGGCGCacga                <  1:3226082/36‑1 (MQ=255)
                      cTGGCGGCAGGGGTGGTTATTGCGCTGGCGCacgac               <  1:6095714/36‑1 (MQ=255)
                      cTGGCGGCAGGGGTGGTTATTGCGCTGGCGCacgac               <  1:3048379/36‑1 (MQ=255)
                      cTGGCGGCAGGGGTGGTTATTGCGCTGGCGCacgac               <  1:1824306/36‑1 (MQ=255)
                       tGGCGGCAGGGGTGGTTATTGCGCTGGCGCacgacg              <  1:2657839/36‑1 (MQ=255)
                       tGGCGGCAGGGGTGGTTATTGCGCTGGCGCacgacg              <  1:5051158/36‑1 (MQ=255)
                        ggcggcAGGGGTGGTTATTGCGCTGGCGCACGACGt             <  1:1162430/36‑1 (MQ=255)
                        ggcggcAGGGGTGGTTATTGCGCTGGCGCACGACGt             <  1:4574346/36‑1 (MQ=255)
                         gcggcAGGGGTGGTTATTGCGCTGGCGCACGACGtt            <  1:4682039/36‑1 (MQ=255)
                         gcggcAGGGGGGGGTATTGCGCTGGCGCacgacggt            >  1:3927277/1‑34 (MQ=21)
                           ggcAGGGGTGGTTATTTCGCTGGCGCACGGCGtttt          >  1:1889914/1‑34 (MQ=21)
                           ggcAGGGGTGGTTATTGCGCTGGCGCACGACGTTAt          >  1:6179985/1‑36 (MQ=255)
                            gcAGGGGTGGTTATTGCGCTGGCGCACGACGTTATc         <  1:6100071/36‑1 (MQ=255)
                             cAGGGGTGGTTATTGCGCTGGCGCGCGACGTTATCa        >  1:5344772/1‑36 (MQ=255)
                              aGGGGTGGTTATTTCGCTGGCGGACGACGTTATCAt       >  1:867496/1‑36 (MQ=25)
                              aGGGGTGGTTATTGCGCTGGCGCACGACGTTATCAt       <  1:1406869/36‑1 (MQ=255)
                              aGGGGTGGTTATTGCGCTGGCGCACGACGTTATCAt       <  1:1917151/36‑1 (MQ=255)
                               ggggTGGTTATTGCGCTGGCGCACGACGTTATCAtt      >  1:1906208/1‑36 (MQ=255)
                               ggggTGGTTATTGCGCTGGCGCACGACGTTATCAtt      <  1:62740/36‑1 (MQ=255)
                                gggTGGTTATTGCGCTGGCGCACGACGTTATCATTa     >  1:3874052/1‑36 (MQ=255)
                                gggTGGTTATTGCGCTGGCGCACGACGTTATCATTa     <  1:6432684/36‑1 (MQ=255)
                                gggTGGTTATTGCGCTGGCGCACGACGTTATCATTa     >  1:1698234/1‑36 (MQ=255)
                                gggTGGTTATTGCGCTGGCGCACGACGTTATCATTa     <  1:6030614/36‑1 (MQ=255)
                                 ggtggtTATTGCGCTGGCGCCCGACGTTACCAttcc    >  1:140538/1‑34 (MQ=21)
                                 ggtggtTATTGCGCTGGCGCACGACGTTATCATTAc    <  1:3243357/36‑1 (MQ=255)
                                 ggtggtTATTGCGCTGGCGCACGACGTTATCATTAc    <  1:4299472/36‑1 (MQ=255)
                                 ggtggtTATTGCGCTGGCGCACGACGTTATCATTAc    >  1:4618373/1‑36 (MQ=255)
                                  gtggtTATTGCGCTGGCGCACGACGTTATCCTTACg   >  1:6681272/1‑36 (MQ=255)
                                  gtggtTATTGCGCTGGCGCACGACGTTATCATTACg   <  1:2976070/36‑1 (MQ=255)
                                  gtggtTATTGCGCTGGCGCACGACGTTATCATTACg   <  1:1914010/36‑1 (MQ=255)
                                   tggtTATTTCGGTGGCGCACGACGTTATTATTActc  >  1:6554553/1‑34 (MQ=21)
                                   tggtTATTGCGCTGGCGCACGACGTTATCATTACGc  <  1:6935558/36‑1 (MQ=255)
                                   tggtTATTGCGCTGGCGCACGACGTTATCATTACGc  <  1:4035180/36‑1 (MQ=255)
                                   |                                   
AGGTTTCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCACGACGTTATCATTACGC  >  REL606/398410‑398480

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 28 ≤  ATCG  < 32 ≤  ATCG