Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 837,087 G→A 9.8% 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,0870GA 9.8% 118.1 / 3.9 41A90A (GCG→GCAnarJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base G (18/18);  new base A (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.51e-01

CCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCAT  >  minE/837024‑837156
                                                               |                                                                     
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:595993/70‑1 (MQ=255)
            aCTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGgcg                                                     <  1:3956713/70‑1 (MQ=255)
               gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg                                                  >  1:3776751/1‑70 (MQ=255)
               gTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAg                                                  >  1:1805778/1‑70 (MQ=255)
                ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                                                                   <  1:242214/52‑1 (MQ=255)
                ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                                                                   <  1:2545857/52‑1 (MQ=255)
                ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGAcc                                                             <  1:1417426/58‑1 (MQ=255)
                ttGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTa                                                >  1:3485711/1‑71 (MQ=255)
                       ttGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGa                                                                      >  1:3951436/1‑42 (MQ=255)
                        tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtggtg                                                                 >  1:1407091/1‑46 (MQ=255)
                        tGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagcag                                         <  1:923048/70‑1 (MQ=255)
                           aCATGTGCACGGCGAATCCCGCGACCGCGTCCAGGCGATGGTGGACCTGCTGGCGCAGTACGagca                                          >  1:3913568/1‑66 (MQ=255)
                           aCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                                                                   >  1:4076199/1‑41 (MQ=255)
                           aCATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGAtgg                                                                   >  1:3273629/1‑41 (MQ=255)
                             aTGTGCACGGCGAATCCCGCGACCTCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGc                                   <  1:3654186/71‑1 (MQ=255)
                              tgtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCt                                  >  1:3275405/1‑71 (MQ=255)
                               gtgCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCtt                                 >  1:3032360/1‑71 (MQ=255)
                                      gCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTa                          <  1:2712548/71‑1 (MQ=255)
                                       cGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACg                                              <  1:4042621/50‑1 (MQ=255)
                                       cGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAg                             >  1:1611202/1‑67 (MQ=255)
                                       cGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTaa                         <  1:2757249/71‑1 (MQ=255)
                                          aTCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACa                      <  1:1775296/71‑1 (MQ=255)
                                           tCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAg                     >  1:4025480/1‑71 (MQ=255)
                                            cccGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGc                    >  1:1330406/1‑71 (MQ=255)
                                             ccGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGcc                   >  1:3590066/1‑71 (MQ=255)
                                               gcgACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCcgc                 >  1:513325/1‑71 (MQ=255)
                                                 gACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTCAGTTAAACAGCCGCGaa              <  1:1569580/71‑1 (MQ=255)
                                                  aCCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGaa              >  1:2019091/1‑71 (MQ=255)
                                                        gCCAGGCTATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCg        <  1:3581183/71‑1 (MQ=255)
                                                        gCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCg        <  1:1094242/71‑1 (MQ=255)
                                                             gCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAg                             <  1:32513/45‑1 (MQ=255)
                                                             gCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGAcc    <  1:996935/70‑1 (MQ=255)
                                                              cGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCa   >  1:1275008/1‑70 (MQ=255)
                                                              cGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCAt  >  1:1270231/1‑71 (MQ=255)
                                                               |                                                                     
CCGCGCCACGTCACTGTTGCTGTTTGAACATGTGCACGGCGAATCCCGCGACCGCGGCCAGGCGATGGTGGACCTGCTGGCGCAGTACGAGCAGCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCAT  >  minE/837024‑837156

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: