Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,424,8630A→T27.9% 13.244intergenic (+81/‑45)nirD/nirCnitrite reductase small subunit/nitrite transporter
Reads supporting (aligned to +/- strand):  new base (0/12):  ref base (8/23):  total (8/36)

TAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTG.TTT.TATT.TTTTGTTTTTATTTTTTAAA  >  REL606/3424828‑3424897
                                   |                                     
taataaTCTTCGGTATATTGCAGGACATTTTTTaaa                                       <  1:1111134/36‑1 (MQ=255)
 aataatCTTCGGTATATTGCAGGACATTTTTTAAAc                                      <  1:1953553/36‑1 (MQ=255)
  ataatCTTCGGTATATTGCAGGACATTTTTTAAACt                                     >  1:1398758/1‑36 (MQ=255)
  ataatCTTCGGTATATTGCAGGACATTTTTTAAACt                                     <  1:1878744/36‑1 (MQ=255)
  ataatCTTCGGTATATTGCAGGACATTTTTTAAACt                                     >  2:559017/1‑36 (MQ=255)
   taatCTTCGGTATATTGCAGGACATTTTTTAAACtt                                    <  2:1956697/36‑1 (MQ=255)
    aatCTTCGGTATATTGCAGGACATTTTTTAAACttt                                   >  2:3379751/1‑36 (MQ=255)
    aatCTTCGGTATATTGCAGGACATTTTTTAAACttt                                   <  2:662677/36‑1 (MQ=255)
      tCTTCGGTATATTGCAGGACATTTTTTAAACttttt                                 <  1:1029447/36‑1 (MQ=255)
      tCTTCGGTATATTGCAGGACATTTTTTAAACtttat                                 >  1:429933/1‑34 (MQ=255)
       cTTCGGTATATTGCAGGACATTTTTTAAACTTTTTg                                <  1:3135749/36‑1 (MQ=255)
           ggTATATTGCAGGACATTTTTTAAACTTTTTG.ttt.t                          <  2:538291/36‑1 (MQ=255)
             tataTTGCAGGACATTTTTTAAACTTTTTG.TTT.TAt                        <  1:3247150/36‑1 (MQ=255)
             tataTTGCAGGACATTTTTTAAACTTTTTG.TTT.TAt                        <  1:1236580/36‑1 (MQ=255)
              atatTGCAGGACATTTTTTAAACTTTTTG.TTT.TAtt                       >  1:1183741/1‑36 (MQ=255)
              atatTGCAGGACATTTTTTAAACTTTTTG.TTT.TAtt                       <  1:1282723/36‑1 (MQ=255)
               tatTGCAGGACATTTTTTAAACTTTTTG.TTT.TAtt.t                     <  2:1067867/36‑1 (MQ=255)
               tatTGCAGGACATTTTTTAAACTTTTTG.TTT.TAtt.t                     <  2:1067011/36‑1 (MQ=255)
                 ttGCAGGACATTTTTTAAACTTTTTG.TTT.TAtt.ttt                   >  2:3535968/1‑36 (MQ=255)
                  tGCAGGACATTTTTTAAACTTTTTG.TTT.TAtt.tttt                  <  1:2046500/36‑1 (MQ=255)
                        tttttttttTATGTTTTTG.TTT.TTCT.TTTTTttttt            <  2:2870420/33‑1 (MQ=9)
                        ttcttttttACTTTTTTTG.TTT.TTTT.TTTTTttttt            <  1:3016463/33‑1 (MQ=2)
                        ttaTTTTTTCTTCTTTTTT.TTT.TTTT.TGTTTttttt            <  2:3281751/34‑1 (MQ=0)
                        tctttttttTTACTTTATT.TTT.TGTT.TTTTTttttt            <  2:3366227/33‑1 (MQ=12)
                        aCATTTTTTAAACTTTTTG.TTT.TATT.TTTTGttttt            >  1:3573144/1‑36 (MQ=255)
                         ttttttttATTTTTTTTG.TTT.TTTTGTTTTGttttt            <  2:3468814/34‑1 (MQ=0)
                         tcttttttTCATTTTTTG.TTT.TTTT.TTTTTtttttt           <  1:3543043/34‑2 (MQ=2)
                         cATTTTTTAAACTTTTTG.TTT.TATT.TTTTGTTTTTa           <  2:61866/36‑1 (MQ=255)
                         aaTTTTTTTTTCTTTTTT.TCT.TTTT.TTTTGttattt           <  2:3553285/35‑2 (MQ=2)
                          tttttttAATTTTTGGG.TTT.TATT.TTGTTTTCTTAt          <  2:2709673/35‑1 (MQ=0)
                          gttttttATTTTTTTTG.TTT.TTTT.ATTTTttttttt          <  2:3734748/35‑3 (MQ=0)
                          gttttttATTCTTTATT.TTT.TTTT.TTTTTttttttt          <  2:3671910/35‑3 (MQ=2)
                          cttttttTTACTTTTTT.TTT.TGTT.TTTTTctttttt          <  2:2781403/35‑3 (MQ=2)
                          cttttttATTTTTTCTT.TTT.TTTT.TTTTGttttttt          <  2:3644893/35‑3 (MQ=9)
                          cttttttATTCTTTTTG.TTT.TCTT.TTTTTttttttt          <  1:3198821/35‑3 (MQ=12)
                           ttttttTATTTTTTTG.TTTGTTTT.TTTTGttttttt          <  2:3742731/36‑3 (MQ=0)
                           ttttttATTTTTTTTT.TTT.TGTT.TTTTTttttttgg         <  2:2901300/36‑4 (MQ=2)
                            tttttATATCTTTTT.TTT.TATT.TTTTTTTTTTTttt        <  1:3112836/36‑1 (MQ=9)
                            tttttAAACTTTTTG.TTT.TATT.TTTTGTTTTTAttt        >  2:807886/1‑36 (MQ=255)
                            tttttAAACTTTTTG.TTT.TATT.TTTTGTTTTTAttt        <  1:2516067/36‑1 (MQ=255)
                             ttttATGCTTTTTT.TTT.TGTT.TTTTCTTTTCTtttt       <  2:2751020/36‑1 (MQ=0)
                             ttttATATTTTTTT.TGT.TATT.TTATTTTTTTTtttt       <  2:3281590/36‑1 (MQ=2)
                             ttttAAACTTTTTG.TTT.TATT.TTTTGTTTTTAtttt       <  1:1804578/36‑1 (MQ=255)
                             ttttAAACTTTTTG.TTT.TATT.TTTTGTTTTTAtttt       <  2:3006371/36‑1 (MQ=255)
                               ttAAACTTTTTG.TTT.TATT.TTTTGTTTTTAtttttt     <  2:829907/36‑1 (MQ=255)
                                 aaaCTTTTTG.TTT.TATT.TTTTGTTTTTATTTTTTaa   <  2:3640737/36‑1 (MQ=255)
                                 aaaCTTTTTG.TTT.TATT.TTTTGTTTTTATTTTTTaa   <  2:3298524/36‑1 (MQ=255)
                                  taCTGTTTGTTTT.TATT.TTTTTTTATTTtttttttt   <  2:3385489/35‑3 (MQ=0)
                                  aaCTTTTTG.TTT.TATT.ATTTGTTTTTATTTTTTaaa  <  1:2694106/36‑1 (MQ=255)
                                   aCTTTTTCTTTT.TTTT.TTTTGTTTTTTttttttttt  <  2:3289803/36‑4 (MQ=2)
                                   |                                     
TAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTG.TTT.TATT.TTTTGTTTTTATTTTTTAAA  >  REL606/3424828‑3424897

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