Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,296 C→G 20.9% A918A (GCC→GCG) ‡ narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,2960CG20.9% 139.9 / 30.0 67A918A (GCC→GCG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (22/31);  new base G (10/4);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 7.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

ACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGA  >  Exported/834228‑834363
                                                                    |                                                                   
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)
              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:5435845/58‑1 (MQ=255)
                            cGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCga                                                    <  1:3581008/58‑1 (MQ=255)
                             gCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg                                                                  <  1:1579095/43‑1 (MQ=255)
                              cTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgt                                                                  <  1:4322964/42‑2 (MQ=255)
                                gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                                                >  1:5430540/1‑42 (MQ=16)
                                gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGaa                                                                >  1:2891107/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:515111/1‑46 (MQ=255)
                                gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                                            >  1:2966700/1‑46 (MQ=255)
                                gATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATAcc                                                            >  1:2030043/1‑46 (MQ=255)
                                  tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                                           <  1:4887448/45‑1 (MQ=255)
                                  tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                                           <  1:2677410/45‑1 (MQ=255)
                                  tGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACAccc                                                           <  1:1775653/45‑1 (MQ=255)
                                     agaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt                               >  1:998717/1‑70 (MQ=255)
                                     agaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCt                               >  1:5344887/1‑70 (MQ=255)
                                      gaAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCa                             <  1:4824569/71‑1 (MQ=255)
                                        aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc                           <  1:205472/71‑1 (MQ=255)
                                        aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc                           >  1:1411784/1‑71 (MQ=255)
                                        aaaTCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAAc                           <  1:792991/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:3679147/1‑69 (MQ=255)
                                                 aTGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCa                    >  1:2032368/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:2595004/1‑62 (MQ=255)
                                                                 tGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAg           >  1:1857751/1‑62 (MQ=255)
                                                                 tGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCg   >  1:1114380/1‑70 (MQ=255)
                                                                  gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa  <  1:5102534/70‑1 (MQ=255)
                                                                  gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa  <  1:1887641/70‑1 (MQ=255)
                                                                  gCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGa  <  1:827078/70‑1 (MQ=255)
                                                                    |                                                                   
ACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATGGATCTGCTGCGTAAGCTCAACTACACCAAAGCGGAAGGTCCGGCGA  >  Exported/834228‑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: