Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,712 T→G 27.9% S390S (TCT→TCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,7120TG27.9% 93.7 / 52.2 68S390S (TCT→TCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (31/18);  new base G (3/16);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 8.70e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCG  >  minE/832642‑832776
                                                                      |                                                                
aGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCt                                                                  <  1:486419/71‑1 (MQ=255)
   aaaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAtt                                                               >  1:497840/1‑71 (MQ=255)
    aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAt                                                                >  1:5679711/1‑69 (MQ=255)
    aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAt                                                                >  1:5661294/1‑69 (MQ=255)
    aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTg                                                              >  1:3755832/1‑71 (MQ=255)
    aaaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTg                                                              >  1:4676397/1‑71 (MQ=255)
     aaCAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTgg                                                             >  1:309927/1‑71 (MQ=255)
     aaCAATCCGCAGTGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGAtcg                                                              >  1:1829555/1‑68 (MQ=255)
     aaCAATCCGCAGTGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGAtc                                                               >  1:1287352/1‑68 (MQ=255)
             ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCt                                                     <  1:3559558/71‑1 (MQ=255)
               agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                                   <  1:4568169/69‑1 (MQ=255)
               agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                                   <  1:5462661/69‑1 (MQ=255)
               agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                                   <  1:2090708/69‑1 (MQ=255)
               aaTGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTgg                                                   <  1:1443698/71‑1 (MQ=255)
                aTGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCt                                                                  <  1:2674713/55‑1 (MQ=255)
                 tGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTAtt                                                               <  1:3854500/57‑1 (MQ=255)
                   gAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTgg                                                             >  1:4543940/1‑57 (MQ=255)
                   gAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGc                                                >  1:3160069/1‑70 (MQ=255)
                      aaCTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg                                                  <  1:4571973/65‑1 (MQ=255)
                       aCTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                                      >  1:2670162/1‑60 (MQ=255)
                             gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCaa                                     <  1:3167222/71‑1 (MQ=255)
                             gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCa                                      >  1:3800650/1‑70 (MQ=255)
                             gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCa                                      >  1:2257604/1‑70 (MQ=255)
                                  tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                                      >  1:5345206/1‑49 (MQ=255)
                                  tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                                      >  1:4182093/1‑49 (MQ=255)
                                  tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                                      >  1:1573237/1‑49 (MQ=255)
                                  tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATgg                                 >  1:4799910/1‑70 (MQ=255)
                                     tACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCg                                               <  1:1178768/53‑1 (MQ=255)
                                     tACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATc                             <  1:4637853/71‑1 (MQ=255)
                                      aCCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCg                                               <  1:2057519/52‑1 (MQ=255)
                                        cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCtt                           <  1:1573428/70‑1 (MQ=255)
                                        cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTg                          <  1:4891240/71‑1 (MQ=255)
                                        cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTg                          <  1:1672421/71‑1 (MQ=255)
                                             gtgAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAATCt                            <  1:3070023/62‑1 (MQ=255)
                                             gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                                  >  1:2187086/1‑58 (MQ=255)
                                             gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                                  >  1:148097/1‑58 (MQ=255)
                                             gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                                  >  1:4714306/1‑58 (MQ=255)
                                               gAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTGGGGAGaaaaa                                          >  1:5324752/1‑47 (MQ=16)
                                                aaaTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg                                                  >  1:1987587/1‑39 (MQ=255)
                                                aaaTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg                                                  >  1:4203655/1‑39 (MQ=255)
                                                  aTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGaa                                           >  1:76384/1‑44 (MQ=255)
                                                  aTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGa                <  1:1441733/71‑1 (MQ=255)
                                                   tGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGAc               >  1:5576061/1‑71 (MQ=255)
                                                    ggTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGaa                                           >  1:2464821/1‑42 (MQ=255)
                                                       tGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCtt                           >  1:2226379/1‑55 (MQ=255)
                                                       tGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCtt                           >  1:1557746/1‑55 (MQ=255)
                                                        gcgcCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGagca                      <  1:1120922/59‑1 (MQ=255)
                                                         tgcCGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGa                                <  1:2862507/47‑1 (MQ=21)
                                                         cgcCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGaa                                           <  1:215322/37‑1 (MQ=255)
                                                         cgcCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCaaa         <  1:173592/71‑1 (MQ=255)
                                                         cgcCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCaa          <  1:1460577/70‑1 (MQ=255)
                                                          gcCGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCa                                      <  1:3965806/41‑1 (MQ=16)
                                                           ccGAACGGTTCTATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAATCtg                           >  1:2078621/1‑50 (MQ=37)
                                                           ccGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCa                                      <  1:2794787/40‑1 (MQ=16)
                                                           ccGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGa                                <  1:1606618/46‑1 (MQ=21)
                                                           ccGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGa                                <  1:879464/46‑1 (MQ=21)
                                                           ccGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGa                                <  1:3684995/46‑1 (MQ=21)
                                                            cGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACt      >  1:1844029/1‑71 (MQ=255)
                                                              aaCGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTgg    >  1:4364605/1‑71 (MQ=255)
                                                              aaCGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTgg    >  1:4600882/1‑71 (MQ=255)
                                                              aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATg                                  <  1:2084685/41‑1 (MQ=16)
                                                              aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATg                                  <  1:5512708/41‑1 (MQ=16)
                                                              aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAAt                              <  1:1231966/45‑1 (MQ=21)
                                                              aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAAt                              <  1:4563994/45‑1 (MQ=21)
                                                              aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAAt                              <  1:3798485/45‑1 (MQ=21)
                                                                cGGTTCTATTGGCTTCCGCTGGGGCGAGAACGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCg  <  1:61479/71‑1 (MQ=255)
                                                                cGGTTCTATTGGCTTACGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCg  <  1:20356/71‑1 (MQ=255)
                                                                   ttCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCa           >  1:2032609/1‑59 (MQ=255)
                                                                   ttCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAAc       >  1:5263683/1‑63 (MQ=255)
                                                                   ttCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGaa                               <  1:622770/39‑1 (MQ=16)
                                                                      |                                                                
AGGAAAACAATCCGGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCG  >  minE/832642‑832776

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: