Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,413 C→T 11.8% L291L (CTG→TTG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,4130CT11.8% 296.7 / 22.0 111L291L (CTG→TTG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (43/54);  new base T (8/5);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 118 total aligned reads displayed.

CTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGT  >  minE/832343‑832480
                                                                      |                                                                   
cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                                                     >  1:270543/1‑71 (MQ=255)
cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                                                     >  1:3395995/1‑71 (MQ=255)
cttTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                                                     >  1:1454583/1‑71 (MQ=255)
   tACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtg                                                                   <  1:2003354/70‑1 (MQ=255)
    aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgtg                                                                 >  1:5422094/1‑71 (MQ=255)
    aCTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                                                  >  1:88103/1‑70 (MQ=255)
     cTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGc                                                                <  1:5286533/71‑1 (MQ=255)
       gAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGa                                                              >  1:5670754/1‑71 (MQ=255)
        aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg                                                               <  1:1851614/69‑1 (MQ=255)
        aaGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCg                                                               <  1:1822553/69‑1 (MQ=255)
         aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc                                                            >  1:1294859/1‑71 (MQ=255)
         aGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATc                                                            <  1:1794494/71‑1 (MQ=255)
           tGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGa                                                              <  1:517313/67‑1 (MQ=255)
               ttACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACtgt                                                                  <  1:617241/59‑1 (MQ=255)
                 aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggct                                                     <  1:2914552/70‑1 (MQ=255)
                   aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                                   >  1:4894099/1‑70 (MQ=255)
                   aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                                   >  1:2249693/1‑70 (MQ=255)
                    aaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                                   <  1:4086787/69‑1 (MQ=255)
                     aGGGACCGAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAc                                                <  1:572442/71‑1 (MQ=255)
                       ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                                               >  1:1196310/1‑70 (MQ=255)
                       ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                                               >  1:4125620/1‑70 (MQ=255)
                        gACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                             <  1:403498/71‑1 (MQ=255)
                        gACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                             <  1:2448142/71‑1 (MQ=255)
                         aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaa                                            <  1:4646452/71‑1 (MQ=255)
                         aCCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaa                                            <  1:2999605/71‑1 (MQ=255)
                          ccAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                                           <  1:4162029/71‑1 (MQ=255)
                            aaaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                                           >  1:1259197/1‑69 (MQ=255)
                             aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACa                                         >  1:5118343/1‑70 (MQ=255)
                              aaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAAc                                                                     >  1:2172988/1‑41 (MQ=255)
                                  gTTGCCGTCACACCAGACTACGCTGAAAACGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc                                   >  1:4448240/1‑71 (MQ=255)
                                   ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                       <  1:3690683/66‑1 (MQ=255)
                                   ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                  <  1:1910169/71‑1 (MQ=255)
                                      ccGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGAt                               <  1:3196267/71‑1 (MQ=255)
                                         tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                                         <  1:3756818/58‑1 (MQ=255)
                                         tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                                         <  1:894360/58‑1 (MQ=255)
                                         tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGa                                <  1:3866823/67‑1 (MQ=255)
                                          cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                                          >  1:2355399/1‑40 (MQ=255)
                                          cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCtg                                                          >  1:2958692/1‑40 (MQ=255)
                                          cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGc                                     >  1:2585835/1‑61 (MQ=255)
                                          cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                           >  1:4275793/1‑71 (MQ=255)
                                          cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                           >  1:2149322/1‑71 (MQ=255)
                                          cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                           >  1:2155976/1‑71 (MQ=255)
                                           acacCAGACTACGCTGAAATCGCCAAACTGTGCGAt                                                             >  1:4251374/1‑36 (MQ=255)
                                           acacCAGACTACGCTGAAATCGCCAAACTGTGCGAt                                                             >  1:3032130/1‑36 (MQ=255)
                                           acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                                           <  1:2837562/54‑1 (MQ=255)
                                           acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGaaa                                           <  1:2830678/54‑1 (MQ=255)
                                           acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                       >  1:592428/1‑58 (MQ=255)
                                           acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGc                          <  1:5383525/71‑1 (MQ=255)
                                             acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACa                                         <  1:2392515/54‑1 (MQ=255)
                                             acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCaa                        <  1:3964059/71‑1 (MQ=255)
                                              ccAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                                   <  1:4590889/43‑1 (MQ=255)
                                              ccAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAAc                                          <  1:5030950/52‑1 (MQ=255)
                                               cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc                                                      >  1:169297/1‑39 (MQ=255)
                                               cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggc                                                      >  1:5126283/1‑39 (MQ=255)
                                               cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGTACCGAAAc                                          >  1:532683/1‑51 (MQ=255)
                                               cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAAc                                          >  1:3821236/1‑51 (MQ=255)
                                               cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAAc                                          >  1:4357905/1‑51 (MQ=255)
                                               cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                      >  1:4531944/1‑71 (MQ=255)
                                                aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                     >  1:1044184/1‑71 (MQ=255)
                                                 gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                    <  1:2195511/71‑1 (MQ=255)
                                                  aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                  <  1:4579405/56‑1 (MQ=255)
                                                  aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                  <  1:2506091/56‑1 (MQ=255)
                                                  aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                  <  1:2641041/56‑1 (MQ=255)
                                                  aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                   >  1:3295111/1‑71 (MQ=255)
                                                  aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                   >  1:4369683/1‑71 (MQ=255)
                                                  aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                   >  1:2937726/1‑71 (MQ=255)
                                                  aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                    >  1:5640309/1‑70 (MQ=255)
                                                    tACTCTGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCAct                                  <  1:2780895/54‑2 (MQ=19)
                                                    tACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                                   <  1:2699977/37‑1 (MQ=255)
                                                    tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                             <  1:712533/43‑1 (MQ=255)
                                                    tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                             <  1:667677/43‑1 (MQ=255)
                                                    tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                 >  1:1171134/1‑71 (MQ=255)
                                                      cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg               >  1:4875849/1‑71 (MQ=255)
                                                       tcTGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAAc                                          <  1:1324339/42‑1 (MQ=16)
                                                       gCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc              >  1:4841478/1‑71 (MQ=255)
                                                        cTGAATTCGCCAAAGTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCg             <  1:2510370/71‑1 (MQ=255)
                                                        cTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGc                          <  1:771621/58‑1 (MQ=255)
                                                        cTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGc                          <  1:3317750/58‑1 (MQ=255)
                                                         tGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAAc                                          >  1:1197339/1‑41 (MQ=16)
                                                         tGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAAc                                          >  1:5243065/1‑41 (MQ=16)
                                                         tGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAAc                                          >  1:5341849/1‑41 (MQ=16)
                                                         tGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCAct                                  >  1:2389869/1‑48 (MQ=19)
                                                         tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                     <  1:2759535/62‑1 (MQ=255)
                                                         tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                 <  1:471326/66‑1 (MQ=255)
                                                           aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                  <  1:4716820/47‑1 (MQ=255)
                                                           aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                  <  1:1202196/47‑1 (MQ=255)
                                                           aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa            <  1:30606/69‑1 (MQ=255)
                                                           aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt           <  1:2656352/70‑1 (MQ=255)
                                                           aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg          <  1:4550437/71‑1 (MQ=255)
                                                           aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg          >  1:4135535/1‑71 (MQ=255)
                                                            aaTCACCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCg                                 <  1:216450/47‑1 (MQ=255)
                                                             aTCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggt        >  1:1041489/1‑70 (MQ=255)
                                                             aTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggg        >  1:1459041/1‑71 (MQ=255)
                                                                gCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCa     >  1:5162204/1‑71 (MQ=255)
                                                                 ccAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCat    <  1:4956967/71‑2 (MQ=255)
                                                                   aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt  <  1:3292348/71‑1 (MQ=255)
                                                                   aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt  <  1:4691403/71‑1 (MQ=255)
                                                                   aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt  <  1:647448/71‑1 (MQ=255)
                                                                      cTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                     >  1:726344/1‑49 (MQ=255)
                                                                      cTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                     >  1:234138/1‑49 (MQ=255)
                                                                      |                                                                   
CTTTACTGAAGTGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGT  >  minE/832343‑832480

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