Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 835,820 T→C 31.3% L180L (TTG→CTG)  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,8200TC31.3% 118.8 / 75.3 98L180L (TTG→CTG) narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (24/43);  new base C (21/9);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 2.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 101 total aligned reads displayed.

AAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAT  >  Exported/835751‑835884
                                                                     |                                                                 
aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                                                  >  1:44376/1‑71 (MQ=255)
aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt                                                                  <  1:3531289/71‑1 (MQ=255)
  aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                                                >  1:2123703/1‑71 (MQ=255)
  aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc                                                                >  1:4756403/1‑71 (MQ=255)
        caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCgcgc                                                           >  1:594367/2‑70 (MQ=255)
        caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc                                                          >  1:3923764/2‑71 (MQ=255)
           cTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc                                                             <  1:4191217/65‑1 (MQ=255)
             tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgtg                                                      >  1:773307/1‑70 (MQ=255)
             tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                                                     >  1:2789212/1‑71 (MQ=255)
              cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc                                                           <  1:3201256/64‑1 (MQ=255)
              cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc                                                     >  1:631081/1‑70 (MQ=255)
                acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGag                                                  <  1:686382/71‑2 (MQ=255)
                acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                                                  <  1:1707990/71‑1 (MQ=255)
                acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa                                                  <  1:5528516/71‑1 (MQ=255)
                  aacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAaca                                                <  1:4207407/71‑1 (MQ=255)
                        caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt                                          >  1:432562/4‑71 (MQ=255)
                        caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt                                          >  1:2853531/4‑71 (MQ=255)
                        cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc                                           <  1:1291312/70‑1 (MQ=255)
                        cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                                          >  1:5462064/1‑71 (MQ=255)
                        cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                                          <  1:5149149/71‑1 (MQ=255)
                         aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                                          >  1:4348567/1‑70 (MQ=255)
                         aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                         <  1:1663632/71‑1 (MQ=255)
                         aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                         <  1:1781324/71‑1 (MQ=255)
                         aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                         <  1:3186099/71‑1 (MQ=255)
                          aaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCtca                                        <  1:3666103/70‑3 (MQ=255)
                            aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat                                      >  1:3079316/1‑70 (MQ=255)
                            aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat                                      >  1:3935238/1‑70 (MQ=255)
                              gagATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAAtcc                                   <  1:5372620/69‑1 (MQ=255)
                              gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAccc                                   <  1:560082/71‑1 (MQ=255)
                              gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc                                    <  1:691822/70‑1 (MQ=255)
                               cGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt                                          >  1:3961053/1‑64 (MQ=255)
                                gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                                                        >  1:2120516/1‑49 (MQ=255)
                                gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                                 <  1:2563743/71‑1 (MQ=255)
                                gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                                 <  1:5643594/71‑1 (MQ=255)
                                gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                                 <  1:5176742/71‑1 (MQ=255)
                                gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg                                 <  1:1687524/71‑1 (MQ=255)
                                    tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                                                        <  1:698975/45‑1 (MQ=255)
                                       aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt                          >  1:1436486/1‑71 (MQ=255)
                                           aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc                       >  1:2746859/1‑70 (MQ=255)
                                            gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGa                     >  1:4997169/1‑71 (MQ=255)
                                            gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa                               <  1:2105302/61‑1 (MQ=255)
                                            gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa                               <  1:370947/61‑1 (MQ=255)
                                            gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                      <  1:4585217/70‑1 (MQ=255)
                                            gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                      <  1:4811165/70‑1 (MQ=255)
                                            gTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCg                      <  1:1012600/70‑1 (MQ=255)
                                             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                                                        <  1:1620059/36‑1 (MQ=255)
                                             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCtg                                                        <  1:1471726/36‑1 (MQ=255)
                                             ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGc                                <  1:4370896/59‑1 (MQ=255)
                                              ttgAAAACAACTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGAcc                   >  1:3261578/3‑71 (MQ=255)
                                              tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACt                                              >  1:2540103/1‑45 (MQ=255)
                                              tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                                       <  1:5537493/52‑1 (MQ=255)
                                              tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                                       <  1:149834/52‑1 (MQ=255)
                                              tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGaa                                       <  1:3212023/52‑1 (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)
                                               tgAAAACACCTTCATGATGTACATGCCGCGCCTGTGCGAACACTGCCTCaa                                       >  1:1353259/2‑51 (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:4604916/1‑71 (MQ=255)
                                               cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt                  >  1:83831/1‑71 (MQ=255)
                                                gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc                                    >  1:2313407/1‑52 (MQ=255)
                                                gAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc                                    >  1:4653343/1‑52 (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)
                                                  aaaCACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaa                                       >  1:1943985/1‑48 (MQ=38)
                                                  aaaCACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaa                                       >  1:5518181/1‑48 (MQ=38)
                                                   aaCACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGc                >  1:3584047/1‑69 (MQ=255)
                                                   aaCACCTTCATGATGTATTTGCCGCGCCTCTGCGAACACTGCCt                                          >  1:4406981/1‑44 (MQ=255)
                                                    acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCg             <  1:3126183/71‑1 (MQ=255)
                                                    acacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCg             <  1:1157043/71‑1 (MQ=255)
                                                    acacCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGccca             <  1:3879190/71‑2 (MQ=255)
                                                    acacCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGccca             <  1:3838918/71‑2 (MQ=255)
                                                     cacCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGccc              >  1:2062836/1‑69 (MQ=255)
                                                      acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg           <  1:2239291/71‑1 (MQ=255)
                                                      acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg           <  1:2812730/71‑1 (MQ=255)
                                                      acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg           <  1:3395560/71‑1 (MQ=255)
                                                      acCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAg           <  1:1122157/71‑1 (MQ=255)
                                                       ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                         <  1:5573785/41‑1 (MQ=255)
                                                       ccTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg                                         <  1:1034980/41‑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)
                                                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaa                                       >  1:3309811/1‑39 (MQ=25)
                                                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCggc       <  1:4020327/70‑2 (MQ=255)
                                                           cATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCgg        <  1:860058/69‑1 (MQ=255)
                                                            atgatgTACCTGCCGCGCCTGTGCGAACACTGCCTCaa                                       >  1:2073067/1‑38 (MQ=25)
                                                            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:5454100/71‑1 (MQ=255)
                                                               atgTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAt  <  1:4593942/71‑1 (MQ=255)
                                                               atgTACCTGCCGCGCCTGTGCGAACACTGCCTCAATCCCAGCTGCGTGGCGACCTGCCCAAGCGgcgcta   <  1:1889141/70‑3 (MQ=255)
                                                                     |                                                                 
AAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCTGCCCGAGCGGTGCGAT  >  Exported/835751‑835884

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: