Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,424,8620A→T30.9% 20.744intergenic (+80/‑46)nirD/nirCnitrite reductase small subunit/nitrite transporter
Reads supporting (aligned to +/- strand):  new base (0/13):  ref base (8/21):  total (8/36)

GTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTT.TATT.TTTTGTTTTTATTTTTTAAA  >  REL606/3424827‑3424897
                                   |                                     
gTAATAATCTTCGGTATATTGCAGGACATTTTTTaa                                       <  2:902961/36‑1 (MQ=255)
 taataaTCTTCGGTATATTGCAGGACATTTTTTaaa                                      <  1:1111134/36‑1 (MQ=255)
  aataatCTTCGGTATATTGCAGGACATTTTTTAAAc                                     <  1:1953553/36‑1 (MQ=255)
   ataatCTTCGGTATATTGCAGGACATTTTTTAAACt                                    >  2:559017/1‑36 (MQ=255)
   ataatCTTCGGTATATTGCAGGACATTTTTTAAACt                                    >  1:1398758/1‑36 (MQ=255)
   ataatCTTCGGTATATTGCAGGACATTTTTTAAACt                                    <  1:1878744/36‑1 (MQ=255)
    taatCTTCGGTATATTGCAGGACATTTTTTAAACtt                                   <  2:1956697/36‑1 (MQ=255)
     aatCTTCGGTATATTGCAGGACATTTTTTAAACttt                                  <  2:662677/36‑1 (MQ=255)
     aatCTTCGGTATATTGCAGGACATTTTTTAAACttt                                  >  2:3379751/1‑36 (MQ=255)
       tCTTCGGTATATTGCAGGACATTTTTTAAACttttt                                <  1:1029447/36‑1 (MQ=255)
       tCTTCGGTATATTGCAGGACATTTTTTAAACtttat                                >  1:429933/1‑34 (MQ=255)
        cTTCGGTATATTGCAGGACATTTTTTAAACTTTTTg                               <  1:3135749/36‑1 (MQ=255)
            ggTATATTGCAGGACATTTTTTAAACTTTTTGttt.t                          <  2:538291/36‑1 (MQ=255)
              tataTTGCAGGACATTTTTTAAACTTTTTGTTT.TAt                        <  1:1236580/36‑1 (MQ=255)
              tataTTGCAGGACATTTTTTAAACTTTTTGTTT.TAt                        <  1:3247150/36‑1 (MQ=255)
               atatTGCAGGACATTTTTTAAACTTTTTGTTT.TAtt                       >  1:1183741/1‑36 (MQ=255)
               atatTGCAGGACATTTTTTAAACTTTTTGTTT.TAtt                       <  1:1282723/36‑1 (MQ=255)
                tatTGCAGGACATTTTTTAAACTTTTTGTTT.TAtt.t                     <  2:1067867/36‑1 (MQ=255)
                tatTGCAGGACATTTTTTAAACTTTTTGTTT.TAtt.t                     <  2:1067011/36‑1 (MQ=255)
                  ttGCAGGACATTTTTTAAACTTTTTGTTT.TAtt.ttt                   >  2:3535968/1‑36 (MQ=255)
                   tGCAGGACATTTTTTAAACTTTTTGTTT.TAtt.tttt                  <  1:2046500/36‑1 (MQ=255)
                         tttttttttTATGTTTTTGTTT.TTCT.TTTTTttttt            <  2:2870420/33‑1 (MQ=9)
                         ttcttttttACTTTTTTTGTTT.TTTT.TTTTTttttt            <  1:3016463/33‑1 (MQ=2)
                         ttaTTTTTTCTTCTTTTTTTTT.TTTT.TGTTTttttt            <  2:3281751/34‑1 (MQ=0)
                         tctttttttTTACTTTATTTTT.TGTT.TTTTTttttt            <  2:3366227/33‑1 (MQ=12)
                         aCATTTTTTAAACTTTTTGTTT.TATT.TTTTGttttt            >  1:3573144/1‑36 (MQ=255)
                          ttttttttATTTTTTTTGTTT.TTTTGTTTTGttttt            <  2:3468814/34‑1 (MQ=0)
                          tcttttttTCATTTTTTGTTT.TTTT.TTTTTtttttt           <  1:3543043/34‑2 (MQ=2)
                          cATTTTTTAAACTTTTTGTTT.TATT.TTTTGTTTTTa           <  2:61866/36‑1 (MQ=255)
                          aaTTTTTTTTTCTTTTTTTCT.TTTT.TTTTGttattt           <  2:3553285/35‑2 (MQ=2)
                           tttttttAATTTTTGGGTTT.TATT.TTGTTTTCTTAt          <  2:2709673/35‑1 (MQ=0)
                           gttttttATTTTTTTTGTTT.TTTT.ATTTTttttttt          <  2:3734748/35‑3 (MQ=0)
                           gttttttATTCTTTATTTTT.TTTT.TTTTTttttttt          <  2:3671910/35‑3 (MQ=2)
                           cttttttTTACTTTTTTTTT.TGTT.TTTTTctttttt          <  2:2781403/35‑3 (MQ=2)
                           cttttttATTTTTTCTTTTT.TTTT.TTTTGttttttt          <  2:3644893/35‑3 (MQ=9)
                           cttttttATTCTTTTTGTTT.TCTT.TTTTTttttttt          <  1:3198821/35‑3 (MQ=12)
                            ttttttTATTTTTTTGTTTGTTTT.TTTTGttttttt          <  2:3742731/36‑3 (MQ=0)
                            ttttttATTTTTTTTTTTT.TGTT.TTTTTttttttgg         <  2:2901300/36‑4 (MQ=2)
                             tttttATATCTTTTTTTT.TATT.TTTTTTTTTTTttt        <  1:3112836/36‑1 (MQ=9)
                             tttttAAACTTTTTGTTT.TATT.TTTTGTTTTTAttt        <  1:2516067/36‑1 (MQ=255)
                             tttttAAACTTTTTGTTT.TATT.TTTTGTTTTTAttt        >  2:807886/1‑36 (MQ=255)
                              ttttATGCTTTTTTTTT.TGTT.TTTTCTTTTCTtttt       <  2:2751020/36‑1 (MQ=0)
                              ttttATATTTTTTTTGT.TATT.TTATTTTTTTTtttt       <  2:3281590/36‑1 (MQ=2)
                              ttttAAACTTTTTGTTT.TATT.TTTTGTTTTTAtttt       <  1:1804578/36‑1 (MQ=255)
                              ttttAAACTTTTTGTTT.TATT.TTTTGTTTTTAtttt       <  2:3006371/36‑1 (MQ=255)
                                ttAAACTTTTTGTTT.TATT.TTTTGTTTTTAtttttt     <  2:829907/36‑1 (MQ=255)
                                  aaaCTTTTTGTTT.TATT.TTTTGTTTTTATTTTTTaa   <  2:3298524/36‑1 (MQ=255)
                                  aaaCTTTTTGTTT.TATT.TTTTGTTTTTATTTTTTaa   <  2:3640737/36‑1 (MQ=255)
                                   aaCTTTTTGTTT.TATT.ATTTGTTTTTATTTTTTaaa  <  1:2694106/36‑1 (MQ=255)
                                   |                                     
GTAATAATCTTCGGTATATTGCAGGACATTTTTTAAACTTTTTGTTT.TATT.TTTTGTTTTTATTTTTTAAA  >  REL606/3424827‑3424897

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