Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,040 T→C 23.9% L75L (TTG→CTG)  narJ → molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE837,0400TC23.9% 112.8 / 35.4 67L75L (TTG→CTG) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (20/31);  new base C (8/8);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAG  >  minE/836979‑837108
                                                             |                                                                    
tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                                              >  1:1243331/1‑70 (MQ=255)
tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                             <  1:2925923/71‑1 (MQ=255)
tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                             <  1:3830050/71‑1 (MQ=255)
  cGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                           <  1:2634176/71‑1 (MQ=255)
  cGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                           <  1:4486640/71‑1 (MQ=255)
   gCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                           <  1:3021745/70‑1 (MQ=255)
    cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTg                                                                    <  1:3115401/60‑1 (MQ=255)
    cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTg                                                                    <  1:3730991/60‑1 (MQ=255)
    cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                           <  1:679740/69‑1 (MQ=255)
    cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa                                                         <  1:5285151/71‑1 (MQ=255)
       gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgt                                                      >  1:4743326/1‑71 (MQ=255)
           cccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCa                                                   <  1:2142867/70‑1 (MQ=255)
           cccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCa                                                   <  1:5129755/70‑1 (MQ=255)
            ccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACg                                                 >  1:4304606/1‑71 (MQ=255)
            ccAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACg                                                 <  1:3982208/71‑1 (MQ=255)
             cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCAc                                                  >  1:1774486/1‑69 (MQ=255)
             cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCAc                                                  >  1:528943/1‑69 (MQ=255)
             cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACgg                                                >  1:4123467/1‑71 (MQ=255)
             cAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACgg                                                >  1:5328383/1‑71 (MQ=255)
               ggCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                              <  1:1142602/71‑1 (MQ=255)
                gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                              >  1:1429854/1‑70 (MQ=255)
                gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                              >  1:1968832/1‑70 (MQ=255)
                   cAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATc                                          <  1:1273041/71‑1 (MQ=255)
                   cAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATc                                          <  1:3038748/71‑1 (MQ=255)
                      tACAGCGAACTGTTCGACCGTGGCCGCTCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcg                                       >  1:5331585/1‑71 (MQ=255)
                       aCAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                      >  1:4351798/1‑71 (MQ=255)
                        cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                      <  1:3180810/70‑1 (MQ=255)
                        cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                      >  1:3275029/1‑70 (MQ=255)
                        cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                      <  1:4655449/70‑1 (MQ=255)
                        cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                      >  1:1954183/1‑70 (MQ=255)
                          gCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                             <  1:5411941/45‑1 (MQ=255)
                            gAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                      <  1:3267845/66‑1 (MQ=255)
                             aaCTGTTCGACCGTGGCCGCGCCACGTCACTGTTg                                                                    <  1:4796508/35‑1 (MQ=255)
                                 gTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGt                                                                <  1:2861544/35‑1 (MQ=255)
                                  ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGa                                                            >  1:2390564/1‑38 (MQ=255)
                                  ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGa                                    >  1:3726861/1‑62 (MQ=255)
                                  ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGa                                    >  1:1313003/1‑62 (MQ=255)
                                     gACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                        <  1:347441/71‑1 (MQ=255)
                                     gACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAACGCGGCCAGGc                        <  1:5569837/71‑1 (MQ=255)
                                      aCCGCGGGCGCACCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGc                        <  1:248931/70‑1 (MQ=255)
                                       ccGCGGGCGCACCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGcaa                      <  1:2033718/71‑3 (MQ=255)
                                       ccGCGGGCGCACCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGcaa                      <  1:4219983/71‑3 (MQ=255)
                                         gTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                    <  1:3691372/71‑1 (MQ=255)
                                            gCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggt                  <  1:1477056/70‑1 (MQ=255)
                                             gcgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg                <  1:2844976/70‑1 (MQ=255)
                                             gcgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg                <  1:4929275/70‑1 (MQ=255)
                                              cgcgCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                              >  1:4555090/1‑40 (MQ=255)
                                                  ccACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct            >  1:5264352/1‑70 (MQ=255)
                                                   cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtgttgg                <  1:5638980/65‑1 (MQ=255)
                                                   cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg                <  1:4664664/65‑1 (MQ=255)
                                                   cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtgg                   >  1:3892172/1‑62 (MQ=255)
                                                   cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc          >  1:4472078/1‑71 (MQ=255)
                                                   cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc          >  1:4174436/1‑71 (MQ=255)
                                                   cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc          >  1:3472542/1‑71 (MQ=255)
                                                   cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc          >  1:2444782/1‑71 (MQ=255)
                                                   cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtgg                <  1:3683966/65‑1 (MQ=255)
                                                    aCGTCCCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGcaa                      >  1:3736215/1‑56 (MQ=255)
                                                     cGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgctg        >  1:3029830/1‑71 (MQ=255)
                                                      gTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGc                        <  1:5464509/54‑1 (MQ=255)
                                                       tcgctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGgc      <  1:1388791/68‑1 (MQ=255)
                                                        cgctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGgcg     >  1:4414028/3‑71 (MQ=255)
                                                         gctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtg                    >  1:3958707/2‑55 (MQ=255)
                                                         gctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGgcg     <  1:1021282/69‑1 (MQ=255)
                                                         gctgctgctgTTCGAACATATTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGgcg     <  1:3803419/69‑1 (MQ=255)
                                                         aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                        <  1:900444/51‑1 (MQ=255)
                                                         aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg     >  1:4511865/1‑70 (MQ=255)
                                                          cTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg     >  1:955244/1‑69 (MQ=255)
                                                           tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtgg                <  1:352813/57‑1 (MQ=255)
                                                           tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtgg                <  1:4333687/57‑1 (MQ=255)
                                                           tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg  >  1:2359348/1‑71 (MQ=255)
                                                            gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                        <  1:2335640/48‑1 (MQ=255)
                                                            gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                        <  1:973012/48‑1 (MQ=255)
                                                             ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                    >  1:1548587/1‑51 (MQ=255)
                                                             |                                                                    
TCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAG  >  minE/836979‑837108

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: