Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,452 G→C 6.2% L390L (CTG→CTC narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,4520GC 6.2% 195.0 / 3.0 64L390L (CTG→CTCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base C (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGC  >  minE/836383‑836519
                                                                     |                                                                   
cgcgGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTgg                                                                    >  1:2122209/1‑71 (MQ=255)
 gcgGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGc                                                                   >  1:3698168/1‑71 (MQ=255)
   gggTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCg                                                                  >  1:3883790/1‑70 (MQ=255)
   gggTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGa                                                                 <  1:2927223/71‑1 (MQ=255)
             ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                         <  1:1047054/69‑1 (MQ=255)
             ggTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                         <  1:1042115/69‑1 (MQ=255)
               tAGCAACGGCATTCTGCCAGACGTAGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATc                                                              >  1:2414971/1‑62 (MQ=255)
                 gCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctg                                                            <  1:3048575/62‑1 (MQ=255)
                   aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                               >  1:3515285/1‑57 (MQ=255)
                   aaCGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAt                                                               >  1:3933056/1‑57 (MQ=255)
                     cgtcATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTg                                               <  1:3695159/68‑1 (MQ=255)
                      ggCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGa                                              <  1:3049400/71‑1 (MQ=255)
                      ggCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGa                                              <  1:2825019/71‑1 (MQ=255)
                        cATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATa                                            >  1:324108/1‑71 (MQ=255)
                            cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTg                                               >  1:2542316/1‑64 (MQ=255)
                            cTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCa                                         >  1:1357481/1‑70 (MQ=255)
                             tGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                           >  1:1358649/1‑67 (MQ=255)
                             tGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCaaa                                       >  1:1501839/1‑71 (MQ=255)
                               ccAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                      <  1:763241/70‑1 (MQ=255)
                                 aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                                >  1:2182173/1‑58 (MQ=255)
                                 aGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                                   <  1:3690589/71‑1 (MQ=255)
                                   aCGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAAtcttct                                                          <  1:1505454/46‑1 (MQ=255)
                                      tCGAAAGCCTGTGTATTCCGGTGCAATATCTCGCCaa                                                                <  1:759919/37‑1 (MQ=16)
                                      tCGAAAGCCTGCGTATTCCGGTGCAATATCTCGCCaa                                                                <  1:468029/37‑1 (MQ=21)
                                      tCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGt                                                <  1:900946/53‑1 (MQ=255)
                                        gAAAGCCTGCGTATTCCGGTGCAATATCTCGCAAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTgcg                            >  1:305060/1‑71 (MQ=255)
                                        gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGaa                                                                >  1:1727636/1‑35 (MQ=255)
                                        gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctgctg                                                         >  1:2621622/1‑42 (MQ=255)
                                        gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATctg                                                            <  1:573089/39‑1 (MQ=255)
                                        gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                                   <  1:594163/48‑1 (MQ=255)
                                        gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                                   <  1:599095/48‑1 (MQ=255)
                                        gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                           >  1:2178356/1‑56 (MQ=255)
                                        gAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcg                            >  1:1760449/1‑71 (MQ=255)
                                          aaGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGAccgcc                                                   >  1:3962629/1‑46 (MQ=255)
                                          aaGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                          >  1:70665/1‑71 (MQ=255)
                                           aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcgc                         <  1:2966477/71‑1 (MQ=255)
                                           aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgcg                          >  1:3622679/1‑70 (MQ=255)
                                           aGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGATGATACCAAACCGGTACTgcgcg                          >  1:3739535/1‑70 (MQ=255)
                                             ccTGCGTATTCCGGTACAGTATCTGGCGAATCTGGTGACCGCCGGTGATACCAAACCGGTACTGCGCGCAc                       <  1:60223/71‑1 (MQ=255)
                                              cTGCGTATTACGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACt                      <  1:2159068/71‑1 (MQ=255)
                                                gCGTATTCCGGTGCAATATCTCGCCAATATGTTGAGTGCCGGCGATACCGGTCCGGTACTGCGGGCGCTg                     >  1:1951945/1‑70 (MQ=255)
                                                 cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCaaa                                       >  1:3933348/1‑51 (MQ=255)
                                                 cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                                   <  1:3096074/55‑1 (MQ=255)
                                                 cGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCgg                                   <  1:1218914/55‑1 (MQ=255)
                                                     ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTgcgc                           >  1:840668/1‑59 (MQ=255)
                                                     ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg                >  1:578007/1‑70 (MQ=255)
                                                     ttCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACg                >  1:1156982/1‑70 (MQ=255)
                                                       ccGGTACAGTATCTGGCGCATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTa              >  1:1809069/1‑70 (MQ=255)
                                                       ccGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATAc                                           >  1:175459/1‑41 (MQ=255)
                                                         ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc           <  1:172290/71‑1 (MQ=255)
                                                         ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc           <  1:36380/71‑1 (MQ=255)
                                                         ggTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGc           <  1:412767/71‑1 (MQ=255)
                                                           tACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTg         <  1:632821/71‑1 (MQ=255)
                                                             cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc       <  1:2347025/71‑1 (MQ=255)
                                                             cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc       <  1:158314/71‑1 (MQ=255)
                                                             cAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc       <  1:153496/71‑1 (MQ=255)
                                                              aGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc       <  1:3795209/70‑1 (MQ=255)
                                                               gTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTAt             <  1:3923861/63‑1 (MQ=255)
                                                                 aTCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGa     <  1:3393651/69‑1 (MQ=255)
                                                                  tCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                      <  1:823780/35‑1 (MQ=255)
                                                                  tCTGGCGAATCTGCTGACCGCCGGTGATACCAAAc                                      <  1:2338074/35‑1 (MQ=255)
                                                                  tCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGc  >  1:1989677/1‑71 (MQ=255)
                                                                  tCTGGCGAATCTGCTGACCGCCGGTGAAACCAAACCGGTACTGCGCGCACTGAAACGTa              >  1:2015765/1‑59 (MQ=255)
                                                                   cTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCt          >  1:4065899/1‑62 (MQ=255)
                                                                   cTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCt          >  1:1698519/1‑62 (MQ=255)
                                                                     ggCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGc       >  1:306660/1‑63 (MQ=255)
                                                                     |                                                                   
CGCGGGTGAGCTGGGTAGCAACGGCATTCTGCCAGACGTCGAAAGCCTGCGTATTCCGGTACAGTATCTGGCGAATCTGCTGACCGCCGGTGATACCAAACCGGTACTGCGCGCACTGAAACGTATGCTGGCGATGC  >  minE/836383‑836519

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: