Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,067 T→C 11.1% V175V (GTT→GTC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,0670TC11.1% 160.0 / 9.9 63V175V (GTT→GTCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (33/23);  new base C (3/4);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

CGCGTGGACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTT  >  minE/831997‑832136
                                                                      |                                                                     
cgcgTGGACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGtt                                                                       <  1:1212092/71‑1 (MQ=255)
      gACGCGGTGGATTTGTTCGTTCATCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTacac                                                                  >  1:651955/1‑70 (MQ=255)
                 tttGTTCGTTCTTCCTGGCATGAGGTGAACTAACTGATCGCCTCATCTAACGtttt                                                                     <  1:5017982/56‑2 (MQ=255)
                 tttGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTa                                                                     <  1:3144920/56‑1 (MQ=255)
                 tttGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTa                                                                     <  1:741651/56‑1 (MQ=255)
                       cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTac                                                                    >  1:2212690/1‑51 (MQ=255)
                       cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTac                                                                    >  1:5235618/1‑51 (MQ=255)
                       cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCaaaaa                                                         >  1:2992601/1‑62 (MQ=255)
                       cgttcTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAAc                                                        >  1:2268023/1‑63 (MQ=255)
                         ttcttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGtt                                                                       <  1:5099408/46‑1 (MQ=255)
                         ttcttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCaaa                                                           <  1:340700/58‑1 (MQ=255)
                           cttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGttt                                                                      >  1:4538960/1‑45 (MQ=255)
                           cttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTa                                                                     >  1:4580473/1‑46 (MQ=255)
                           cttcCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGAc                                            >  1:349094/1‑71 (MQ=255)
                              ccTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGAc                                            <  1:2145428/68‑1 (MQ=255)
                               cTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCg                                          <  1:4029823/69‑1 (MQ=255)
                                  gCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTa                                                                     >  1:1536165/1‑39 (MQ=255)
                                       aggTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGt                                <  1:215467/71‑1 (MQ=255)
                                           gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTAc                                                     >  1:352943/1‑46 (MQ=255)
                                           gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                            >  1:493275/1‑71 (MQ=255)
                                           gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                            >  1:4941474/1‑71 (MQ=255)
                                           gaacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                            >  1:4393584/1‑71 (MQ=255)
                                            aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAATTACGGCCCGGACCGTGTCGCTGGCTtctc                           >  1:4040059/1‑71 (MQ=255)
                                            aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAATTACGGCCCGGACCGTGTCGCTGGCTtctc                           >  1:2628453/1‑71 (MQ=255)
                                            aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtct                            <  1:131151/70‑1 (MQ=255)
                                            aacgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtc                             <  1:5481344/69‑1 (MQ=255)
                                              cgaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGc                         >  1:3485223/1‑71 (MQ=255)
                                               gaacTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGcc                        <  1:4976961/71‑1 (MQ=255)
                                                  cTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAAt                     <  1:4964170/71‑1 (MQ=255)
                                                    gATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGAc                                            >  1:1732775/1‑46 (MQ=255)
                                                     attgccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATccc                  <  1:3152349/68‑1 (MQ=255)
                                                     aTCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCgtgt                                       <  1:4243162/50‑1 (MQ=255)
                                                     aTCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTcc                  >  1:936587/1‑71 (MQ=255)
                                                      tCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCgtg                                        <  1:4168738/48‑1 (MQ=255)
                                                      tCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAAt                     >  1:1866958/1‑67 (MQ=255)
                                                      tCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCg                 >  1:2305049/1‑71 (MQ=255)
                                                       cgccgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTg                                     >  1:3150246/1‑50 (MQ=255)
                                                        gccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGc               >  1:634764/1‑71 (MQ=255)
                                                        gccgccGCTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGc               >  1:4053174/1‑71 (MQ=255)
                                                        gccgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTCGCTGGTTTCTCGCCAATTCCGGc               <  1:2273519/71‑1 (MQ=255)
                                                          cgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa             >  1:123709/1‑71 (MQ=255)
                                                           gcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCgtg                                        <  1:3474836/43‑1 (MQ=255)
                                                            cATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGtt                               >  1:2335341/1‑51 (MQ=255)
                                                            cATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATg           <  1:5015116/71‑1 (MQ=255)
                                                              tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTc                   <  1:4360545/61‑1 (MQ=255)
                                                              tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGc               <  1:993397/65‑1 (MQ=255)
                                                              tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGc               <  1:2610061/65‑1 (MQ=255)
                                                              tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc         >  1:4559810/1‑71 (MQ=255)
                                                              gcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGc               <  1:183876/64‑1 (MQ=255)
                                                              gcTAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTc         <  1:5170591/70‑1 (MQ=255)
                                                                tAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTcc                  >  1:5432170/1‑60 (MQ=255)
                                                                tAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGc               <  1:2670855/63‑1 (MQ=255)
                                                                tAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGtt                               <  1:2404767/47‑1 (MQ=255)
                                                                tAACGTCTGGACCATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTtct                            >  1:5250444/1‑50 (MQ=255)
                                                                 aaCGTTTACACCATCAAAAACTACGGCCCGGACCgtgt                                       >  1:1888533/1‑38 (MQ=255)
                                                                 aaCGTTTACACCATCAAAAACTACGGCCCGGACCgtgt                                       >  1:1158981/1‑38 (MQ=255)
                                                                 aaCGTTTACACCATCAAAAACTACGGCCCGGACCgtg                                        >  1:3044162/1‑37 (MQ=255)
                                                                 aaCGTTTACACCATCAAAAACTACGGCCCGGACCgtg                                        >  1:1733049/1‑37 (MQ=255)
                                                                 aaCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCg        >  1:517407/1‑69 (MQ=255)
                                                                 aaCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGAt      >  1:317158/1‑71 (MQ=255)
                                                                  aCGTTTACACCATCAAAAACTACGGCCCGGACCgtgt                                       <  1:2798108/37‑1 (MQ=255)
                                                                  aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCg                 <  1:3344827/59‑1 (MQ=255)
                                                                  aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa             >  1:2961017/1‑63 (MQ=255)
                                                                  aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc         <  1:2022767/67‑1 (MQ=255)
                                                                  aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc         <  1:3704031/67‑1 (MQ=255)
                                                                   cGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTtctc                           >  1:2013047/1‑48 (MQ=255)
                                                                    gTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGt   >  1:1435766/1‑71 (MQ=255)
                                                                     tttACACCATCAAAAACTACGGCCCGGACCGTGTTGc                                    <  1:1025657/37‑1 (MQ=255)
                                                                     tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCt                                   <  1:5180480/38‑1 (MQ=255)
                                                                     tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCt                                   <  1:3931721/38‑1 (MQ=255)
                                                                     tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCt                                   <  1:3938579/38‑1 (MQ=255)
                                                                      ttACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt  >  1:3071241/1‑70 (MQ=255)
                                                                      |                                                                     
CGCGTGGACGCGGTGGATTTGTTCGTTCTTCCTGGCAGGAGGTGAACGAACTGATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTT  >  minE/831997‑832136

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: