Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,087 G→A 20.4% A90A (GCG→GCA 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,0870GA20.4% 101.2 / 20.2 49A90A (GCG→GCAnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base G (13/26);  new base A (5/5);  total (18/31)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGAC  >  minE/837020‑837153
                                                                   |                                                                  
gTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                                                 <  1:2108400/71‑1 (MQ=255)
  gggcgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggt                                                               <  1:874864/68‑1 (MQ=255)
     cgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGa                                                            <  1:2391791/71‑1 (MQ=255)
        accACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct                                                         >  1:3603164/2‑71 (MQ=255)
        accACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct                                                         >  1:1722484/2‑71 (MQ=255)
          cACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc                                                       >  1:1785592/1‑71 (MQ=255)
          cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                                                       >  1:2225904/1‑71 (MQ=255)
          cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                                                       >  1:2221230/1‑71 (MQ=255)
          cACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgc                                                       <  1:1465211/71‑1 (MQ=255)
            cGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtg                                                              <  1:1637635/62‑1 (MQ=255)
                gctgctgctgTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggt                                                               >  1:1731759/2‑57 (MQ=255)
                aCTGTTGCTGTTTGAACATGTGCACGGCTAATCCCGCGACCGCGGCCAGGCGATGGTGGACctgctg                                                     <  1:2037465/67‑1 (MQ=255)
                  tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggt                                                               <  1:2800032/55‑1 (MQ=255)
                  tGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGAcc                                                          >  1:1267926/1‑60 (MQ=255)
                    ctGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTa                                             <  1:29106/70‑1 (MQ=255)
                       cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                                                 >  1:1634932/1‑48 (MQ=255)
                       cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtg                                                                 >  1:1804337/1‑48 (MQ=255)
                       cTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGa                                          <  1:3321884/71‑1 (MQ=255)
                       cTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAAtggtgg                                                             <  1:537710/52‑1 (MQ=255)
                         gTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                                         >  1:3108591/1‑70 (MQ=255)
                          tttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                                         <  1:3216655/69‑1 (MQ=255)
                          ttcgAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGa                                          >  1:3425150/4‑68 (MQ=255)
                           ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGag                                         <  1:2047565/68‑1 (MQ=255)
                            tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag                                      >  1:2668388/1‑70 (MQ=255)
                            tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGAGCAGTACAagcag                                      >  1:326620/1‑70 (MQ=255)
                             gAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcagc                                     <  1:2181591/70‑1 (MQ=255)
                                cATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACgg                                 <  1:930769/71‑1 (MQ=255)
                                  tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt                               <  1:2562968/71‑1 (MQ=255)
                                  tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt                               >  1:1476639/1‑71 (MQ=255)
                                       aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa                                                            <  1:277858/37‑1 (MQ=255)
                                       aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGa                                                            <  1:1546293/37‑1 (MQ=255)
                                       aCGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCa                           <  1:2165574/70‑1 (MQ=255)
                                          ccGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTGGCGGAGTATGAAAAGGTCGGCCTGCAgct                        <  1:1233270/69‑3 (MQ=255)
                                           cGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaa                      <  1:2922964/71‑1 (MQ=255)
                                            gAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAg                          <  1:145029/66‑1 (MQ=255)
                                            gAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaaa                     <  1:3263003/71‑1 (MQ=255)
                                             aaTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaaa                     <  1:1018712/70‑1 (MQ=255)
                                              aTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagc                                        >  1:2055620/1‑50 (MQ=255)
                                              aTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACa                   <  1:2831789/71‑1 (MQ=255)
                                              aTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACa                   <  1:406679/71‑1 (MQ=255)
                                               tCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaaa                     >  1:2930074/1‑68 (MQ=255)
                                                 ccGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAg                  <  1:3461646/69‑1 (MQ=255)
                                                    cgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGc                 >  1:1244298/1‑67 (MQ=255)
                                                    cgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcgcg             <  1:581286/71‑1 (MQ=255)
                                                      aCCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCATTTAAACAGCCGCGa            <  1:2929158/70‑1 (MQ=255)
                                                          cgGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGc       <  1:3897475/71‑1 (MQ=255)
                                                           ggCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGcc      >  1:2983916/1‑71 (MQ=255)
                                                               aGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaa                      <  1:1308766/51‑1 (MQ=255)
                                                               aGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGAc  <  1:2483109/71‑1 (MQ=255)
                                                                   |                                                                  
GTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGAC  >  minE/837020‑837153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: