Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,190 C→T 11.9% S216S (AGC→AGT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1900CT11.9% 231.1 / 16.0 84S216S (AGC→AGTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base T (6/4);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGG  >  minE/832122‑832260
                                                                    |                                                                      
gcgaTGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGttt                                                                      >  1:3752547/4‑71 (MQ=255)
gCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                                                                      >  1:2286859/1‑71 (MQ=255)
 cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                                                                      >  1:3806575/1‑70 (MQ=255)
      tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                                                                <  1:1302671/71‑1 (MQ=255)
         aTGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTg                                                             <  1:1520559/71‑1 (MQ=255)
         aTGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGAc                                                               <  1:2846263/69‑1 (MQ=255)
               tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                                                        <  1:1480589/70‑1 (MQ=255)
               tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg                                                       >  1:2595552/1‑71 (MQ=255)
                 ttACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                                                                     <  1:585165/55‑1 (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:1830429/71‑1 (MQ=255)
                    cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt                                                  <  1:2150949/71‑1 (MQ=255)
                     gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                 >  1:917353/1‑71 (MQ=255)
                     gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                 >  1:1258687/1‑71 (MQ=255)
                     gCACCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                                                                  >  1:3469226/1‑54 (MQ=255)
                      cATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa                                                          <  1:1785476/61‑1 (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:2526776/71‑1 (MQ=255)
                          gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc                                            <  1:3345869/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)
                               cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                                                                >  1:3525635/1‑46 (MQ=255)
                               cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt                                                 >  1:1317014/1‑61 (MQ=255)
                               cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg                                       <  1:2010019/71‑1 (MQ=255)
                               cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg                                       >  1:1408826/1‑71 (MQ=255)
                               cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg                                       >  1:2843832/1‑71 (MQ=255)
                               cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg                                       >  1:2483430/1‑71 (MQ=255)
                                 cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                                                                      >  1:89341/1‑38 (MQ=255)
                                 cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCt                                                                       >  1:1523585/1‑37 (MQ=255)
                                 cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                     <  1:1525087/71‑1 (MQ=255)
                                 cGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGtc                                     >  1:1970190/1‑71 (MQ=255)
                                 cGCTATCTCTAGCTGATTGGCGGTACTTGCTTAAGCtt                                                                      >  1:3912387/1‑38 (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:2645165/70‑1 (MQ=255)
                                     aTCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGcc                                  >  1:3292544/1‑70 (MQ=255)
                                     aTCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGAc                                                               <  1:1050555/41‑1 (MQ=255)
                                      tctgtcGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGcc                                  <  1:3875853/69‑1 (MQ=255)
                                       ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                >  1:1034909/1‑70 (MQ=255)
                                       ctctcGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGc                                >  1:1599583/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:3029009/71‑1 (MQ=255)
                                                 ttGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGc                     >  1:2810613/1‑71 (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:2858709/1‑40 (MQ=255)
                                                                   gCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTcc                                  >  1:81729/1‑40 (MQ=255)
                                                                    cTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCTCAAACCTGGGGCGAGCAAACTGACGTACCgg  <  1:375734/71‑1 (MQ=255)
                                                                    cTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTg                         <  1:1304482/48‑1 (MQ=255)
                                                                    |                                                                      
GCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGG  >  minE/832122‑832260

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: