Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,599 A→T 17.8% E1019D (GAA→GAT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,5990AT17.8% 144.0 / 25.3 67E1019D (GAA→GATnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base T (6/6);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGT  >  Exported/834534‑834668
                                                                 |                                                                     
cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACg                                                                   <  1:1001419/70‑1 (MQ=255)
cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGt                                                                  >  1:4343050/1‑71 (MQ=255)
cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGt                                                                  >  1:4414918/1‑71 (MQ=255)
   gATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTc                                                               <  1:1515092/71‑1 (MQ=255)
    atatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCt                                                              <  1:4720274/71‑1 (MQ=255)
       ttCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa                                                            >  1:3980592/1‑70 (MQ=255)
         cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACaa                                                         >  1:913106/1‑71 (MQ=255)
         cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                                          <  1:3276936/70‑1 (MQ=255)
         cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                                          <  1:19916/70‑1 (MQ=255)
             cacaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGcc                                                     >  1:5485810/1‑71 (MQ=255)
              gcaGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGcgg                                                    <  1:1426495/70‑3 (MQ=255)
              acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                                    >  1:3431265/1‑71 (MQ=255)
              acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                                    >  1:775879/1‑71 (MQ=255)
              acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                                    >  1:2810182/1‑71 (MQ=255)
                aGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGt                                                  <  1:3199465/71‑1 (MQ=255)
                aGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGt                                                  <  1:4438733/71‑1 (MQ=255)
                  ccGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa                                                                       >  1:314518/1‑48 (MQ=255)
                  ccGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa                                                                       >  1:1951803/1‑48 (MQ=255)
                      gcAAAATTATCTCCAGCCCGCCCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac                                             >  1:1316838/1‑70 (MQ=255)
                      gcAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac                                             >  1:4136311/1‑70 (MQ=255)
                       taaaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACAcc                                            >  1:2529333/2‑70 (MQ=255)
                       cAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                                          >  1:1576026/1‑56 (MQ=255)
                        aaaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCggg                                                   >  1:2979045/1‑62 (MQ=255)
                         aaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAAc                                         >  1:692219/1‑71 (MQ=255)
                          aaTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACg                                        >  1:5147618/1‑71 (MQ=255)
                            ttATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGtt                                      <  1:4566101/71‑1 (MQ=255)
                               tctcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCAc                                   >  1:422689/1‑71 (MQ=255)
                               tctcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCAc                                   >  1:3109015/1‑71 (MQ=255)
                                   cAGCCCGACCTGGTCTGGTCTGGAAGATGAACATGTTTCTTACAACGCCGGTTACACCAACGTTCAc                                   <  1:3874263/67‑1 (MQ=255)
                                   cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAga                               <  1:3325402/71‑2 (MQ=255)
                                   cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc                               >  1:2279472/1‑71 (MQ=255)
                                   cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc                               >  1:3771642/1‑71 (MQ=255)
                                   cAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAgt                               >  1:3245122/1‑70 (MQ=255)
                                    aGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAgtt                              <  1:3154329/71‑3 (MQ=255)
                                        cGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa                                 >  1:5564691/1‑64 (MQ=255)
                                          aCCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc                         >  1:3053275/1‑70 (MQ=255)
                                           ccTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                       <  1:2859978/71‑1 (MQ=255)
                                            ctggtccggtCTGGAAAGCTATCACGTCTCCTATAACGCGGGATACACCAACGtt                                      <  1:1773061/55‑1 (MQ=255)
                                            ctggtccggtCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGtt                                      <  1:4670662/55‑1 (MQ=255)
                                             tggtctggtcTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATg                     >  1:4345503/1‑71 (MQ=255)
                                               gtctggtctgGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTg                             <  1:4221275/61‑1 (MQ=255)
                                               gtctggtctgGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGc                   >  1:1666562/1‑71 (MQ=255)
                                               gtccggtctgGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTgg                    <  1:4443603/70‑1 (MQ=255)
                                                 ctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg                  >  1:1296049/1‑70 (MQ=255)
                                                 ccggtctggAAAGAGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGc                   >  1:2991480/3‑69 (MQ=255)
                                                  tggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                       <  1:1256267/64‑1 (MQ=255)
                                                  tggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                       <  1:1249383/64‑1 (MQ=255)
                                                  tggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                       <  1:1082990/64‑1 (MQ=255)
                                                  tggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                       <  1:997608/64‑1 (MQ=255)
                                                   ggtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCAc               <  1:4005560/71‑1 (MQ=255)
                                                    gtctggAAGATGAACACGTTTCTTACAACGCCGGTTaca                                              >  1:4015644/1‑39 (MQ=255)
                                                    gtctggAAGATGAACACGTTTCTTACAACGCCGGTTaca                                              >  1:3717928/1‑39 (MQ=255)
                                                    gtctggAAGATGAACACGTTTCTTACAACGCCGGTTa                                                >  1:4715456/1‑37 (MQ=255)
                                                    gtctggAAGATGAACACGTTTCTTACAACGCCGGTTa                                                >  1:1806213/1‑37 (MQ=255)
                                                    gtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGt                 <  1:191628/68‑1 (MQ=255)
                                                    gtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTAc               <  1:3967120/70‑1 (MQ=255)
                                                    gtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTAc               <  1:4367756/70‑1 (MQ=255)
                                                    gtctggAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACg              >  1:4241321/1‑71 (MQ=255)
                                                     tctggATGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTg                             <  1:3257601/55‑1 (MQ=255)
                                                      ctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGct            >  1:2064287/1‑71 (MQ=255)
                                                       tggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGc             <  1:4808438/69‑1 (MQ=255)
                                                         gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctctc         >  1:5173350/1‑71 (MQ=255)
                                                         gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctctc         <  1:1764458/71‑1 (MQ=255)
                                                         gAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctctc         <  1:1157887/71‑1 (MQ=255)
                                                           aGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTg       <  1:1013436/71‑1 (MQ=255)
                                                            gATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc                         <  1:5523829/52‑1 (MQ=255)
                                                              tGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGctcac         >  1:223304/1‑64 (MQ=255)
                                                               gAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc                         <  1:4209999/49‑1 (MQ=255)
                                                               gAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcg   >  1:1469327/1‑71 (MQ=255)
                                                                aaCACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt  <  1:2843237/71‑1 (MQ=255)
                                                                 acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                       <  1:2962516/49‑1 (MQ=255)
                                                                 acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa                       <  1:2550197/49‑1 (MQ=255)
                                                                 acacGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGgtcgt  >  1:4333227/1‑70 (MQ=255)
                                                                 |                                                                     
CGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGT  >  Exported/834534‑834668

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: