Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 832,457 A→C 6.7% A305A (GCA→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported832,4570AC 6.7% 238.2 / 9.5 90A305A (GCA→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (33/51);  new base C (3/3);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

ACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCA  >  Exported/832388‑832524
                                                                     |                                                                   
acCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCaa                                                                    <  1:3964059/71‑1 (MQ=255)
  cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                                  >  1:4531944/1‑71 (MQ=255)
  cAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATg                                                                  >  1:5291713/1‑71 (MQ=255)
   aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                 >  1:1044184/1‑71 (MQ=255)
   aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                 <  1:4082993/71‑1 (MQ=255)
   aGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                 <  1:3169226/71‑1 (MQ=255)
    gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                                <  1:2195511/71‑1 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                               >  1:3295111/1‑71 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                               >  1:3504627/1‑71 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                               >  1:2937726/1‑71 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                               >  1:4369683/1‑71 (MQ=255)
     aCTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                                >  1:5640309/1‑70 (MQ=255)
       tACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                             >  1:1171134/1‑71 (MQ=255)
         cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                                           >  1:4875849/1‑71 (MQ=255)
         cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                                           <  1:2331287/71‑1 (MQ=255)
          gCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                                          >  1:4841478/1‑71 (MQ=255)
           cTGAATTCGCCAAAGTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCg                                                         <  1:2510370/71‑1 (MQ=255)
            tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                 <  1:2759535/62‑1 (MQ=255)
            tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCt                                                             <  1:471326/66‑1 (MQ=255)
            tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                        >  1:4293305/1‑71 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgcg                                                               >  1:4746651/1‑62 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                        <  1:30606/69‑1 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAt                                                       <  1:2656352/70‑1 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg                                                      <  1:4550437/71‑1 (MQ=255)
              aaaTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg                                                      >  1:4135535/1‑71 (MQ=255)
                aTCGCCAAACTGTGCGATCTGTGGCTGGCCCCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggt                                                    >  1:1041489/1‑70 (MQ=255)
                aTCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggg                                                    >  1:1459041/1‑71 (MQ=255)
                   gCCAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCa                                                 >  1:5162204/1‑71 (MQ=255)
                    ccAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCat                                                <  1:4956967/71‑2 (MQ=255)
                     cAAATTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCatg                                               >  1:972398/1‑69 (MQ=255)
                      aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt                                              <  1:4691403/71‑1 (MQ=255)
                      aaaTTGTGCGACCAGTGGCTGGCACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATGGGCCATGt                                              <  1:3292348/71‑1 (MQ=255)
                      aaaCTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                              <  1:647448/71‑1 (MQ=255)
                         cTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                 >  1:726344/1‑49 (MQ=255)
                         cTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATgg                                                                 >  1:234138/1‑49 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                            >  1:711005/1‑49 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                            >  1:4667873/1‑49 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                            >  1:1951876/1‑49 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTg                                                            >  1:1306486/1‑49 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                                          <  1:3870600/51‑1 (MQ=255)
                              cGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGc                                                          <  1:950805/51‑1 (MQ=255)
                                 tCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGt                                    >  1:1595232/1‑70 (MQ=255)
                                  cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                                >  1:4393556/1‑41 (MQ=255)
                                  cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGgc                                                                >  1:2328187/1‑41 (MQ=255)
                                  cTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGa                                                        >  1:5161605/1‑49 (MQ=255)
                                     tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                   <  1:4393558/67‑1 (MQ=255)
                                     tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                                  <  1:3415584/68‑1 (MQ=255)
                                      ggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                   <  1:5047364/66‑1 (MQ=255)
                                       gctggcACCGAAACAAGGCACTGATAGCGCCCTGGCGATGGCAATg                                                      >  1:3183114/1‑46 (MQ=255)
                                            cACCGAAACAGGGCACCGATGCGGCAATGGCGCTgg                                                           <  1:5106141/36‑1 (MQ=255)
                                             aCCGAAACAGGGCTCCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                                  <  1:545220/60‑1 (MQ=255)
                                              ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                              <  1:3193421/47‑1 (MQ=255)
                                              ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                              <  1:2098150/47‑1 (MQ=255)
                                              ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGt                                              <  1:689026/47‑1 (MQ=255)
                                              ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                        <  1:2653088/69‑1 (MQ=255)
                                               cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGa                      >  1:991128/1‑70 (MQ=255)
                                               cGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                     <  1:2245183/71‑1 (MQ=255)
                                                gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtg                                          >  1:1505257/1‑49 (MQ=255)
                                                gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                     >  1:2577609/1‑70 (MQ=255)
                                                  aaCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACaa                   >  1:5637108/1‑70 (MQ=255)
                                                   aCAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGa                                  <  1:4273547/54‑1 (MQ=255)
                                                    cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                     <  1:2510651/66‑1 (MQ=255)
                                                    cAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAAccc                <  1:3708554/71‑1 (MQ=255)
                                                     aGGGCACCGATGCGTCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                     <  1:3677352/65‑1 (MQ=255)
                                                       ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                      <  1:5364858/46‑1 (MQ=255)
                                                       ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                      <  1:2955774/46‑1 (MQ=255)
                                                       ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                      <  1:3481310/46‑1 (MQ=255)
                                                       ggCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAAtgctgc                                      <  1:4930043/46‑1 (MQ=255)
                                                        gCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGc            <  1:903663/71‑1 (MQ=255)
                                                         cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                   <  1:3703196/47‑1 (MQ=255)
                                                         cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                   <  1:3710010/47‑1 (MQ=255)
                                                         cACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                   <  1:249938/47‑1 (MQ=255)
                                                           ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                   <  1:1822337/45‑1 (MQ=255)
                                                           ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTg                                   <  1:116302/45‑1 (MQ=255)
                                                           ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg         <  1:860809/71‑1 (MQ=255)
                                                           ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg         <  1:1412590/71‑1 (MQ=255)
                                                           ccGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGATGCGTGAATTCCACCTCGACAACCCAAGCCAg         <  1:1297774/71‑1 (MQ=255)
                                                            cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg         >  1:2784953/1‑70 (MQ=255)
                                                            cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt        <  1:4327533/71‑1 (MQ=255)
                                                            cGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGt        <  1:5093007/71‑1 (MQ=255)
                                                             gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                        <  1:4033357/54‑1 (MQ=255)
                                                             gATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAg         <  1:5112555/69‑1 (MQ=255)
                                                              aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt      <  1:1572162/71‑1 (MQ=255)
                                                              aTGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt      <  1:481577/71‑1 (MQ=255)
                                                              aTGCGGCAATGGCGCTGGCGATGGGCCACCTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAt      <  1:2553456/71‑1 (MQ=255)
                                                                  ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                        <  1:2450968/49‑1 (MQ=255)
                                                                  ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTc                        <  1:4976215/49‑1 (MQ=255)
                                                                  ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa  <  1:400012/71‑1 (MQ=255)
                                                                  ggCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa  <  1:4598708/71‑1 (MQ=255)
                                                                    cAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTAttt    >  1:4337798/1‑67 (MQ=255)
                                                                    cAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCAAGTAttt    >  1:4649254/1‑67 (MQ=255)
                                                                     aaTGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGAc                     <  1:868249/49‑1 (MQ=255)
                                                                     |                                                                   
ACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCA  >  Exported/832388‑832524

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: