Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,442 G→A 12.3% Q300Q (CAG→CAA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,4420GA12.3% 220.4 / 17.1 82Q300Q (CAG→CAAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (39/32);  new base A (5/5);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

AACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAA  >  minE/832373‑832512
                                                                     |                                                                      
aaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                                       >  1:1602109/1‑71 (MQ=255)
aaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                                       >  1:2216480/1‑71 (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)
             acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGc                                                             >  1:1565246/1‑68 (MQ=255)
             acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                                           >  1:1494829/1‑70 (MQ=255)
                 cAGACTACGCTGCAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAGACAGGGCACCGATGCGGCAATg                                                      >  1:1818749/1‑71 (MQ=255)
                 cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                      >  1:1469432/1‑71 (MQ=255)
                 cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                      <  1:2466817/71‑1 (MQ=255)
                 cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                      >  1:384313/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:2850585/1‑46 (MQ=255)
                             aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc                                                                   >  1:1329575/1‑46 (MQ=255)
                             aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                           <  1:2973687/70‑1 (MQ=255)
                             aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                           <  1:2522704/70‑1 (MQ=255)
                             aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                           <  1:3213067/70‑1 (MQ=255)
                             aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                           <  1:1942153/70‑1 (MQ=255)
                             aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                           <  1:966293/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:1055514/71‑1 (MQ=255)
                                     aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                  <  1:3310625/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:1569310/1‑70 (MQ=255)
                                      aaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                  >  1:3090243/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:3799497/1‑71 (MQ=255)
                                               aTCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt                        >  1:3609478/1‑71 (MQ=255)
                                                tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                 <  1:3353623/45‑1 (MQ=255)
                                                tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                 <  1:1039792/45‑1 (MQ=255)
                                                tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                 <  1:3653872/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)
                                                    tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCg                         <  1:2416428/65‑1 (MQ=255)
                                                    tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                      <  1:986865/68‑1 (MQ=255)
                                                    tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGACACGTAATGCTGCGTGa                      <  1:2115255/68‑1 (MQ=255)
                                                    tggctgacACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAAtt                   >  1:802/1‑71 (MQ=255)
                                                     ggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTc                  >  1:2049401/1‑71 (MQ=255)
                                                        tggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCAc               >  1:2446409/1‑71 (MQ=255)
                                                          gcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCt             <  1:2313349/71‑1 (MQ=255)
                                                          gcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCt             >  1:765298/1‑71 (MQ=255)
                                                             ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                            <  1:2477854/37‑1 (MQ=255)
                                                               gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa        <  1:2194359/71‑1 (MQ=255)
                                                               gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACa        <  1:843796/71‑1 (MQ=255)
                                                                aaaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa       >  1:1622538/1‑71 (MQ=255)
                                                                aaaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa       >  1:1302991/1‑71 (MQ=255)
                                                                 aaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCt             >  1:2583947/1‑64 (MQ=255)
                                                                   cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc    <  1:3215379/71‑1 (MQ=255)
                                                                   cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAcc     >  1:2465597/1‑70 (MQ=255)
                                                                     gggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCaa  <  1:918481/71‑1 (MQ=255)
                                                                     |                                                                      
AACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAA  >  minE/832373‑832512

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: