Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,294 G→T 21.2% A918S (GCC→TCC) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,2940GT21.2% 161.9 / 25.4 66A918S (GCC→TCC) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (22/30);  new base T (10/4);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

GGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGA  >  Exported/834224‑834363
                                                                      |                                                                     
ggCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTg                                                                       <  1:2770132/71‑1 (MQ=255)
    aCTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGAAAGGGATTGCCTgg                                                                  <  1:419840/71‑1 (MQ=255)
        aCGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGaa                                                                <  1:291860/70‑1 (MQ=255)
           aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg                                                                  <  1:267508/65‑1 (MQ=255)
           aaCGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                                            >  1:4518783/1‑71 (MQ=255)
              gCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTgg                                                                  <  1:3178416/62‑1 (MQ=255)
               cTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGCTAATGGCGGTAAAGGGATTGCCTGGAACACCCaga                                                        <  1:761213/71‑1 (MQ=255)
                 ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                                                                      <  1:1780629/55‑1 (MQ=255)
                 ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc                                                                      <  1:2732616/55‑1 (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:2660199/70‑1 (MQ=255)
                    ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                                    <  1:3160474/70‑1 (MQ=255)
                     ctctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga                                                  <  1:4396547/71‑1 (MQ=255)
                       ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc                                                      <  1:210801/65‑1 (MQ=255)
                        tATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATg                                                >  1:1766222/1‑70 (MQ=255)
                         aTCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGa                                              >  1:894272/1‑71 (MQ=255)
                         aTCGGCCCGCTGATGGAGAAAATAGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGa                                              >  1:4372801/1‑71 (MQ=255)
                          tCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                                                  <  1:4481599/50‑1 (MQ=255)
                          tCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGAt                                             >  1:3797728/1‑71 (MQ=255)
                           cGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATgg                                               <  1:802142/68‑1 (MQ=255)
                                cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                                                  <  1:2551366/44‑1 (MQ=255)
                                cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                                                  <  1:4663557/44‑1 (MQ=255)
                                cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                                                  <  1:2845189/44‑1 (MQ=255)
                                cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                                    <  1:3581008/58‑1 (MQ=255)
                                cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                                    <  1:5435845/58‑1 (MQ=255)
                                 gCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                                                  <  1:1579095/43‑1 (MQ=255)
                                  cTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgt                                                                  <  1:4322964/42‑2 (MQ=255)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                                                >  1:2891107/1‑42 (MQ=16)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                                                >  1:5430540/1‑42 (MQ=16)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                                                >  1:5658164/1‑42 (MQ=16)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                                                >  1:935884/1‑42 (MQ=16)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGa                                                                 >  1:198260/1‑41 (MQ=16)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGa                                                                 >  1:1536115/1‑41 (MQ=16)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                                            >  1:2966700/1‑46 (MQ=255)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                                            >  1:515111/1‑46 (MQ=255)
                                    gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                                            >  1:2030043/1‑46 (MQ=255)
                                      tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                                           <  1:2677410/45‑1 (MQ=255)
                                      tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                                           <  1:4887448/45‑1 (MQ=255)
                                      tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                                           <  1:1775653/45‑1 (MQ=255)
                                         agaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt                               >  1:5344887/1‑70 (MQ=255)
                                         agaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt                               >  1:998717/1‑70 (MQ=255)
                                          gaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCa                             <  1:4824569/71‑1 (MQ=255)
                                            aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc                           <  1:792991/71‑1 (MQ=255)
                                            aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc                           >  1:1411784/1‑71 (MQ=255)
                                            aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc                           <  1:205472/71‑1 (MQ=255)
                                             aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga                                                  >  1:3745241/1‑47 (MQ=255)
                                             aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga                                                  >  1:2346196/1‑47 (MQ=255)
                                             aaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                                    <  1:3094082/45‑1 (MQ=255)
                                              aTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTa                         <  1:5024182/71‑1 (MQ=255)
                                              aTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCAAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTa                         >  1:665690/1‑71 (MQ=255)
                                                 ggTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTacac                      >  1:3989957/1‑71 (MQ=255)
                                                     aTGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCgaga                                                  >  1:4239199/1‑39 (MQ=255)
                                                     aTGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCa                    >  1:2032368/1‑69 (MQ=255)
                                                     aTGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCa                    >  1:3679147/1‑69 (MQ=255)
                                                       ggCGGTAAAGGGATTTCGTGGAATACCCAGGATGAAATCGATTTCCTCGGCAAACTCAATTACACCaagc                 <  1:1823331/70‑3 (MQ=255)
                                                       ggCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTaca                       >  1:1385535/1‑64 (MQ=255)
                                                        gCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGc                <  1:3435614/70‑1 (MQ=255)
                                                             aaaGGGATTGCCTGGAACACCCAGAGCGAGAGGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAgg          >  1:4798302/1‑71 (MQ=255)
                                                              aaGGGATTGCCTGGAACACCCAGAGCGAGATGGAt                                             >  1:1149982/1‑35 (MQ=255)
                                                                gggATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTcc       <  1:3989496/71‑1 (MQ=255)
                                                                gggATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTcc       >  1:761562/1‑71 (MQ=255)
                                                                   aTTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGaa            >  1:2290267/1‑63 (MQ=255)
                                                                     tGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAg           >  1:1857751/1‑62 (MQ=255)
                                                                     tGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAg           >  1:2595004/1‑62 (MQ=255)
                                                                     tGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCg   >  1:1114380/1‑70 (MQ=255)
                                                                      gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa  <  1:1887641/70‑1 (MQ=255)
                                                                      gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa  <  1:827078/70‑1 (MQ=255)
                                                                      gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa  <  1:5102534/70‑1 (MQ=255)
                                                                      |                                                                     
GGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGA  >  Exported/834224‑834363

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: