Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,442 G→A 10.4% 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,4420GA10.4% 235.9 / 13.3 86Q300Q (CAG→CAAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base A (4/5);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCC  >  minE/832377‑832510
                                                                 |                                                                    
gTTGCCGTCACACCAGACTACGCTGAAAACGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAc                                                                 >  1:4448240/1‑71 (MQ=255)
 ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                                     <  1:669984/66‑1 (MQ=255)
 ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                                     <  1:3690683/66‑1 (MQ=255)
 ttGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                                                <  1:1910169/71‑1 (MQ=255)
    ccGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGAt                                                             <  1:3196267/71‑1 (MQ=255)
       tCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGa                                                              <  1:3866823/67‑1 (MQ=255)
        cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGc                                                                   >  1:2585835/1‑61 (MQ=255)
        cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                                         >  1:2149322/1‑71 (MQ=255)
        cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                                         >  1:2155976/1‑71 (MQ=255)
        cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCgg                                                         >  1:4275793/1‑71 (MQ=255)
         acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAgg                                                                     >  1:592428/1‑58 (MQ=255)
         acacCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGc                                                        <  1:5383525/71‑1 (MQ=255)
           acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCaa                                                      <  1:3964059/71‑1 (MQ=255)
             cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                    >  1:4531944/1‑71 (MQ=255)
             cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                    >  1:5291713/1‑71 (MQ=255)
              aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                   >  1:1044184/1‑71 (MQ=255)
              aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                   <  1:3169226/71‑1 (MQ=255)
              aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                   <  1:4082993/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:2056672/56‑1 (MQ=255)
                aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                                                <  1:2641041/56‑1 (MQ=255)
                aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                 >  1:3504627/1‑71 (MQ=255)
                aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                 >  1:3295111/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)
                  tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                               >  1:1171134/1‑71 (MQ=255)
                    cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                             >  1:4875849/1‑71 (MQ=255)
                    cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                             <  1:2331287/71‑1 (MQ=255)
                     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)
                       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:1202196/47‑1 (MQ=255)
                         aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCAcc                                                                <  1:4716820/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:4135535/1‑71 (MQ=255)
                         aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg                                        <  1:4550437/71‑1 (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)
                               ccAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCat                                  <  1:4956967/71‑2 (MQ=255)
                                cAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCatg                                 >  1:972398/1‑69 (MQ=255)
                                 aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt                                <  1:3292348/71‑1 (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:1951876/1‑49 (MQ=255)
                                         cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                              >  1:1306486/1‑49 (MQ=255)
                                         cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                              >  1:4667873/1‑49 (MQ=255)
                                         cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                              >  1:711005/1‑49 (MQ=255)
                                         cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                            <  1:950805/51‑1 (MQ=255)
                                         cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                            <  1:3870600/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)
                                                tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                     <  1:4393558/67‑1 (MQ=255)
                                                tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                    <  1:3415584/68‑1 (MQ=255)
                                                 ggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                     <  1:5047364/66‑1 (MQ=255)
                                                  gctggcACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATg                                        >  1:3183114/1‑46 (MQ=255)
                                                       cACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                             <  1:5106141/36‑1 (MQ=255)
                                                        aCCGAAACAGGGCTCCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                    <  1:545220/60‑1 (MQ=255)
                                                         ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                <  1:3193421/47‑1 (MQ=255)
                                                         ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                <  1:689026/47‑1 (MQ=255)
                                                         ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                <  1:2098150/47‑1 (MQ=255)
                                                         ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc          <  1:2653088/69‑1 (MQ=255)
                                                          cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGa        >  1:991128/1‑70 (MQ=255)
                                                          cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc       <  1:2245183/71‑1 (MQ=255)
                                                           gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtg                            >  1:1505257/1‑49 (MQ=255)
                                                           gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc       >  1:2577609/1‑70 (MQ=255)
                                                             aaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa     >  1:5637108/1‑70 (MQ=255)
                                                              aCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                    <  1:4273547/54‑1 (MQ=255)
                                                               cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc       <  1:2510651/66‑1 (MQ=255)
                                                               cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc  <  1:3708554/71‑1 (MQ=255)
                                                                aGGGCACCGATGCGTCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc       <  1:3677352/65‑1 (MQ=255)
                                                                 |                                                                    
GTTGCCGTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCC  >  minE/832377‑832510

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: