Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,718 C→T 31.3% G392G (GGC→GGT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,7180CT31.3% 99.8 / 46.9 67G392G (GGC→GGTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (2/19);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 4.57e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCGAAGAAACCGAGC  >  minE/832655‑832788
                                                               |                                                                      
ggAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCt                                                                 <  1:3559558/71‑1 (MQ=255)
  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                                               <  1:2090708/69‑1 (MQ=255)
  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                                               <  1:5462661/69‑1 (MQ=255)
  agtGGAAAACCGTAGCAGTTAATACCGCAGGTGAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTgg                                                               <  1:4568169/69‑1 (MQ=255)
  aaTGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTgg                                                               <  1:1443698/71‑1 (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:2257604/1‑70 (MQ=255)
                gCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCa                                                  >  1:3800650/1‑70 (MQ=255)
                     tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                                                  >  1:5345206/1‑49 (MQ=255)
                     tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                                                  >  1:1573237/1‑49 (MQ=255)
                     tAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGc                                                                  >  1:4182093/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:1672421/71‑1 (MQ=255)
                           cAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTg                                      <  1:4891240/71‑1 (MQ=255)
                                gtgAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAATCt                                        <  1:3070023/62‑1 (MQ=255)
                                gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                                              >  1:4714306/1‑58 (MQ=255)
                                gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                                              >  1:2187086/1‑58 (MQ=255)
                                gCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATg                                              >  1:148097/1‑58 (MQ=255)
                                  gAATTGGTAGTGCCGAACGGTTCGATTGGTTTCCGCTGGGGAGaaaaa                                                      >  1:5324752/1‑47 (MQ=16)
                                   aaaTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg                                                              >  1:4203655/1‑39 (MQ=255)
                                   aaaTGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTggg                                                              >  1:1987587/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:1557746/1‑55 (MQ=255)
                                          tGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCtt                                       >  1:2226379/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:3684995/46‑1 (MQ=21)
                                              ccGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGa                                            <  1:879464/46‑1 (MQ=21)
                                              ccGAACGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGa                                            <  1:1606618/46‑1 (MQ=21)
                                               cGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACt                  >  1:1844029/1‑71 (MQ=255)
                                                 aaCGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTgg                >  1:4600882/1‑71 (MQ=255)
                                                 aaCGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTgg                >  1:4364605/1‑71 (MQ=255)
                                                 aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATg                                              <  1:2084685/41‑1 (MQ=16)
                                                 aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATg                                              <  1:5512708/41‑1 (MQ=16)
                                                 aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAAt                                          <  1:4563994/45‑1 (MQ=21)
                                                 aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAAt                                          <  1:3798485/45‑1 (MQ=21)
                                                 aaCGGTTCGATTGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAAt                                          <  1:1231966/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)
                                                           ttGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAATCTGGaa                                    <  1:1539148/41‑2 (MQ=16)
                                                           ttGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAATCTGGa                                     <  1:2229447/40‑1 (MQ=16)
                                                           ttGGTTTCCGCTGGGGAGAAAAAGGCAAATGGAATCTGGa                                     <  1:5375076/40‑1 (MQ=16)
                                                           ttGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAgcgc                              <  1:915561/47‑1 (MQ=255)
                                                           ttGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAgcgc                              <  1:4067535/47‑1 (MQ=255)
                                                             ggCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGag                                    >  1:2284561/1‑39 (MQ=255)
                                                             ggCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCGAAGAAAcc      >  1:1501415/1‑69 (MQ=255)
                                                             ggCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCGAAGAAACCg     <  1:4592070/70‑1 (MQ=255)
                                                              gCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCGAAGAAACCGAg   <  1:599192/71‑1 (MQ=255)
                                                               cTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCGAAGAAACCGAGc  <  1:2787474/71‑1 (MQ=255)
                                                               |                                                                      
GGAATGGAAAACTGTCGCCTTTAATACCAATGGCGAAATGGTTGCGCCGAACGGTTCTATTGGCTTCCGCTGGGGCGAGAAGGGCAAATGGAATCTTGAGCAGCGCGACGGCAAAACTGGCGAAGAAACCGAGC  >  minE/832655‑832788

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: