Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,250,824 T→C H640R (CAC→CGC)  napA ← nitrate reductase, periplasmic, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,250,8240TC99.5% 820.9 / ‑5.0 200H640R (CAC→CGC) napAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (91/108);  total (91/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 268 total aligned reads displayed.

TTACCGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGTGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGT  >  REL606/2250752‑2250899
                                                                        |                                                                           
ttACCGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATcgtgaccgcgac                                                                          <  2:706578/76‑1 (MQ=255)
ttACCGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGGGCCCGATcgtgaacgcgac                                                                          >  1:3460828/1‑76 (MQ=255)
 tACCGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTTTGGTAGTCATCGAACGGGGGCCGATCgtgaccccgacc                                                                         >  2:3219257/1‑76 (MQ=255)
 tACCGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCgtgaccgcgacc                                                                         <  2:7770117/76‑1 (MQ=255)
  aCCGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGtgaccgcgaccg                                                                        <  1:3618249/76‑1 (MQ=255)
   ccGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCCTCGAACGGTGCCAGGTCGtgaccgcgagcgg                                                                       >  1:7284274/1‑75 (MQ=255)
    cGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGaccgcgaccgcg                                                                      <  1:8294748/76‑1 (MQ=255)
     gTTCACCACCGGCCAGCGCAGCCCGCGCGCTTTGTGGTAGTCATCGCACGGTGCCAGAACCTGGCCCCGGCCGCGa                                                                     >  2:2372998/1‑76 (MQ=255)
     gTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGaccgcgaccgcgc                                                                     >  1:6657756/1‑75 (MQ=255)
     gTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGAccgcgaccgcga                                                                     >  1:1839755/1‑76 (MQ=255)
      ttCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGGGCCAGATCGTGACcgcgaccgcgac                                                                    >  1:9623546/1‑76 (MQ=255)
      ttCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGGCATCGAAAGGGGCCAGATCGTGACcgcgaccgcgac                                                                    >  2:6295974/1‑76 (MQ=255)
        caccacCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCa                                                                  <  2:5291667/76‑1 (MQ=255)
        caccacCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCa                                                                  <  1:5934167/76‑1 (MQ=255)
         cccaccGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCaa                                                                 <  1:9913432/75‑1 (MQ=255)
          ccaccGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCaaa                                                                <  1:4690841/76‑1 (MQ=255)
            accGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAAcc                                                              <  2:7364467/76‑1 (MQ=255)
            accGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAAcc                                                              <  2:8852462/76‑1 (MQ=255)
             ccGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCa                                                             <  1:2578812/76‑1 (MQ=255)
             ccGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCa                                                             >  2:1456320/1‑76 (MQ=255)
             ccGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCCTGACCGCGACCGCGACCAAACCa                                                             >  1:6138179/1‑76 (MQ=255)
              cGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGCCCAAACCAt                                                            >  1:5696366/1‑76 (MQ=255)
              cGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCAt                                                            <  2:6056828/76‑1 (MQ=255)
               ggCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATg                                                           >  2:3322110/1‑76 (MQ=255)
               gcccAGCGCGAACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGAGCGTGACCGCGACCGCGACCAAACCATg                                                           <  1:4264992/74‑1 (MQ=255)
                gCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGc                                                          >  1:4573602/1‑76 (MQ=255)
                 ccAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCg                                                         <  2:3744898/76‑1 (MQ=255)
                 ccAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCg                                                         <  1:631758/76‑1 (MQ=255)
                 ccAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCg                                                         <  2:5896518/76‑1 (MQ=255)
                  cAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGt                                                        >  2:4928900/1‑76 (MQ=255)
                  cAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGt                                                        >  2:4724244/1‑76 (MQ=255)
                  cAGAGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGt                                                        >  2:8063608/1‑76 (MQ=255)
                   aGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTa                                                       >  2:1943033/1‑76 (MQ=255)
                   aGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTa                                                       <  2:697191/76‑1 (MQ=255)
                   aGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACAATGCGTa                                                       >  2:4444985/1‑76 (MQ=255)
                    gcgcAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTAc                                                      >  1:9309988/1‑76 (MQ=255)
                    gcgcAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTAc                                                      <  1:1133083/76‑1 (MQ=255)
                     cgcAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTAct                                                     >  1:1296132/1‑76 (MQ=255)
                     cgcAGACCGCGCGCTTTGTGCTAGTCATCGAACGGTGCCAGATCGTGACCGCGCCCGCGACCAAACCATGCGTAct                                                     <  1:3762300/76‑1 (MQ=255)
                     cgcAGACCGCGCGCTTGGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTAct                                                     <  1:6620204/76‑1 (MQ=255)
                      gcAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGAACCAACCATGAGTActc                                                    >  2:9426091/1‑76 (MQ=255)
                       cAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTActct                                                   <  2:10333519/76‑1 (MQ=255)
                       cAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACAAAACCATGCGTActct                                                   >  2:5069853/1‑76 (MQ=255)
                        aGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCtt                                                  >  2:6021387/1‑76 (MQ=255)
                         gACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTc                                                 >  2:4573271/1‑76 (MQ=255)
                         gACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTc                                                 <  1:9734193/76‑1 (MQ=255)
                            cgcgcgCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGaa                                              <  1:2716153/76‑1 (MQ=255)
                              cgcgcTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACa                                            >  2:2834048/1‑76 (MQ=255)
                              cgcgcTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACa                                            <  2:1284088/76‑1 (MQ=255)
                              cgcgcTTGATAGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAAcg                                            <  2:5723511/76‑2 (MQ=255)
                               gcgcTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAg                                           <  1:5065651/76‑1 (MQ=255)
                                cgcTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGc                                          >  2:6878886/1‑76 (MQ=255)
                                  cTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGccc                                        >  2:3645094/1‑76 (MQ=255)
                                  cTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGccc                                        <  2:2952398/76‑1 (MQ=255)
                                   tttGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCt                                       >  1:7453259/1‑76 (MQ=255)
                                    ttGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCtt                                      <  2:2916911/76‑1 (MQ=255)
                                       tgGTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCttttt                                   <  2:5801747/76‑1 (MQ=255)
                                        ggTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTg                                  <  1:10343352/76‑1 (MQ=255)
                                        ggTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTg                                  >  2:9768324/1‑76 (MQ=255)
                                        ggTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTg                                  <  1:35348/76‑1 (MQ=255)
                                        ggTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTg                                  <  1:5409448/76‑1 (MQ=255)
                                        ggTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTg                                  <  1:8379038/76‑1 (MQ=255)
                                         gTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGc                                 <  2:10138254/76‑1 (MQ=255)
                                         gTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGc                                 <  2:524293/76‑1 (MQ=255)
                                         gTAGTCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGc                                 >  1:8720745/1‑76 (MQ=255)
                                         gTAGTCATCGAACGGTACCAGATCGTGCCCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGc                                 <  1:7909554/76‑1 (MQ=255)
                                             tCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGAt                             <  1:2026811/76‑1 (MQ=255)
                                             tCATCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGAt                             <  1:8063608/76‑1 (MQ=255)
                                               aTCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAg                           >  1:8869196/1‑76 (MQ=255)
                                               aTCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAg                           <  1:5901665/76‑1 (MQ=255)
                                               aTCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAg                           <  2:9112568/76‑1 (MQ=255)
                                               aTCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAg                           <  1:605065/76‑1 (MQ=255)
                                               aTCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAg                           >  2:9598542/1‑76 (MQ=255)
                                                tCGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGa                          <  2:7984106/76‑1 (MQ=255)
                                                 cGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGaa                         <  2:40908/76‑1 (MQ=255)
                                                 cGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGaa                         <  2:2133935/76‑1 (MQ=255)
                                                 cGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGaa                         >  2:8128303/1‑76 (MQ=255)
                                                 cGAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGaa                         >  2:3829568/1‑76 (MQ=255)
                                                  gAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAg                        >  1:9652314/1‑76 (MQ=255)
                                                  gAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAg                        <  2:4306715/76‑1 (MQ=255)
                                                  gAACGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAg                        >  2:7683248/1‑76 (MQ=255)
                                                   aaCGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGc                       >  2:6577760/1‑76 (MQ=255)
                                                    aCGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGcc                      <  1:7484758/76‑1 (MQ=255)
                                                    aCGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGcc                      <  2:6717182/76‑1 (MQ=255)
                                                     cGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGccc                     >  1:6971379/1‑76 (MQ=255)
                                                     cGGTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGccc                     <  1:906214/76‑1 (MQ=255)
                                                     cGGGGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGccc                     <  1:6525112/76‑1 (MQ=255)
                                                       gTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAg                   <  1:6431179/76‑1 (MQ=255)
                                                       gTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAg                   >  1:2936547/1‑76 (MQ=255)
                                                       gTGCCAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAg                   >  1:4453579/1‑76 (MQ=255)
                                                           cAGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTcg               >  1:8209295/1‑76 (MQ=255)
                                                            aGATCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGaa                         <  1:4772233/65‑1 (MQ=255)
                                                               tCGTGACCGCGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCggg           >  1:242051/1‑76 (MQ=255)
                                                                 gtgaccgcgaccGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGAt         >  2:7931928/1‑76 (MQ=255)
                                                                   gaccgcgaccgcGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTc       >  2:2587426/1‑76 (MQ=255)
                                                                   gaccgcgaccgcGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTc       >  1:2345336/1‑76 (MQ=255)
                                                                   gaccgcgaccgcGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTc       >  2:4773948/1‑76 (MQ=255)
                                                                   gaccgcgaccgcGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTc       >  2:335435/1‑76 (MQ=255)
                                                                    accgcgaccgcgACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCa      >  1:4102801/1‑76 (MQ=255)
                                                                     ccgcgaccgcgaCCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCAt     <  1:1511690/76‑1 (MQ=255)
                                                                        |                                                                           
TTACCGTTCACCACCGGCCAGCGCAGACCGCGCGCTTTGTGGTAGTCATCGAACGGTGCCAGATCGTGACCGTGACCGCGACCAAACCATGCGTACTCTTCGAACAGCCCTTTTTGCAGATAGAAGCCCAGCTCGCGGGATTCATCGT  >  REL606/2250752‑2250899

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: