Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,858 T→C 41.2% C192C (TGT→TGC narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported835,8580TC41.2% 40.3 / 94.1 85C192C (TGT→TGCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base C (22/13);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

AGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAT  >  Exported/835790‑835927
                                                                     |                                                                     
aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt                                                                     >  1:1436486/1‑71 (MQ=255)
    aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc                                                                  >  1:2746859/1‑70 (MQ=255)
     gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGa                                                                >  1:4997169/1‑71 (MQ=255)
     gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                                                                 <  1:1012600/70‑1 (MQ=255)
     gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                                                                 <  1:4811165/70‑1 (MQ=255)
     gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                                                                 <  1:4585217/70‑1 (MQ=255)
       ttgAAAACAACTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGAcc                                                              >  1:3261578/3‑71 (MQ=255)
       tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                                                                      >  1:2505592/1‑63 (MQ=255)
        tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGCACACTGCCTCAATCCCAGCTGCGTGGCGACCt                                                             >  1:2191295/2‑71 (MQ=255)
        tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCt                                                             >  1:664389/2‑71 (MQ=255)
        tgAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCt                                                             >  1:727956/2‑71 (MQ=255)
        cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg                                                                      >  1:5626662/1‑62 (MQ=255)
        cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgt                                                                       <  1:2833587/61‑1 (MQ=255)
        cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt                                                             >  1:139118/1‑71 (MQ=255)
        cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt                                                             >  1:83831/1‑71 (MQ=255)
        cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt                                                             >  1:4604916/1‑71 (MQ=255)
         gAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTTCCTCAATCCCAGCTGCGTGGCGACCTg                                                            >  1:2362692/1‑71 (MQ=255)
         gAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTg                                                            >  1:4134574/1‑71 (MQ=255)
          aaaaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGc                                                           <  1:2847320/71‑1 (MQ=255)
          aaaaCACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGc                                                           >  1:559860/1‑71 (MQ=255)
            aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGc                                                           >  1:3584047/1‑69 (MQ=255)
             acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCg                                                        <  1:3126183/71‑1 (MQ=255)
             acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCg                                                        <  1:1157043/71‑1 (MQ=255)
             acacCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGccca                                                        <  1:3838918/71‑2 (MQ=255)
             acacCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGccca                                                        <  1:3879190/71‑2 (MQ=255)
              cacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGccc                                                         >  1:2062836/1‑69 (MQ=255)
               acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                      <  1:2812730/71‑1 (MQ=255)
               acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                      <  1:2239291/71‑1 (MQ=255)
               acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                      <  1:3395560/71‑1 (MQ=255)
               acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                      <  1:1122157/71‑1 (MQ=255)
                ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg                                                      <  1:5009832/70‑1 (MQ=255)
                ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGc                                                     >  1:1982968/1‑71 (MQ=255)
                   tCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGt                                                  >  1:4422061/1‑71 (MQ=255)
                   tCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggc                                                  >  1:1641382/1‑70 (MQ=255)
                    cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggc                                                  <  1:4020327/70‑2 (MQ=255)
                    cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCgg                                                   <  1:860058/69‑1 (MQ=255)
                     atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCgg                                                   <  1:73443/68‑1 (MQ=255)
                     atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgc                                                >  1:3207135/1‑71 (MQ=255)
                      tgatgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCg                                               <  1:5125726/71‑1 (MQ=255)
                      tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgct                                               <  1:5492776/71‑2 (MQ=255)
                      tgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCACAGCTGCGTGGCGACCTGCCCAAGCGgcgc                                                >  1:3424952/1‑70 (MQ=255)
                        atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                             <  1:4593942/71‑1 (MQ=255)
                        atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt                                             <  1:5454100/71‑1 (MQ=255)
                        atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgcta                                              <  1:1889141/70‑3 (MQ=255)
                             cctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTAc                                         >  1:2959457/3‑70 (MQ=255)
                             cctGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTAc                                         >  1:965612/3‑70 (MQ=255)
                               tGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACATGCCCGAGCGGTGCGATTTACAAg                                      >  1:4295982/1‑71 (MQ=255)
                                gCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaac                                     >  1:5474769/1‑69 (MQ=255)
                                gCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaac                                     >  1:1914269/1‑69 (MQ=255)
                                 ccGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg                                    >  1:4131440/1‑71 (MQ=255)
                                 ccGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGCGACCTGCCCGAGCGGTGCGATTTACAAGCGt                                   >  1:4347390/1‑71 (MQ=255)
                                   gcgcCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc                                                                  >  1:3401348/1‑39 (MQ=255)
                                   gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                                       >  1:3909760/1‑66 (MQ=255)
                                   gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                                       >  1:2057523/1‑66 (MQ=255)
                                   gcgcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaa                                       >  1:4792923/1‑66 (MQ=255)
                                     gcCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTga                                 <  1:1424948/70‑1 (MQ=255)
                                     gcCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaa                                <  1:4008683/71‑1 (MQ=255)
                                       cTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgct                                               <  1:1982515/54‑2 (MQ=255)
                                        tgtgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCt                                                             >  1:863524/1‑39 (MQ=255)
                                        tgtgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaag                               <  1:5114094/69‑1 (MQ=255)
                                         gtgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaa                                >  1:3075944/1‑67 (MQ=255)
                                          tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTAc                                         <  1:2743499/57‑1 (MQ=255)
                                          tgCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                             <  1:5674081/69‑1 (MQ=255)
                                          tgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCg                                                    >  1:2189811/1‑46 (MQ=12)
                                          tgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa                             <  1:5420271/69‑1 (MQ=255)
                                          tgCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCgaagaa                             <  1:1690358/69‑1 (MQ=255)
                                           gCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCg                                               >  1:2820836/1‑50 (MQ=255)
                                           gCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCAGCGGTGCg                                               >  1:1934623/1‑49 (MQ=255)
                                           gCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaa                                      >  1:298998/1‑58 (MQ=255)
                                           gCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACaaa                                      >  1:4783294/1‑58 (MQ=255)
                                            cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCg                                    <  1:4920128/60‑1 (MQ=255)
                                            cGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTg                                  >  1:1109417/1‑62 (MQ=255)
                                             gAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAagaaga                           >  1:5406577/1‑68 (MQ=255)
                                              aaCACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATgg                        >  1:2365420/1‑70 (MQ=255)
                                              aaCACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGGCGCTATCTACAAACGCGAagaagacg                         <  1:809713/69‑3 (MQ=255)
                                                  cTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                             <  1:4544790/61‑1 (MQ=255)
                                                  cTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCg                   <  1:3433784/71‑1 (MQ=255)
                                                       tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATTGCATCGTCCt               >  1:854048/1‑70 (MQ=255)
                                                       tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATc                    >  1:1163641/1‑65 (MQ=255)
                                                       tGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCt               >  1:5067144/1‑70 (MQ=255)
                                                         aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaag                               >  1:3137770/1‑52 (MQ=255)
                                                         aaCCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGGGAAGAAGATGGCATCGTCCTGAt            <  1:4859008/71‑1 (MQ=255)
                                                            cccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCg          <  1:4446315/71‑1 (MQ=255)
                                                             ccGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                             <  1:998695/51‑1 (MQ=255)
                                                              cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                            <  1:3814880/51‑1 (MQ=255)
                                                              cGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAc        <  1:2858466/71‑1 (MQ=255)
                                                               ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                            <  1:2387615/50‑1 (MQ=255)
                                                               ggCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGAcc       <  1:4683953/71‑1 (MQ=255)
                                                                 cATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGaagaag                            <  1:4593058/48‑1 (MQ=255)
                                                                   tgtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                             <  1:2313650/45‑1 (MQ=255)
                                                                    gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                             <  1:543383/44‑1 (MQ=255)
                                                                    gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                             <  1:5629327/44‑1 (MQ=255)
                                                                    gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                             <  1:1390108/44‑1 (MQ=255)
                                                                    gtgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTgaagaa                             <  1:1287185/44‑1 (MQ=255)
                                                                    gtgtgGCGACCTGCCCGAGCGGGGCGATTTACAAGCGTgaagaa                             <  1:1740031/44‑1 (MQ=255)
                                                                     tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt  >  1:4396349/1‑70 (MQ=255)
                                                                     tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt  >  1:2016248/1‑70 (MQ=255)
                                                                     tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt  >  1:1775170/1‑70 (MQ=255)
                                                                     tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt  >  1:534635/1‑70 (MQ=255)
                                                                     tgtgGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAt  >  1:456002/1‑70 (MQ=255)
                                                                     |                                                                     
AGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGAT  >  Exported/835790‑835927

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