Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,040 T→C 10.0% 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,0400TC10.0% 107.6 / 4.5 40L75L (TTG→CTG) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base C (2/2);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

GATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTA  >  minE/836977‑837110
                                                               |                                                                      
gATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGt                                                                  <  1:2362405/70‑1 (MQ=255)
  tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                                                >  1:4065062/1‑70 (MQ=255)
  tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGt                                                                  >  1:4127202/1‑68 (MQ=255)
  tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                               <  1:80472/71‑1 (MQ=255)
   ccGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTg                                                               >  1:590696/1‑70 (MQ=255)
    cGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGaa                                                             >  1:149827/1‑71 (MQ=255)
     gCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAAc                                                            <  1:2005756/71‑1 (MQ=255)
      cTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACa                                                           >  1:3852910/1‑71 (MQ=255)
         gATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAt                                                          <  1:2311282/69‑1 (MQ=255)
          aTGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgtg                                                       <  1:422146/71‑1 (MQ=255)
          aTGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtgtg                                                       <  1:1875212/71‑1 (MQ=255)
                aGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc                                                 >  1:2037232/1‑71 (MQ=255)
                aGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc                                                 >  1:3799140/1‑71 (MQ=255)
                  gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                <  1:3288558/70‑1 (MQ=255)
                  gcgcAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                >  1:954371/1‑70 (MQ=255)
                        tACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCg                                                <  1:519187/64‑1 (MQ=255)
                        tACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATccc                                          >  1:3704445/1‑70 (MQ=255)
                        tACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcg                                         >  1:1823996/1‑71 (MQ=255)
                          cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACAtg                                                         >  1:2546614/1‑53 (MQ=255)
                          cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                        <  1:1710408/70‑1 (MQ=255)
                          cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                        >  1:1744669/1‑70 (MQ=255)
                             cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                        <  1:3057701/67‑1 (MQ=255)
                             cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc                                    >  1:160268/1‑71 (MQ=255)
                             cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                     >  1:2736397/1‑70 (MQ=255)
                             cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                     >  1:2800016/1‑70 (MQ=255)
                                aCTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                     >  1:1340696/1‑67 (MQ=255)
                                aCTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATCTACACGGCGAATCCCGCGAc                                     >  1:528105/1‑67 (MQ=255)
                                    ttCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAAc                                                            <  1:1029069/40‑1 (MQ=255)
                                          cGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgcg                                       <  1:404767/55‑1 (MQ=255)
                                               ccGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtgg                  <  1:3914957/71‑1 (MQ=255)
                                                cgcaCCACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGa                 <  1:23359/71‑1 (MQ=255)
                                                   gcCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACct              <  1:2252290/71‑1 (MQ=255)
                                                   accACGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACct              >  1:857270/2‑71 (MQ=255)
                                                      aCGTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACctgc            >  1:46547/1‑70 (MQ=255)
                                                       cGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGCACctgctg          <  1:3545037/71‑1 (MQ=255)
                                                        gTCGCTGCTGCTGTTCGAACATGTTCATGCCGAGTCCCGCGATCGCGGCCAGGCAATGGTGGACCTGCTgg         <  1:3787014/71‑1 (MQ=255)
                                                           aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg       <  1:3956713/70‑1 (MQ=255)
                                                           aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg       <  1:595993/70‑1 (MQ=255)
                                                              gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg    >  1:1805778/1‑70 (MQ=255)
                                                              gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg    >  1:3776751/1‑70 (MQ=255)
                                                               ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                     <  1:2545857/52‑1 (MQ=255)
                                                               ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                     <  1:242214/52‑1 (MQ=255)
                                                               ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGAcc               <  1:1417426/58‑1 (MQ=255)
                                                               ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTa  >  1:3485711/1‑71 (MQ=255)
                                                               |                                                                      
GATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTA  >  minE/836977‑837110

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: