Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,287,091 C→T 31.6% Q126* (CAA→TAA)  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
*W3110S.gb1,287,0910CT31.6% 50.3 / 30.1 38Q126* (CAA→TAA) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base C (12/14);  new base T (9/3);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTG  >  W3110S.gb/1287024‑1287156
                                                                   |                                                                 
gCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGTaa                                                                 <  1:1423997/70‑1 (MQ=255)
gCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCaa                                                                 <  1:2758806/70‑1 (MQ=255)
    ttAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCaa                                                                 <  1:2604898/66‑1 (MQ=255)
     tAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCg                                                             >  1:1030850/1‑69 (MQ=255)
      aaaCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAg                                                          >  1:1735985/1‑71 (MQ=255)
            ccGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTg                                                     >  1:737404/1‑70 (MQ=255)
             cgcgAATTGCCGGACCATCTGCCGCTGTAACTGGAGTACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTgg                                                    >  1:2075302/1‑70 (MQ=255)
                   ttGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGttt                                             >  1:4632783/1‑71 (MQ=255)
                     gCCGGACCATCTGCCGCTGTATCTGGAGTATCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTg                                            >  1:1064095/1‑70 (MQ=255)
                       cGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGaaa                                         <  1:4236711/71‑1 (MQ=255)
                         gACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGa                                       <  1:2716320/71‑1 (MQ=255)
                          aCCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGat                                      <  1:2734829/71‑1 (MQ=255)
                                   cGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCg                                                       >  1:2169593/1‑45 (MQ=255)
                                       gTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGa                                       >  1:3174802/1‑57 (MQ=255)
                                        tATCTGGAGTACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTc                         <  1:1274052/70‑1 (MQ=255)
                                         aTCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTg                       >  1:1684623/1‑71 (MQ=255)
                                            tGGATTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGcc                              >  1:702443/1‑61 (MQ=255)
                                             ggAGTACCTGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGAtt                          >  1:1054357/1‑64 (MQ=255)
                                              gAGTACCTGGCGCAGCTGCCGCAAAGTGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTgg                      <  1:4173479/67‑1 (MQ=255)
                                              gAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTgg                      <  1:3680028/67‑1 (MQ=255)
                                              gAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTgg                      <  1:3987280/67‑1 (MQ=255)
                                              gAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTgg                      <  1:3993282/67‑1 (MQ=255)
                                              gAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGc                     <  1:857712/68‑1 (MQ=255)
                                                gTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGc                >  1:2363425/1‑71 (MQ=255)
                                                    cTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTgg                      <  1:2166163/61‑1 (MQ=255)
                                                     tGGCGCAGCTGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAgc           >  1:1952014/1‑71 (MQ=255)
                                                      ggCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGatat                                    >  1:4950852/1‑45 (MQ=255)
                                                         gcAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAgcgcgc       <  1:499200/71‑1 (MQ=255)
                                                              tGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTg  >  1:2186482/1‑71 (MQ=255)
                                                              tGCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTg  >  1:4861081/1‑71 (MQ=255)
                                                               gCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTc                         >  1:832481/1‑47 (MQ=255)
                                                               gCCGTAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTc                         >  1:1724560/1‑47 (MQ=255)
                                                               gCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCAtt                  >  1:2475510/1‑54 (MQ=255)
                                                               gCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGa             >  1:4724691/1‑59 (MQ=255)
                                                               gCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGa             >  1:2616757/1‑59 (MQ=255)
                                                                 cGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCAt                   <  1:61477/51‑1 (MQ=255)
                                                                 cGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCAt                   <  1:537460/51‑1 (MQ=255)
                                                                 cGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCAt                   <  1:4796252/51‑1 (MQ=255)
                                                                   |                                                                 
GCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTG  >  W3110S.gb/1287024‑1287156

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: