Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,248 T→G 27.6% S902S (TCT→TCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,2480TG27.6% 61.0 / 42.9 58S902S (TCT→TCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (17/25);  new base G (5/11);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

CAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATG  >  Exported/834178‑834317
                                                                      |                                                                     
cAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct                                                                       <  1:2961830/71‑1 (MQ=255)
 aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                                                      >  1:4131261/1‑71 (MQ=255)
     aaaaCCTCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                  >  1:3373877/1‑71 (MQ=255)
     aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                  <  1:537943/71‑1 (MQ=255)
     aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                  >  1:2028724/1‑71 (MQ=255)
     aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                  >  1:4059875/1‑71 (MQ=255)
     aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                                  >  1:3264816/1‑71 (MQ=255)
           gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                                                            <  1:151129/71‑1 (MQ=255)
                  acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct                                                                       <  1:38355/53‑1 (MQ=255)
                  acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg                                                        >  1:2989700/1‑68 (MQ=255)
                   caTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                                                      >  1:3775983/1‑53 (MQ=255)
                     ttATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGaaa                                                  <  1:276779/71‑1 (MQ=255)
                        tGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAt                                                                     >  1:1417611/1‑49 (MQ=255)
                          gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                                             <  1:512139/71‑1 (MQ=255)
                          gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                                             <  1:2863443/71‑1 (MQ=255)
                           tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGGGAAAATCGGt                                            <  1:918393/71‑1 (MQ=255)
                           tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                            <  1:2213556/71‑1 (MQ=255)
                           tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                            <  1:4618116/71‑1 (MQ=255)
                             gTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:2423121/1‑71 (MQ=255)
                               agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                        >  1:1718139/1‑71 (MQ=255)
                               agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg                                        >  1:2292926/1‑71 (MQ=255)
                                gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg                                       <  1:2770130/71‑1 (MQ=255)
                                gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg                                       <  1:2331647/71‑1 (MQ=255)
                                 agCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGc                                      <  1:2732854/71‑1 (MQ=255)
                                   cgcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAAt                                                <  1:2180099/59‑1 (MQ=255)
                                     cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt                                   >  1:2191912/1‑70 (MQ=255)
                                       aTTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                                             <  1:3378964/58‑1 (MQ=255)
                                        ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGc                                                                 >  1:3829941/1‑37 (MQ=255)
                                        ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:5053776/1‑60 (MQ=255)
                                        ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt                                   >  1:2175296/1‑67 (MQ=255)
                                        ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAg                               >  1:4055910/1‑71 (MQ=255)
                                          aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                                                            <  1:4876335/40‑1 (MQ=255)
                                          aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAgg                              <  1:3738083/70‑1 (MQ=255)
                                             cGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAgg                              <  1:4685043/67‑1 (MQ=255)
                                              ggCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTg                         <  1:2770132/71‑1 (MQ=255)
                                                  aCTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGAAAGGGATTGCCTgg                    <  1:419840/71‑1 (MQ=255)
                                                    ttACAAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAggg                             <  1:3928344/61‑1 (MQ=255)
                                                      aCGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGaa                  <  1:291860/70‑1 (MQ=255)
                                                       tgAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAAt                                         <  1:271263/45‑1 (MQ=255)
                                                         aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg                    <  1:267508/65‑1 (MQ=255)
                                                         aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc              >  1:4518783/1‑71 (MQ=255)
                                                           cGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                            <  1:3049915/39‑1 (MQ=255)
                                                           cGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt                                            <  1:1646611/39‑1 (MQ=255)
                                                           cGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAggg                             <  1:584154/54‑1 (MQ=255)
                                                            gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGa                            <  1:2556674/54‑1 (MQ=255)
                                                            gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGa                            <  1:2753986/54‑1 (MQ=255)
                                                            gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg                    <  1:3178416/62‑1 (MQ=255)
                                                             cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:4420579/1‑39 (MQ=255)
                                                             cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:3443090/1‑39 (MQ=255)
                                                             cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa                                          >  1:2996780/1‑39 (MQ=255)
                                                             cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGCTAATGGCGGTAAAGGGATTGCCTGGAACACCCaga          <  1:761213/71‑1 (MQ=255)
                                                              tttACCTCTATCGGCCCGCGGATGGAGAAAATCGGTa                                           >  1:2821064/1‑37 (MQ=255)
                                                               ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                        <  1:1780629/55‑1 (MQ=255)
                                                               ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                        <  1:2732616/55‑1 (MQ=255)
                                                                tACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa                                  >  1:2957626/1‑44 (MQ=255)
                                                                tACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc        <  1:3115838/70‑1 (MQ=255)
                                                                tACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATACCCaggat        <  1:4847850/70‑4 (MQ=255)
                                                                 aCCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCg       <  1:5329022/70‑1 (MQ=255)
                                                                  ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga      <  1:3160474/70‑1 (MQ=255)
                                                                  ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga      <  1:2660199/70‑1 (MQ=255)
                                                                   ctctATCGGCCCGCTGATGGAGAAAATCGGTAATgg                                       <  1:2364127/36‑1 (MQ=255)
                                                                   ctctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga    <  1:4396547/71‑1 (MQ=255)
                                                                     ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc        <  1:210801/65‑1 (MQ=255)
                                                                      tATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATg  >  1:1766222/1‑70 (MQ=255)
                                                                      |                                                                     
CAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATG  >  Exported/834178‑834317

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: