Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,230 C→A 19.0% Q230K (CAA→AAA) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,2300CA19.0% 175.8 / 20.5 68Q230K (CAA→AAA) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (20/35);  new base A (6/7);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTG  >  minE/832161‑832300
                                                                     |                                                                      
ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                                        >  1:1034909/1‑70 (MQ=255)
ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                                        >  1:2059903/1‑70 (MQ=255)
ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                                                        >  1:1599583/1‑70 (MQ=255)
   tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                                      <  1:2128473/69‑1 (MQ=255)
   tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                                                    >  1:2177040/1‑71 (MQ=255)
   tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                                                    >  1:3680238/1‑71 (MQ=255)
      cTGATTGGCTGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTg                                                                 <  1:3590913/71‑1 (MQ=255)
      cTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTg                                                                 <  1:1659438/71‑1 (MQ=255)
        gATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                                                               >  1:1951126/1‑71 (MQ=255)
         aTTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTgggg                                                              >  1:3435915/1‑71 (MQ=255)
          ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                                                             >  1:2810613/1‑71 (MQ=255)
          ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                                                             <  1:3029009/71‑1 (MQ=255)
           ttgCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCg                                                            <  1:3865011/69‑1 (MQ=255)
             gCGGTACTTGCTTAAGCTTCTACGACTGGTACAGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAg                                                          <  1:2912485/71‑1 (MQ=255)
             gCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAg                                                          >  1:1652358/1‑71 (MQ=255)
               ggTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCa                                                                       >  1:2643338/1‑56 (MQ=255)
                gTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCa                                                        >  1:3821945/1‑70 (MQ=255)
                 tACTTGCTTAAGCTTCTACGACTCGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTggg                                                               <  1:869607/62‑1 (MQ=255)
                   cTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACt                                                    <  1:3917703/71‑1 (MQ=255)
                    ttGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaaa                                                                     >  1:2709945/1‑53 (MQ=255)
                      gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTg                                                   <  1:700467/69‑1 (MQ=255)
                       cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTg                                                   <  1:1064268/68‑1 (MQ=255)
                       cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                                                <  1:1482933/71‑1 (MQ=255)
                         tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa                                              >  1:1842024/1‑71 (MQ=255)
                             cTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCTCAAACCTGGGGCGAGCAAACTGACGTACCgg                                          <  1:375734/71‑1 (MQ=255)
                             cTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTg                                                                 <  1:1304482/48‑1 (MQ=255)
                               tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGc                                                             <  1:1338106/50‑1 (MQ=255)
                               tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCa                                                        <  1:3075661/55‑1 (MQ=255)
                               tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGa                                         >  1:2231834/1‑70 (MQ=255)
                               tCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGa                                         <  1:3881838/70‑1 (MQ=255)
                               tCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGaa                                        >  1:902386/1‑71 (MQ=255)
                                cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                                      <  1:3819244/40‑1 (MQ=255)
                                cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                                                      <  1:2584082/40‑1 (MQ=255)
                                cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                         <  1:877986/69‑1 (MQ=255)
                                cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                       >  1:1349517/1‑71 (MQ=255)
                                 tACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGa                                                           >  1:2636449/1‑50 (MQ=255)
                                 tACGACTGGTATTGCGACCTGCCGCCCGCATCGCCGATGACCTGGGGCg                                                            <  1:2179961/49‑1 (MQ=255)
                                  aCGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                       >  1:3906156/1‑69 (MQ=255)
                                   cGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATctg                                    >  1:2893140/1‑71 (MQ=255)
                                   cGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATct                                     <  1:3163374/70‑1 (MQ=255)
                                   cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                       <  1:3273029/68‑1 (MQ=255)
                                      cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCg                                                            <  1:2402401/44‑1 (MQ=255)
                                      cTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctgctg                                 <  1:2482226/71‑1 (MQ=255)
                                           attGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATctgc                                   <  1:3019215/62‑1 (MQ=255)
                                           attGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTgg                            >  1:1164402/3‑71 (MQ=255)
                                             tGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                         <  1:3892563/56‑1 (MQ=255)
                                              gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                         <  1:3828782/55‑1 (MQ=255)
                                              gCGACTTGCCTCATGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTa                          >  1:1000755/1‑70 (MQ=255)
                                               cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                         <  1:1217644/70‑1 (MQ=255)
                                               cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                         <  1:3645032/70‑1 (MQ=255)
                                               cGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTat                         <  1:3768923/70‑2 (MQ=255)
                                                   ttGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc                    >  1:2335376/1‑71 (MQ=255)
                                                     gCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTgg                            <  1:2292016/61‑1 (MQ=255)
                                                     gCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTgg                            <  1:313955/61‑1 (MQ=255)
                                                      cctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa                  >  1:1639362/1‑70 (MQ=255)
                                                        tcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa                  <  1:2042996/68‑1 (MQ=255)
                                                        tcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa                  <  1:2338139/68‑1 (MQ=255)
                                                          ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAc             <  1:2468480/71‑1 (MQ=255)
                                                            gCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat           >  1:2493073/1‑71 (MQ=255)
                                                             cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat           >  1:2353916/1‑70 (MQ=255)
                                                             cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAca            <  1:3653117/69‑1 (MQ=255)
                                                              gTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatca         <  1:784003/71‑1 (MQ=255)
                                                              gTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatc          >  1:3764119/1‑70 (MQ=255)
                                                               tctcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTgg                            <  1:971659/51‑1 (MQ=255)
                                                                  ccGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCt               >  1:2402554/1‑61 (MQ=255)
                                                                  ccGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACatcatc       <  1:2893811/69‑1 (MQ=255)
                                                                  ccGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc     >  1:3814001/1‑71 (MQ=255)
                                                                   cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctgctg                                 <  1:678182/42‑1 (MQ=255)
                                                                    gCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc     >  1:907856/1‑69 (MQ=255)
                                                                    gCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCt   <  1:700540/71‑1 (MQ=255)
                                                                    gCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCt   <  1:844787/71‑1 (MQ=255)
                                                                     cAAACCTGGGGCGAGCAAACTGACGTACCGTAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg  <  1:1783589/71‑1 (MQ=255)
                                                                     cAAACCTGGGGCGAGCAAACTGACGTACCGTAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg  <  1:2627818/71‑1 (MQ=255)
                                                                     cAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg  >  1:457520/1‑71 (MQ=255)
                                                                     |                                                                      
CTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTG  >  minE/832161‑832300

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: