Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,851 C→T 25.0% S1103S (AGC→AGT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,8510CT25.0% 144.6 / 46.3 84S1103S (AGC→AGTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base T (8/13);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

CGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAG  >  minE/834782‑834921
                                                                     |                                                                      
cGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGc                                                                        >  1:4637303/1‑70 (MQ=255)
 ggCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGa                                                                      >  1:2051331/1‑71 (MQ=255)
   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:3647857/71‑1 (MQ=255)
                gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                                       <  1:4230928/71‑1 (MQ=255)
                gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatgct                                                       <  1:778149/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:2037702/3‑70 (MQ=255)
                    ttaaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcg                                                   >  1:2673952/3‑70 (MQ=255)
                      aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGCCtt                                                  >  1:456030/1‑70 (MQ=255)
                      aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                                 <  1:5208427/71‑1 (MQ=255)
                      aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                                 <  1:3100376/71‑1 (MQ=255)
                      aaCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACttt                                                 <  1:2071783/71‑1 (MQ=255)
                      aaCTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCTGCTAATGCTGAcgct                                                 >  1:270191/1‑68 (MQ=255)
                       aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt                                                  >  1:1695595/1‑69 (MQ=255)
                       aCTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACtt                                                  >  1:1076048/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)
                          tCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGaaa                                                                    <  1:4457032/48‑3 (MQ=19)
                            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:1606080/48‑1 (MQ=19)
                                 gccgcACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACctgct                                                             <  1:1119767/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:434018/45‑1 (MQ=255)
                                        cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                                         <  1:3101428/45‑1 (MQ=255)
                                        cAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctgatg                                                         <  1:5671764/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:4990428/1‑39 (MQ=255)
                                                       cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                                >  1:3318768/1‑39 (MQ=255)
                                                       cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                                >  1:5580994/1‑39 (MQ=255)
                                                       cactccacCTACAGCGACAACCTGCTGATGCTGACTTTa                                                >  1:4188080/1‑39 (MQ=255)
                                                       cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTcg                                           >  1:5123155/1‑44 (MQ=255)
                                                       cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg                 >  1:2545110/1‑70 (MQ=255)
                                                       cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg                 >  1:4805053/1‑70 (MQ=255)
                                                       cactccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTg                 <  1:4499858/70‑1 (MQ=255)
                                                          tccacCTACAGCGACAACCTGCTGATGCTGACTTTAGGt                                             >  1:1228432/1‑39 (MQ=255)
                                                            cacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGcc            <  1:5009583/70‑1 (MQ=255)
                                                            cacCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGcc            <  1:4441411/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)
                                                                     |                                                                      
CGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGCTGACTTTAGGTCGCGGTGGTCCGGTGGTCTGGTTGAGTGAAGCCGATGCCAAAG  >  minE/834782‑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: