Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,233 C→T 20.9% Y897Y (TAC→TAT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,2330CT20.9% 120.5 / 25.9 62Y897Y (TAC→TATnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (4/9);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

AGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAA  >  Exported/834163‑834301
                                                                      |                                                                    
aGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc                                                                      <  1:1073281/71‑1 (MQ=255)
aGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc                                                                      <  1:4876532/71‑1 (MQ=255)
 gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc                                                                      >  1:1805810/1‑70 (MQ=255)
      atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc                                                                  <  1:1171900/67‑1 (MQ=255)
      atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc                                                                  <  1:1531459/67‑1 (MQ=255)
      atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc                                                                  <  1:1759713/67‑1 (MQ=255)
         tGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAAc                                                                  <  1:3983531/66‑1 (MQ=255)
           aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTATAGCGCGATTATCCGGCGACTTACGAACGCtt                                                              <  1:1439717/68‑1 (MQ=255)
           aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCtt                                                              <  1:1762731/68‑1 (MQ=255)
           aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                                           >  1:3137129/1‑71 (MQ=255)
           aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                                           >  1:1664951/1‑71 (MQ=255)
               cAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct                                                       <  1:2961830/71‑1 (MQ=255)
                aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                                      >  1:4131261/1‑71 (MQ=255)
                aGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGc                                                                <  1:2274748/61‑1 (MQ=255)
                aGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACctcgc                                                      <  1:1793489/71‑3 (MQ=255)
                 ggTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa                                                                   >  1:4401786/1‑57 (MQ=255)
                  gTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCttt                                                             <  1:3205612/62‑1 (MQ=255)
                    aaaaCCTCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                  >  1:3373877/1‑71 (MQ=255)
                    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                  >  1:2028724/1‑71 (MQ=255)
                    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                  <  1:537943/71‑1 (MQ=255)
                    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                  >  1:3264816/1‑71 (MQ=255)
                    aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg                                                  >  1:4059875/1‑71 (MQ=255)
                          gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                                           >  1:2475642/1‑56 (MQ=255)
                          gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                                           >  1:2682170/1‑56 (MQ=255)
                          gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt                                            <  1:151129/71‑1 (MQ=255)
                              cacacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                                           >  1:1043999/1‑52 (MQ=255)
                                cacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctc                                                        <  1:5323215/53‑1 (MQ=255)
                                 acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct                                                       <  1:38355/53‑1 (MQ=255)
                                 acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg                                        >  1:2989700/1‑68 (MQ=255)
                                  caTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa                                                      >  1:3775983/1‑53 (MQ=255)
                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa                                                                   >  1:4049950/1‑39 (MQ=255)
                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa                                                                   >  1:2688210/1‑39 (MQ=255)
                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa                                                                   >  1:2117296/1‑39 (MQ=255)
                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGa                                                                    >  1:1524131/1‑38 (MQ=255)
                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAAc                                                                  <  1:4409047/40‑1 (MQ=255)
                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                                           >  1:4276509/1‑47 (MQ=255)
                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                                           >  1:1337985/1‑47 (MQ=255)
                                   aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc                                                           >  1:4700031/1‑47 (MQ=255)
                                    ttATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGaaa                                  <  1:276779/71‑1 (MQ=255)
                                       tGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAt                                                     >  1:1417611/1‑49 (MQ=255)
                                         gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                             <  1:2863443/71‑1 (MQ=255)
                                         gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg                             <  1:512139/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:2331647/71‑1 (MQ=255)
                                               gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg                       <  1:2770130/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)
                                                                      |                                                                    
AGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAA  >  Exported/834163‑834301

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: