Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,423 T→A 12.2% L294Q (CTG→CAG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,4230TA12.2% 241.1 / 22.9 90L294Q (CTG→CAG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base A (5/6);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGT  >  minE/832354‑832490
                                                                     |                                                                   
tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCt                                                                     <  1:219463/70‑1 (MQ=255)
  cGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCtgt                                                                   <  1:1111121/70‑1 (MQ=255)
   gTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtgg                                                                 <  1:2952137/71‑1 (MQ=255)
   gTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtgg                                                                 <  1:3665311/71‑1 (MQ=255)
    ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc                                                                <  1:1655177/71‑1 (MQ=255)
    ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc                                                                <  1:3627765/71‑1 (MQ=255)
    ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc                                                                <  1:2074091/71‑1 (MQ=255)
    ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc                                                                <  1:884587/71‑1 (MQ=255)
      aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctg                                                              >  1:2784215/1‑71 (MQ=255)
      aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctg                                                              >  1:278699/1‑71 (MQ=255)
      aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctg                                                              >  1:456901/1‑71 (MQ=255)
         aaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctg                                                              <  1:3156765/68‑1 (MQ=255)
           gggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                                                         <  1:125354/71‑1 (MQ=255)
           gggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                                                         <  1:1716563/71‑1 (MQ=255)
             gACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCg                                                        >  1:106184/1‑70 (MQ=255)
              aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCtgtg                                                                  <  1:2162351/59‑1 (MQ=255)
               ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                                                     <  1:829271/71‑1 (MQ=255)
               ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                                                     <  1:2526951/71‑1 (MQ=255)
                cAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                                             <  1:3288342/62‑1 (MQ=255)
                 aaaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctggc                                                            >  1:934550/1‑62 (MQ=255)
                 aaaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                                                   <  1:1192691/71‑1 (MQ=255)
                  aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCa                                                           >  1:2768099/1‑62 (MQ=255)
                   aaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                 >  1:1602109/1‑71 (MQ=255)
                   aaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                 >  1:2216480/1‑71 (MQ=255)
                    aCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCtgt                                                                   >  1:3613883/1‑52 (MQ=255)
                    aCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCtgt                                                                   >  1:930493/1‑52 (MQ=255)
                      tGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCa                                              >  1:2970708/1‑71 (MQ=255)
                        ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                            >  1:3620646/1‑71 (MQ=255)
                          gCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAg                                                  >  1:2231822/1‑63 (MQ=255)
                            cGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGAt                                         <  1:350066/70‑1 (MQ=255)
                            cGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATg                                        <  1:708412/71‑1 (MQ=255)
                             gTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGc                                       <  1:1945686/71‑1 (MQ=255)
                               cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                                       <  1:2179610/53‑1 (MQ=255)
                                acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGc                                       >  1:1565246/1‑68 (MQ=255)
                                acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                     >  1:1494829/1‑70 (MQ=255)
                                   ccAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAc                                                          >  1:3241691/1‑46 (MQ=255)
                                    cAGACTACGCTGCAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAGACAGGGCACCGATGCGGCAATg                                >  1:1818749/1‑71 (MQ=255)
                                    cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                <  1:2466817/71‑1 (MQ=255)
                                    cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                >  1:384313/1‑71 (MQ=255)
                                    cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                >  1:1469432/1‑71 (MQ=255)
                                     aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               <  1:2809393/71‑1 (MQ=255)
                                      gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               >  1:3039408/1‑70 (MQ=255)
                                       actactCTGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCa                                              <  1:3502691/54‑1 (MQ=19)
                                       actactCTGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCa                                              <  1:816732/54‑1 (MQ=19)
                                         tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                              >  1:2164379/1‑68 (MQ=255)
                                         tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                              >  1:3805538/1‑68 (MQ=255)
                                         tACGCTGAAATCGCCAAACTGTGAGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                              >  1:1175749/1‑68 (MQ=255)
                                           cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                           <  1:2326459/69‑1 (MQ=255)
                                            gCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                        <  1:1074059/71‑1 (MQ=255)
                                              tGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCa                                              >  1:3584605/1‑47 (MQ=19)
                                              tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGc                                               >  1:3692799/1‑46 (MQ=255)
                                              tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCa                                              >  1:1536119/1‑47 (MQ=255)
                                                aaaTCGCCTAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                     <  1:1715310/70‑1 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc                                             >  1:1329575/1‑46 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc                                             >  1:2850585/1‑46 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                     <  1:3213067/70‑1 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                     <  1:2973687/70‑1 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                     <  1:966293/70‑1 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                     <  1:1942153/70‑1 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                     <  1:2522704/70‑1 (MQ=255)
                                                 aaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATgg                   >  1:1098882/1‑71 (MQ=255)
                                                  aTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggg                  >  1:2943493/1‑71 (MQ=255)
                                                   tCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGc                 <  1:3717259/71‑1 (MQ=255)
                                                    cGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                         >  1:1321583/1‑62 (MQ=255)
                                                     gCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGcc                <  1:1701434/70‑1 (MQ=255)
                                                     gCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCa               >  1:3053525/1‑71 (MQ=255)
                                                      ccAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCat              <  1:418280/71‑2 (MQ=255)
                                                       cAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCActg                                           <  1:1583995/41‑3 (MQ=19)
                                                       cAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCatg             >  1:2124485/1‑69 (MQ=255)
                                                        aaaTTGTGCGACCAGTGGCTGGCACCGAAACAagg                                                <  1:2224396/35‑1 (MQ=16)
                                                        aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCa                                              >  1:1999934/1‑37 (MQ=255)
                                                        aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt            <  1:3310625/71‑1 (MQ=255)
                                                        aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt            <  1:1055514/71‑1 (MQ=255)
                                                         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                     >  1:3767962/1‑45 (MQ=255)
                                                         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCAcct            >  1:1642185/1‑68 (MQ=255)
                                                         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACg             >  1:1401032/1‑69 (MQ=255)
                                                         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt            >  1:1703189/1‑70 (MQ=255)
                                                         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt            >  1:3090243/1‑70 (MQ=255)
                                                         aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt            >  1:1569310/1‑70 (MQ=255)
                                                            tgtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               <  1:2401328/48‑1 (MQ=255)
                                                            tgtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               <  1:2405978/48‑1 (MQ=255)
                                                            tgtgCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGTGAtt        >  1:2941500/1‑70 (MQ=255)
                                                             gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgc       >  1:1987441/1‑71 (MQ=255)
                                                             gtgCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGc                 >  1:2269038/1‑61 (MQ=255)
                                                              tgCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGTGAttt       >  1:596674/1‑68 (MQ=255)
                                                                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc    <  1:3913534/71‑1 (MQ=255)
                                                                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc    <  1:2682659/71‑1 (MQ=255)
                                                                 gATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               <  1:1933014/43‑1 (MQ=255)
                                                                  aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt  >  1:3609478/1‑71 (MQ=255)
                                                                  aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt  >  1:3799497/1‑71 (MQ=255)
                                                                   tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                           <  1:1039792/45‑1 (MQ=255)
                                                                   tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                           <  1:3653872/45‑1 (MQ=255)
                                                                   tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                           <  1:3353623/45‑1 (MQ=255)
                                                                   tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCg   >  1:2492217/1‑69 (MQ=255)
                                                                   tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt  >  1:1233680/1‑70 (MQ=255)
                                                                     tgtgGCTGGCACCGAAACATGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAt         >  1:2915741/1‑61 (MQ=255)
                                                                     |                                                                   
TGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGT  >  minE/832354‑832490

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: