Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 834,854 C→A 24.1% D1104E (GAC→GAA narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported834,8540CA24.1% 181.1 / 39.4 83D1104E (GAC→GAAnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (28/35);  new base A (8/12);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAG  >  Exported/834785‑834921
                                                                     |                                                                   
cAACCAGGAAAAAGCGCTCAACTTTCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa                                                                    <  1:3852092/71‑1 (MQ=255)
 aaCCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa                                                                    <  1:2695136/70‑1 (MQ=255)
    ccggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAcc                                                                 <  1:4368715/68‑1 (MQ=255)
    cAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACaa                                                                   >  1:1402299/1‑68 (MQ=255)
     cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                                                  <  1:567457/67‑1 (MQ=255)
     cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                                                  <  1:3157486/67‑1 (MQ=255)
     cggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAc                                                                  <  1:189225/67‑1 (MQ=255)
     aGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctg                                                               >  1:3144394/1‑71 (MQ=255)
       gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                                             <  1:916451/71‑1 (MQ=255)
       gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                                             <  1:1125150/71‑1 (MQ=255)
       gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct                                                             <  1:3603640/71‑1 (MQ=255)
       gAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTAGAGTGAAAACct                                                                >  1:2600314/1‑68 (MQ=255)
         aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctga                                                           >  1:4887880/1‑71 (MQ=255)
          aaaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaat                                                          >  1:4634689/1‑71 (MQ=255)
           aaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatg                                                         <  1:3896675/71‑1 (MQ=255)
           aaGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatg                                                         <  1:5477393/71‑1 (MQ=255)
            aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc                                                        >  1:305566/1‑71 (MQ=255)
             gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                                       <  1:4230928/71‑1 (MQ=255)
             gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                                       <  1:778149/71‑1 (MQ=255)
             gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                                       <  1:3647857/71‑1 (MQ=255)
              cgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                                       >  1:3873271/1‑70 (MQ=255)
                cttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAc                                                    >  1:3654334/4‑71 (MQ=255)
                cttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAc                                                    >  1:132462/4‑71 (MQ=255)
                 ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACcttct                                                             <  1:2322660/59‑1 (MQ=255)
                 ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg                                                   >  1:2673952/3‑70 (MQ=255)
                 ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg                                                   >  1:2037702/3‑70 (MQ=255)
                   aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGCCtt                                                  >  1:456030/1‑70 (MQ=255)
                   aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                                 <  1:5208427/71‑1 (MQ=255)
                   aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                                 <  1:2071783/71‑1 (MQ=255)
                   aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                                 <  1:3100376/71‑1 (MQ=255)
                   aaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                                                 >  1:270191/1‑68 (MQ=255)
                    aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt                                                  >  1:1076048/1‑69 (MQ=255)
                    aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt                                                  >  1:1695595/1‑69 (MQ=255)
                     cTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAg                                               >  1:5416809/1‑71 (MQ=255)
                       tCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGGCa                                                                    <  1:2015197/48‑1 (MQ=255)
                       tCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGt                                             >  1:3679281/1‑71 (MQ=255)
                         cTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTcg                                           <  1:1511231/71‑1 (MQ=255)
                           gACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgctaa                                                           <  1:3006596/53‑3 (MQ=255)
                            aCGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg                                        <  1:4046506/71‑1 (MQ=255)
                            aCGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                                                 <  1:1805862/62‑4 (MQ=255)
                             cgccgcACCAGAAATGGGGCATTCAGTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                                                 <  1:1920594/61‑4 (MQ=255)
                              gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct                                                             <  1:1119767/48‑1 (MQ=19)
                              gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct                                                             <  1:1606080/48‑1 (MQ=19)
                               ccgcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtgg                                     <  1:4735454/71‑1 (MQ=255)
                                cgcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggt                                    >  1:2559416/1‑71 (MQ=255)
                                 gcACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtc                                   >  1:4439869/1‑71 (MQ=255)
                                   aCCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccg                                 <  1:368621/71‑1 (MQ=255)
                                     cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                                         <  1:5671764/45‑1 (MQ=255)
                                     cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                                         <  1:3101428/45‑1 (MQ=255)
                                     cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                                         <  1:434018/45‑1 (MQ=255)
                                      aGAAGTAGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcgg                                        <  1:3939067/61‑1 (MQ=255)
                                       gAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAgg                                              >  1:2987148/1‑54 (MQ=255)
                                          gTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCt                          >  1:1825984/1‑71 (MQ=255)
                                           tGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                                       <  1:5218110/41‑1 (MQ=255)
                                                tATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtc                                   <  1:4076261/56‑1 (MQ=255)
                                                 aTCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGtt                      <  1:825208/68‑1 (MQ=255)
                                                   ccactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccg                                 <  1:3825328/55‑1 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                                >  1:4400425/1‑39 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                                >  1:5580994/1‑39 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                                >  1:4990428/1‑39 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                                >  1:4188080/1‑39 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                                >  1:3318768/1‑39 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTcg                                           >  1:5123155/1‑44 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg                 <  1:4499858/70‑1 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg                 >  1:4805053/1‑70 (MQ=255)
                                                    cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg                 >  1:2545110/1‑70 (MQ=255)
                                                       tccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGt                                             >  1:1228432/1‑39 (MQ=255)
                                                         cacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGcc            <  1:4441411/70‑1 (MQ=255)
                                                         cacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGcc            <  1:5009583/70‑1 (MQ=255)
                                                         cacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGcc            <  1:3818974/70‑1 (MQ=255)
                                                          acCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGccca          <  1:1165771/71‑3 (MQ=255)
                                                          acCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGa          <  1:2058386/71‑1 (MQ=255)
                                                           ccTACAGCGACAACCTGCTGATGCTGACTTTAGGTCTcggt                                       <  1:5643508/41‑1 (MQ=255)
                                                           ccTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGAt         >  1:1094528/1‑71 (MQ=255)
                                                           ccTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGAt         >  1:3040779/1‑71 (MQ=255)
                                                               cAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccg                                 >  1:3555616/1‑43 (MQ=255)
                                                               cAGCGACAACCTGCTGATGCTGACTTTAGGTCGcggtggtccg                                 >  1:3593418/1‑43 (MQ=255)
                                                                aGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTg                     <  1:1167006/54‑1 (MQ=255)
                                                                aGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTg                     <  1:4804939/54‑1 (MQ=255)
                                                                 gCGACACCCTGCTGATGCTGACTTTAGGTCGCGGGGGTCCGGTGGTCTGGTTg                     <  1:815115/53‑1 (MQ=255)
                                                                 gCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCa     >  1:5554046/1‑69 (MQ=255)
                                                                 gCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCa     >  1:2693890/1‑69 (MQ=255)
                                                                  cGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAg  <  1:1988310/71‑1 (MQ=255)
                                                                     |                                                                   
CAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAG  >  Exported/834785‑834921

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: