Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,524,719 T→A 7.7% S37T (TCC→ACC)  fdnH → formate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,524,7190T→A 7.7% 3.2 ‑1.4 ‑0.739109fdnHformate dehydrogenase‑N, Fe‑S (beta) subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (18/18):  total (21/18)
Fisher's exact test strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.14e-01

ATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGC  >  REL606/1524684‑1524753
                                   |                                  
aTTACAAAGCAGAAGTCGCAAAACTTATCGACGttt                                    >  1:1255434/1‑36 (MQ=255)
aTTACAAAGCAGAAGTCGCAAAACTTATCGACGttt                                    >  1:5909200/1‑36 (MQ=255)
aTTACAAAGCAGAAGACGAAAAACTTATCGACGttt                                    >  1:4886220/1‑36 (MQ=37)
  tACAAAGCAGAAGTCGCCAACCTTATCGACGTTTcc                                  <  1:2049748/36‑1 (MQ=37)
  tACAAAGCAGAAGTCGCAAAACTTATCGACGTTTcc                                  >  1:142274/1‑36 (MQ=255)
  tACAAAGCAGAAGTCGCAAAACTTATCGACGTTTcc                                  <  1:1871251/36‑1 (MQ=255)
  tACAAAGCAGAAGTCGCAAAACTTATCGACGTTTcc                                  >  1:2932206/1‑36 (MQ=255)
    aaaaGCAGAAGTCGCAAAACTTATCGACGTTTCCAc                                <  1:3667840/35‑1 (MQ=255)
     aaaGCAGAAGTCGCAAAACTTATCGACGTTTCCAcc                               >  1:1085443/1‑36 (MQ=255)
     aaaGCAGAAGTCGCAAAACTTATCGACGTTTCCAcc                               >  1:3311025/1‑36 (MQ=255)
     aaaGCAGAAGTCGCAAAACTTATCGACGTTTACAcc                               >  1:5293043/1‑36 (MQ=255)
         cAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTa                           >  1:4023379/1‑36 (MQ=255)
          acaaGTCGTAAAACTTATCGACTTTTCCACCTGTAt                          <  1:3389289/34‑1 (MQ=21)
          aGAAGTCGCAAAACTTATCGACGTTTCCACCTGTAt                          <  1:1352526/36‑1 (MQ=255)
          aGAAGTCGCAAAACTTATCGACGTTTCCACCTGTAt                          >  1:3301878/1‑36 (MQ=255)
           gAAGTCGCAAAACTTATCGACGTTTCCACCTGTATc                         <  1:3984765/36‑1 (MQ=255)
           gAAGTCGCAAAACTTATCGACGTTTCCACCTGTATc                         >  1:3166514/1‑36 (MQ=255)
           gAAGTCGCAAAAATTATCGACGTTTCCACCTGTATc                         >  1:4787509/1‑36 (MQ=255)
             aGTCGCAAAACTTATCGACGTTTCCACCTGTATCgg                       >  1:3178719/1‑36 (MQ=255)
               tCGCAAAACTTATCGACGTTTCCACCTGTATCGGCt                     <  1:5599756/36‑1 (MQ=255)
                 gCAAAACTTATCGACGTTTCCACCTGTATCGGCTGt                   <  1:5717456/36‑1 (MQ=255)
                 gCAAAACTTATCGACGTTTCCACCGGTATCGGCTgg                   >  1:112182/1‑35 (MQ=255)
                  cAAAACTTATCGACGTTTCCACCTGTATCGGCTGTa                  <  1:1415941/36‑1 (MQ=255)
                   aaaaCTTATCGACGTTTCCACCTGTATCGGCTGTaa                 <  1:3708440/36‑1 (MQ=255)
                   aaaaCTTATCGACGTTTCCACATGTATCGGCTggta                 >  1:3496344/1‑33 (MQ=18)
                    aaaCTTCTCGACGTTTCCACCTGTATCGGCTGTaaa                <  1:2657801/36‑1 (MQ=255)
                       cTTATCGACGTTTCCACCTGTATCGGCTGTAAAGcc             <  1:1316090/36‑1 (MQ=255)
                        ttATCGACGTTTCCACCTGTATCGGCTGTAAAGccc            >  1:1067450/1‑35 (MQ=255)
                        ttATCGACGTTTCCACCTGTATCGGCTGTAAAGCCt            <  1:4229451/36‑1 (MQ=255)
                        ttATCGACGTTTCCACCTGTATCGGCTGTAAAGCCt            <  1:4896877/36‑1 (MQ=255)
                          aTCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGt          <  1:5277361/36‑1 (MQ=38)
                           tCGACGTTTCCACCTGTATCGGCTGTAAAAACTGtt         >  1:5575386/1‑35 (MQ=14)
                           tCGACGTTACCACCTGTATCGGCTGTAAAGTCtttt         >  1:5513135/1‑33 (MQ=14)
                             gACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAg       <  1:1729964/36‑1 (MQ=38)
                             gACGTTTCCACCTGTATCGGCTGGAAAGCCTGttcg       >  1:3477074/1‑33 (MQ=14)
                             gACGTTTACACCTGTATCGGCTGTAAAGACCGTccg       >  1:854135/1‑34 (MQ=14)
                              aCGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAgg      <  1:2152838/36‑1 (MQ=38)
                              aCGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAgg      >  1:199685/1‑36 (MQ=38)
                              aCGTTTCCACCTGTATCGGCTGTAAAACCCTTCAgg      >  1:5610717/1‑36 (MQ=14)
                              aCGTTACCACCTGTATCGGCTGGAAAGCCTTTCCgg      >  1:4435820/1‑36 (MQ=14)
                               cGTTTCAATCTGTGTCGGCTGCAAAGCCTGTCAGGt     <  1:470729/36‑1 (MQ=9)
                                gTTACCACCTGTATCGGCTGTAAAACCTGTTATGTg    >  1:2494692/1‑36 (MQ=14)
                                  ttCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGc  <  1:3576464/36‑1 (MQ=35)
                                   |                                  
ATTACAAAGCAGAAGTCGCAAAACTTATCGACGTTTCCACCTGTATCGGCTGTAAAGCCTGTCAGGTGGC  >  REL606/1524684‑1524753

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG