Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,677 G→C 17.1% L1045L (CTG→CTC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,6770GC17.1% 174.0 / 26.5 76L1045L (CTG→CTCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (4/9);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

ACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCC  >  minE/834610‑834747
                                                                   |                                                                      
aCAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTa                                                                      <  1:2709250/70‑1 (MQ=255)
 taaCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTa                                                                      <  1:3553339/68‑1 (MQ=255)
  aaCGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt                                                                     >  1:2042326/1‑69 (MQ=255)
   aCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCa                                                                   <  1:965805/70‑1 (MQ=255)
   aCGCGGGATACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAg                                                                  >  1:807224/1‑71 (MQ=255)
    cGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAgg                                                                 <  1:2041313/71‑1 (MQ=255)
       cGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGAt                                                               >  1:1650743/1‑70 (MQ=255)
        ggTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt                                                                     <  1:3104310/63‑1 (MQ=255)
        ggTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTAt                                                                     <  1:2732990/63‑1 (MQ=255)
          ttACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAc                                                            >  1:1316801/1‑70 (MQ=255)
           tACACCAACGTTCATGAGTTAATTCCGTGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCatcc                                                          <  1:3462352/71‑4 (MQ=255)
            acacCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCa                                                          <  1:3838051/70‑1 (MQ=255)
              acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAcc                                                           <  1:7733/67‑1 (MQ=255)
              acCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTg                                                       <  1:3906863/71‑1 (MQ=255)
                cAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                                                     >  1:107192/1‑71 (MQ=255)
                cAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                                                     >  1:1683362/1‑71 (MQ=255)
                    gTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGc                                                  >  1:2720589/1‑70 (MQ=255)
                    gTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCg                                                 >  1:3904524/1‑71 (MQ=255)
                      tCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGt                                                >  1:2369895/1‑70 (MQ=255)
                      tCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGt                                                >  1:1241453/1‑70 (MQ=255)
                       cACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAg                                                         >  1:2554934/1‑60 (MQ=255)
                         cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt                                            <  1:202207/71‑1 (MQ=255)
                         cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt                                            <  1:269066/71‑1 (MQ=255)
                         cGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAtt                                            <  1:2117088/71‑1 (MQ=255)
                          gAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAttt                                           <  1:205397/71‑1 (MQ=255)
                             cTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAttt                                           <  1:906536/68‑1 (MQ=255)
                               gATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCg                                         <  1:519447/68‑1 (MQ=255)
                                aTCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTg                                                       <  1:2906606/53‑1 (MQ=255)
                                aTCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGt                                       <  1:1306503/69‑1 (MQ=255)
                                aTCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGa                                     <  1:3213661/71‑1 (MQ=255)
                                 tCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                                                     >  1:1478975/1‑54 (MQ=255)
                                 tCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGa                                                     >  1:1831850/1‑54 (MQ=255)
                                 tCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGa                                     >  1:3933976/1‑70 (MQ=255)
                                   ccATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAg                                  >  1:3103255/1‑71 (MQ=255)
                                    cATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGc                                 >  1:3604711/1‑71 (MQ=255)
                                     aTGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGcc                                >  1:1294340/1‑71 (MQ=255)
                                      tGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAcc                                                           <  1:631000/43‑1 (MQ=255)
                                      tGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCAcc                                                           <  1:3759842/43‑1 (MQ=255)
                                      tGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCct                               >  1:649083/1‑71 (MQ=255)
                                      tGGCGCACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCct                               <  1:3139982/71‑1 (MQ=255)
                                        gcgcACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGa                                     <  1:598416/63‑1 (MQ=255)
                                        gcgcACGCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGa                                     <  1:3624194/63‑1 (MQ=255)
                                        gCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCctgc                             <  1:1219791/71‑1 (MQ=255)
                                         cGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCctg                              <  1:2218364/69‑1 (MQ=255)
                                          gTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCctgctg                           <  1:2237273/71‑1 (MQ=255)
                                            aCGCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCtggtggc                         >  1:2513776/1‑70 (MQ=255)
                                             cGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGt                                       >  1:1391373/1‑56 (MQ=255)
                                             cGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGtt                        <  1:1766230/71‑1 (MQ=255)
                                             cGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGt                         <  1:1069052/70‑1 (MQ=255)
                                              gCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCtggtggct                        >  1:921693/1‑68 (MQ=255)
                                              gCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCtggtggct                        >  1:1339550/1‑68 (MQ=255)
                                              gCTGTCGGGACGCCAGCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCCTGGTGGctt                       <  1:1594662/71‑1 (MQ=255)
                                              gCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGtt                        >  1:3119623/1‑70 (MQ=255)
                                               ctctctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAttt                                           >  1:1865375/1‑50 (MQ=255)
                                               ctctctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAttt                                           >  1:1459784/1‑50 (MQ=255)
                                               ctctctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGAttt                                           >  1:602307/1‑50 (MQ=255)
                                               ctctctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTa                      <  1:211357/71‑1 (MQ=255)
                                                   ctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGaaa                                   >  1:2364691/1‑54 (MQ=255)
                                                   ctGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGt                  >  1:372095/1‑71 (MQ=255)
                                                    tGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTc                 >  1:1663510/1‑71 (MQ=255)
                                                       tcgtcAGCAACTGTATCAGGATCACCAGTGGATGc                                                  <  1:2082184/35‑1 (MQ=255)
                                                         gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg            <  1:3034171/71‑1 (MQ=255)
                                                         gtcAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgccg            <  1:432688/71‑1 (MQ=255)
                                                            agcagcTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCct                               <  1:463309/49‑1 (MQ=255)
                                                            aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc             >  1:1179273/1‑67 (MQ=255)
                                                            aGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTccgcc             >  1:3103029/1‑67 (MQ=255)
                                                             gCAGCTCTATCAGGATCATCCGTGGATGCGTGCTTTTGGTGAAAGCct                               <  1:2479820/48‑1 (MQ=255)
                                                             gCAACTGTATCAGGATCACCAGTGGATTCGTGATTTCGGTGAAAGCCTGCTGGTTTATCg                   <  1:2242873/60‑1 (MQ=255)
                                                             gCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAg                                  >  1:3530863/1‑45 (MQ=255)
                                                             gCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCg        <  1:1606177/71‑1 (MQ=255)
                                                              cAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCGTGCTGGTTTa                      >  1:3890062/1‑56 (MQ=255)
                                                                aCTGTATCAGGATCACCAGTGGATGGGTGATTTCGGTGAAAGc                                 >  1:2328104/1‑43 (MQ=255)
                                                                aCTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGaca     >  1:817567/1‑71 (MQ=255)
                                                                aCTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGa       <  1:3585401/69‑1 (MQ=255)
                                                                 cTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGacac    >  1:838356/1‑71 (MQ=255)
                                                                  tGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATc         <  1:3039403/65‑1 (MQ=255)
                                                                   gTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACAccc  >  1:719501/1‑71 (MQ=255)
                                                                   |                                                                      
ACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCGTACGCTCTCTGGTCGTCAGCAACTGTATCAGGATCACCAGTGGATGCGTGATTTCGGTGAAAGCCTGCTGGTTTATCGTCCGCCGATCGACACCC  >  minE/834610‑834747

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