Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,248 T→G 14.3% S902S (TCT→TCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,2480TG14.3% 81.5 / 12.2 42S902S (TCT→TCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (26/10);  new base G (4/2);  total (30/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

AGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATG  >  minE/834179‑834317
                                                                     |                                                                     
aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                                                      >  1:378712/1‑71 (MQ=255)
aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                                                      >  1:3951602/1‑71 (MQ=255)
 ggTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAt                                                                     >  1:3084154/1‑71 (MQ=255)
    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                  <  1:19573/71‑1 (MQ=255)
    aaaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCgg                                                                  >  1:2444626/1‑71 (MQ=255)
     aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGc                                                                 >  1:1980575/1‑71 (MQ=255)
     aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGc                                                                 >  1:186387/1‑71 (MQ=255)
      aaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgggc                                                                >  1:1391442/1‑69 (MQ=255)
       aCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGccc                                                               <  1:2561972/71‑1 (MQ=255)
         cgcgCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCg                                                              >  1:1932590/1‑70 (MQ=255)
          gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                                                            >  1:1281191/1‑71 (MQ=255)
              cacacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGAt                                                         <  1:3158329/70‑1 (MQ=255)
                cacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg                                                        >  1:2937784/1‑69 (MQ=255)
                      aTGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAAt                                                >  1:320999/1‑71 (MQ=255)
                        ggTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct                                                                       <  1:3806373/46‑1 (MQ=255)
                              agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                        >  1:603530/1‑71 (MQ=255)
                              agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                        >  1:3275776/1‑71 (MQ=255)
                               gagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATgg                                       <  1:3911881/71‑1 (MQ=255)
                                agCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGc                                      >  1:1367208/1‑71 (MQ=255)
                                  cgcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCgg                                    >  1:2996914/1‑71 (MQ=255)
                                  cgcgACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATg                                                        >  1:3911127/1‑51 (MQ=255)
                                    cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt                                   >  1:558700/1‑70 (MQ=255)
                                     gATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTaa                                 <  1:1628483/71‑1 (MQ=255)
                                     gATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa                                  >  1:1968555/1‑70 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa                                                          >  1:2904507/1‑44 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa                                                          >  1:2842278/1‑44 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg                                                        >  1:102027/1‑46 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg                                                        >  1:132678/1‑46 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:306251/1‑60 (MQ=255)
                                       ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:2718844/1‑60 (MQ=255)
                                         aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa                                                          >  1:75292/1‑42 (MQ=255)
                                              gCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                        <  1:3776216/71‑1 (MQ=255)
                                               cGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGcc                       <  1:103031/71‑1 (MQ=255)
                                                  cTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:1922197/1‑49 (MQ=255)
                                                  cTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                    <  1:746364/71‑1 (MQ=255)
                                                   ttACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:2619149/1‑48 (MQ=255)
                                                        aaCGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTaaa                                <  1:3743074/53‑1 (MQ=255)
                                                        aaCGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTaaa                                <  1:890923/53‑1 (MQ=255)
                                                          cGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGAttt                         <  1:1136861/58‑2 (MQ=255)
                                                            cTTTACCTCTATCGGCCCGCTGATGGAGAATATCGGTAATGGCGGTAAAGGGATTGcc                       <  1:570014/58‑1 (MQ=255)
                                                              ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc        >  1:3054743/1‑71 (MQ=255)
                                                               tACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATACCCaggat        >  1:897088/1‑67 (MQ=255)
                                                                 ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGGCTGGaa                  >  1:3735170/1‑58 (MQ=255)
                                                                    ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATg  >  1:911210/1‑71 (MQ=255)
                                                                     |                                                                     
AGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATG  >  minE/834179‑834317

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: