Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,205 C→T 17.8% Y221Y (TAC→TAT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,2050CT17.8% 207.4 / 33.9 90Y221Y (TAC→TATnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base T (9/7);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTG  >  minE/832137‑832269
                                                                    |                                                                
tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                                                                 <  1:1480589/70‑1 (MQ=255)
tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg                                                                >  1:2595552/1‑71 (MQ=255)
    aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGc                                                               <  1:1296926/68‑1 (MQ=255)
    aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGa                                                             <  1:2233866/70‑1 (MQ=255)
    aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGa                                                             <  1:1987235/70‑1 (MQ=255)
     cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc                                                            <  1:3524280/70‑1 (MQ=255)
     cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt                                                           <  1:2150949/71‑1 (MQ=255)
     cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt                                                           <  1:1830429/71‑1 (MQ=255)
      gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                          >  1:917353/1‑71 (MQ=255)
      gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                          >  1:1258687/1‑71 (MQ=255)
        aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGc                                                        >  1:396847/1‑71 (MQ=255)
          cGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc                                                       >  1:246326/1‑70 (MQ=255)
           gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                                     >  1:2472942/1‑71 (MQ=255)
           gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                                     <  1:3345869/71‑1 (MQ=255)
           gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                                     <  1:2526776/71‑1 (MQ=255)
           gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcccc                                                     <  1:1947401/71‑3 (MQ=255)
            ggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACCTGcctcc                                                    >  1:2531308/1‑71 (MQ=255)
             gTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGTTTCTACGACTGGTACTGCGACTTGcctcct                                                   <  1:2330743/71‑1 (MQ=255)
              tGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc                                                            <  1:3529190/61‑1 (MQ=255)
              tGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt                                                           <  1:3922766/62‑1 (MQ=255)
                cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                          >  1:1317014/1‑61 (MQ=255)
                cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg                                                >  1:1408826/1‑71 (MQ=255)
                cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg                                                >  1:2483430/1‑71 (MQ=255)
                cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg                                                <  1:2010019/71‑1 (MQ=255)
                cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg                                                >  1:2843832/1‑71 (MQ=255)
                  cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                              >  1:1970190/1‑71 (MQ=255)
                  cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                              <  1:1525087/71‑1 (MQ=255)
                   gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                             >  1:3162483/1‑71 (MQ=255)
                   gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                             >  1:3778550/1‑71 (MQ=255)
                   gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                             >  1:510064/1‑71 (MQ=255)
                   gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                             >  1:1003353/1‑71 (MQ=255)
                   gCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                              >  1:2269785/1‑70 (MQ=255)
                    cTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                            <  1:1105305/71‑1 (MQ=255)
                      aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGccgcccgc                                                 >  1:3501843/1‑64 (MQ=255)
                      aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGcc                                           >  1:3292544/1‑70 (MQ=255)
                      aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGcc                                           <  1:2645165/70‑1 (MQ=255)
                       tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGcc                                           <  1:3875853/69‑1 (MQ=255)
                        ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                         >  1:1599583/1‑70 (MQ=255)
                        ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                         >  1:1034909/1‑70 (MQ=255)
                        ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                         >  1:2059903/1‑70 (MQ=255)
                          gtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCgat                                        >  1:214126/2‑67 (MQ=255)
                           tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                       <  1:2128473/69‑1 (MQ=255)
                           tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                     >  1:3680238/1‑71 (MQ=255)
                           tcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAAc                                     >  1:2177040/1‑71 (MQ=255)
                             gCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                            <  1:3811769/62‑1 (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)
                                       ggTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                                     <  1:3037837/43‑1 (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)
                                             tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                             >  1:247265/1‑45 (MQ=255)
                                              gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTg                    <  1:700467/69‑1 (MQ=255)
                                               cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTg                    <  1:1064268/68‑1 (MQ=255)
                                               cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                 <  1:1482933/71‑1 (MQ=255)
                                                 tAAGTTTCTACGACTGGTATTGCGACCTGccgcccgc                                                 >  1:247908/1‑37 (MQ=16)
                                                 tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa               >  1:1842024/1‑71 (MQ=255)
                                                  aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                             <  1:824904/40‑1 (MQ=255)
                                                    gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtctc                                            >  1:824643/1‑39 (MQ=255)
                                                    gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtct                                             <  1:2590855/38‑1 (MQ=255)
                                                    gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                           >  1:81729/1‑40 (MQ=255)
                                                    gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                           >  1:2858709/1‑40 (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:2584082/40‑1 (MQ=255)
                                                        cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCaa                                       <  1:3819244/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)
                                                                    |                                                                
TCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTG  >  minE/832137‑832269

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: