Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,048 T→C 9.7% F77F (TTT→TTC 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,0480TC 9.7% 109.1 / 4.5 41F77F (TTT→TTCnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (18/19);  new base C (2/2);  total (20/21)
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.79e-01

TCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAG  >  minE/836979‑837117
                                                                     |                                                                     
tCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGttt                                                                       >  1:4065062/1‑70 (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:3799140/1‑71 (MQ=255)
              aGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGc                                                        >  1:2037232/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:1744669/1‑70 (MQ=255)
                        cAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                               <  1:1710408/70‑1 (MQ=255)
                           cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCcgc                                               <  1:3057701/67‑1 (MQ=255)
                           cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAcc                                           >  1:160268/1‑71 (MQ=255)
                           cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                            >  1:2800016/1‑70 (MQ=255)
                           cGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGAc                                            >  1:2736397/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:3776751/1‑70 (MQ=255)
                                                            gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg           >  1:1805778/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)
                                                                   tttGAACATGTGCACGGCGAATCCCGCGACCGCGGc                                      <  1:362568/36‑1 (MQ=255)
                                                                    ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGa                               >  1:3951436/1‑42 (MQ=255)
                                                                     tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGc                                 >  1:2153105/1‑39 (MQ=255)
                                                                     tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtg                          >  1:1407091/1‑46 (MQ=255)
                                                                     tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag  <  1:923048/70‑1 (MQ=255)
                                                                     |                                                                     
TCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCGACCGTGGCCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAG  >  minE/836979‑837117

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: