Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606398,4450T→G16.7% 10.148518secFprotein 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                                   >  2:1640082/1‑36 (MQ=37)
  gTTTCCGCTTTGAGTGGCGACTGGCGGCTGGggggg                                   >  1:663628/1‑36 (MQ=37)
  gTTTCCGCTTTGAGTGGCGACTGGCGGCAGGggtgg                                   <  1:3505041/36‑1 (MQ=255)
  gTTTCCGCTTTGAGTGGCGACTGGCGGCAGGggtgg                                   <  1:2428956/36‑1 (MQ=255)
    ttCCGCTTTGAGTGGCGACTGGCGGGCGGggggggt                                 >  2:2015885/1‑34 (MQ=21)
    ttCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGtt                                 <  1:1739083/36‑1 (MQ=255)
    ttCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGtt                                 <  2:3045950/36‑1 (MQ=255)
     tCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGTTa                                <  2:3631768/36‑1 (MQ=255)
      ccGCTTTGAGTGGCGACTGGCGGCCGGGGGGGGTAt                               >  2:272612/1‑36 (MQ=25)
         cTTTGAGTGGCGACTGGCGGCAGGGGTGGTTATTgc                            <  2:3413220/36‑1 (MQ=255)
          tttGAGTGGCGACTGGCGGCAGGGGTGGTTATTgcg                           <  1:3724984/36‑1 (MQ=255)
          tttGAGTGGCGACTGGCGGCAGGGGGGGGTTTTggg                           >  2:1801205/1‑34 (MQ=21)
           ttGAGTGGCGACTGGCGGCAGGGGTGGTTATTgcgc                          <  1:910973/36‑1 (MQ=255)
           ttGAGTGGCGACTGGCGGCAGGGGGGGGTTTTgagt                          >  1:3130999/1‑33 (MQ=19)
            tGAGTGGCGACTGGCGGGAGGGGGGGGGATTGCGCt                         >  2:3374206/1‑36 (MQ=19)
            tGAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCt                         <  1:2418772/36‑1 (MQ=255)
            tGAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCt                         <  2:3732922/36‑1 (MQ=255)
            tGAGTGGCGACTGGCGGCAGGGGGGGTTAATGGGCt                         >  2:2589750/1‑36 (MQ=25)
            ggAGGGGCGACTGGCGGCAGGGGTGGTTATTGCGCt                         <  1:664555/35‑1 (MQ=255)
             gAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTg                        <  1:3311228/36‑1 (MQ=255)
              aGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTgg                       <  1:808025/36‑1 (MQ=255)
              aGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTgg                       <  1:1834616/36‑1 (MQ=255)
               gTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTGgc                      <  1:3737726/36‑1 (MQ=255)
               gTGGCGACTGGCGGCAGGGGGGGGTATTGGGCTggg                      >  2:2648272/1‑35 (MQ=25)
                tGGCGACTGGCGGCAGGGGGGGTTATTCACCTGgcg                     >  2:1638834/1‑36 (MQ=19)
                  gcaaCTGGCGGCAGGGGTGGTTATTGCGCTGGCGCa                   <  1:2423151/33‑1 (MQ=255)
                  gCGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCa                   <  1:2964712/36‑1 (MQ=255)
                   caaCTGGCGGCAGGGGTGGTTATTGCGCTGGCGCac                  <  1:1709791/34‑1 (MQ=255)
                   cGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCac                  <  2:1123614/36‑1 (MQ=255)
                   cGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCac                  <  2:481222/36‑1 (MQ=255)
                     aCTGGCGGCAGGGGTGGTTATTGCGCTGGCGCacga                <  1:3226082/36‑1 (MQ=255)
                      cTGGCGGCAGGGGTGGTTATTGCGCTGGCGCacgac               <  2:2356504/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              <  2:1311948/36‑1 (MQ=255)
                        ggcggcAGGGGTGGTTATTGCGCTGGCGCACGACGt             <  2:835136/36‑1 (MQ=255)
                        ggcggcAGGGGTGGTTATTGCGCTGGCGCACGACGt             <  1:1162430/36‑1 (MQ=255)
                         gcggcAGGGGTGGTTATTGCGCTGGCGCACGACGtt            <  2:942829/36‑1 (MQ=255)
                         gcggcAGGGGGGGGTATTGCGCTGGCGCacgacggt            >  2:188067/1‑34 (MQ=25)
                           ggcAGGGGTGGTTATTTCGCTGGCGCACGGCGtttt          >  1:1889914/1‑34 (MQ=25)
                           ggcAGGGGTGGTTATTGCGCTGGCGCACGACGTTAt          >  2:2440775/1‑36 (MQ=255)
                            gcAGGGGTGGTTATTGCGCTGGCGCACGACGTTATc         <  2:2360861/36‑1 (MQ=255)
                             cAGGGGTGGTTATTGCGCTGGCGCGCGACGTTATCa        >  2:1605562/1‑36 (MQ=255)
                              aGGGGTGGTTATTTCGCTGGCGGACGACGTTATCAt       >  1:867496/1‑36 (MQ=37)
                              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     <  2:2693474/36‑1 (MQ=255)
                                gggTGGTTATTGCGCTGGCGCACGACGTTATCATTa     >  1:1698234/1‑36 (MQ=255)
                                gggTGGTTATTGCGCTGGCGCACGACGTTATCATTa     <  2:2291404/36‑1 (MQ=255)
                                gggTGGTTATTGCGCTGGCGCACGACGTTATCATTa     >  2:134842/1‑36 (MQ=255)
                                 ggtggtTATTGCGCTGGCGCCCGACGTTACCAttcc    >  1:140538/1‑34 (MQ=25)
                                 ggtggtTATTGCGCTGGCGCACGACGTTATCATTAc    <  1:3243357/36‑1 (MQ=255)
                                 ggtggtTATTGCGCTGGCGCACGACGTTATCATTAc    <  2:560262/36‑1 (MQ=255)
                                 ggtggtTATTGCGCTGGCGCACGACGTTATCATTAc    >  2:879163/1‑36 (MQ=255)
                                  gtggtTATTGCGCTGGCGCACGACGTTATCCTTACg   >  2:2942062/1‑36 (MQ=255)
                                  gtggtTATTGCGCTGGCGCACGACGTTATCATTACg   <  1:1914010/36‑1 (MQ=255)
                                  gtggtTATTGCGCTGGCGCACGACGTTATCATTACg   <  1:2976070/36‑1 (MQ=255)
                                   tggtTATTTCGGTGGCGCACGACGTTATTATTActc  >  2:2815343/1‑34 (MQ=21)
                                   tggtTATTGCGCTGGCGCACGACGTTATCATTACGc  <  2:295970/36‑1 (MQ=255)
                                   tggtTATTGCGCTGGCGCACGACGTTATCATTACGc  <  2:3196348/36‑1 (MQ=255)
                                   |                                   
AGGTTTCCGCTTTGAGTGGCGACTGGCGGCAGGGGTGGTTATTGCGCTGGCGCACGACGTTATCATTACGC  >  REL606/398410‑398480

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