Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,423 T→A 13.1% 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
*Exported832,4230TA13.1% 285.5 / 28.0 106L294Q (CTG→CAG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (43/49);  new base A (8/6);  total (51/55)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGT  >  Exported/832354‑832490
                                                                     |                                                                   
      aCAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggct                                                               <  1:2914552/70‑1 (MQ=255)
        aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                                             >  1:2249693/1‑70 (MQ=255)
        aaaGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctgg                                                             >  1:4894099/1‑70 (MQ=255)
          aGGGACCGAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAc                                                          <  1:572442/71‑1 (MQ=255)
            ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                                                         >  1:4125620/1‑70 (MQ=255)
            ggACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCAcc                                                         >  1:1196310/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)
                  aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGtggctggc                                                            <  1:2482576/61‑1 (MQ=255)
                  aaaCTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACa                                                   >  1:5118343/1‑70 (MQ=255)
                       gTTGCCGTCACACCAGACTACGCTGAAAACGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc                                             >  1:4448240/1‑71 (MQ=255)
                        ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                 <  1:3690683/66‑1 (MQ=255)
                        ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                 <  1:669984/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:2958692/1‑40 (MQ=255)
                               cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGc                                               >  1:2585835/1‑61 (MQ=255)
                               cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                     >  1:2155976/1‑71 (MQ=255)
                               cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                     >  1:4275793/1‑71 (MQ=255)
                               cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                     >  1:2149322/1‑71 (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: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)
                                    cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                >  1:5291713/1‑71 (MQ=255)
                                     aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               <  1:4082993/71‑1 (MQ=255)
                                     aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               >  1:1044184/1‑71 (MQ=255)
                                     aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               <  1:3169226/71‑1 (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)
                                       aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                            <  1:2056672/56‑1 (MQ=255)
                                       aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                             >  1:3295111/1‑71 (MQ=255)
                                       aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                             >  1:3504627/1‑71 (MQ=255)
                                       aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                             >  1:2937726/1‑71 (MQ=255)
                                       aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                             >  1:4369683/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:667677/43‑1 (MQ=255)
                                         tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGa                                                       <  1:712533/43‑1 (MQ=255)
                                         tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                           >  1:1171134/1‑71 (MQ=255)
                                           cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                         <  1:2331287/71‑1 (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:5243065/1‑41 (MQ=16)
                                              tGAAGTGGCCAAATTGTGCGACCAGTGGCTGGCACCGAAAc                                                    >  1:1197339/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)
                                              tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                      >  1:4293305/1‑71 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                            <  1:4716820/47‑1 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                            <  1:1202196/47‑1 (MQ=255)
                                                aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                             >  1:4746651/1‑62 (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)
                                                 aaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCg                                           <  1:1224752/47‑1 (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)
                                                        aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt            <  1:4691403/71‑1 (MQ=255)
                                                        aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt            <  1:647448/71‑1 (MQ=255)
                                                          attgtgCGACCAGTGGCTGGCACCGAAACAAGGCACTGAt                                         <  1:2578175/38‑1 (MQ=16)
                                                           cTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               >  1:234138/1‑49 (MQ=255)
                                                           cTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                               >  1:726344/1‑49 (MQ=255)
                                                                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                          >  1:1306486/1‑49 (MQ=255)
                                                                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                          >  1:1951876/1‑49 (MQ=255)
                                                                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                          >  1:4667873/1‑49 (MQ=255)
                                                                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                          >  1:711005/1‑49 (MQ=255)
                                                                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                        <  1:3870600/51‑1 (MQ=255)
                                                                cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                        <  1:950805/51‑1 (MQ=255)
                                                                   tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt  >  1:1595232/1‑70 (MQ=255)
                                                                    cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                              >  1:2328187/1‑41 (MQ=255)
                                                                    cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                              >  1:4393556/1‑41 (MQ=255)
                                                                    cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                      >  1:5161605/1‑49 (MQ=255)
                                                                     |                                                                   
TGCGTTACAAAGGGACCAAAACTGTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGT  >  Exported/832354‑832490

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: