Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,535,870 T→A 8.9% intergenic (‑12/+67) yddE ← / ← narV PhzC‑PhzF family protein/nitrate reductase 2 (NRZ), gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,535,8700TA 8.9% 176.8 / 8.3 56intergenic (‑12/+67)yddE/narVPhzC‑PhzF family protein/nitrate reductase 2 (NRZ), gamma subunit
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base A (2/3);  total (26/30)
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.63e-01

CCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGT  >  NC_000913/1535723‑1536014
                                                                                                                                                   |                                                                                                                                                
ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCtg                                                                                                                                                 >  4:264716/1‑149 (MQ=255)
   aaTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGa                                                                                                                                              <  2:68446/149‑1 (MQ=255)
    aTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGAc                                                                                                                                             <  4:276433/149‑1 (MQ=255)
    aTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGAc                                                                                                                                             <  2:12866/149‑1 (MQ=255)
        ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGga                                                                                                                                         >  4:378506/1‑149 (MQ=255)
              tCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCAGTCACAGGAGAAcgt                                                                                                                                   <  4:528317/149‑3 (MQ=255)
                  gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCAGTCACAGGAGAAcgttaat                                                                                                                               <  3:592674/149‑7 (MQ=255)
                       aTAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCt                                                                                                                          <  3:474916/149‑1 (MQ=255)
                           gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                                                                                                      <  1:127729/149‑1 (MQ=255)
                           gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                                                                                                      >  2:320047/1‑149 (MQ=255)
                           gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                                                                                                      >  4:310930/1‑149 (MQ=255)
                           gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                                                                                                      <  4:576117/149‑1 (MQ=255)
                            cTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGACGCGAAACGGTTGTGAGGTAAAAGCATCGACGTTGTACACCTGCGGGTTCATTAACGTTCTCCTGTGAATGGAGAACTATTATAGCCTGGAAg                                                                                                                     >  2:19122/1‑149 (MQ=255)
                              gCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCCGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTg                                                                                                                   >  2:153913/1‑149 (MQ=255)
                                 tttGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGcc                                                                                                                >  2:516838/1‑149 (MQ=255)
                                      gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGgaga                                                                                                           <  1:19122/149‑1 (MQ=255)
                                      gcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTAGAGAACTATCATAGCCTGCAAGTGGCCGgaga                                                                                                           >  2:408930/1‑149 (MQ=255)
                                           gCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                                                                                                      <  3:378506/149‑1 (MQ=255)
                                             tAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCAATAACGGTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAgg                                                                                                    >  2:609836/1‑149 (MQ=255)
                                                gATTATCGGCGGGGAAAACAACACTGGCAGAATTTCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCt                                                                                                 >  2:374700/1‑149 (MQ=255)
                                                gATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCt                                                                                                 <  1:320047/149‑1 (MQ=255)
                                                 aTTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTa                                                                                                >  2:46935/1‑149 (MQ=255)
                                                           gggAAAACAAAACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg                                                                                      >  3:38464/1‑149 (MQ=255)
                                                                  caacaCCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAt                                                                               >  2:603130/1‑149 (MQ=255)
                                                                     cacCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTAt                                                                            >  2:294572/1‑149 (MQ=255)
                                                                      acCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATc                                                                           <  1:153913/149‑1 (MQ=255)
                                                                       ccGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACCTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcg                                                                          <  4:543544/149‑1 (MQ=255)
                                                                          gCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATcgccg                                                                       >  4:501342/1‑149 (MQ=255)
                                                                                ttGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                                                                                                      >  4:371232/1‑112 (MQ=255)
                                                                                ttGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                                                                                                      <  3:371232/112‑1 (MQ=255)
                                                                                          aCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                       >  2:341299/1‑149 (MQ=255)
                                                                                          aCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                       <  3:501342/149‑1 (MQ=255)
                                                                                          aCGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGa                                                       <  3:265532/149‑1 (MQ=255)
                                                                                                  gAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATcgacg                                               >  3:216353/1‑149 (MQ=255)
                                                                                                   aGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCGTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATcgacga                                              <  1:141490/149‑1 (MQ=255)
                                                                                                   aGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTAAAATTATCGCCGCGAACGCACAATTTGATATcgacga                                              <  1:374700/149‑1 (MQ=255)
                                                                                                        aaaGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTaaa                                         <  4:216353/149‑1 (MQ=255)
                                                                                                           gCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTa                                      <  3:483822/149‑1 (MQ=255)
                                                                                                                 aCGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCaaa                                <  3:544763/149‑1 (MQ=255)
                                                                                                                     ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg                                                                                      >  3:198961/1‑91 (MQ=255)
                                                                                                                     ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg                                                                                      <  4:198961/91‑1 (MQ=255)
                                                                                                                     ggTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGgc                            <  4:321405/149‑1 (MQ=255)
                                                                                                                      gTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaa                                                                                 >  4:238709/1‑95 (MQ=255)
                                                                                                                      gTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaa                                                                                 <  3:238709/95‑1 (MQ=255)
                                                                                                                                 gTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCcacaca                >  1:398421/1‑149 (MQ=255)
                                                                                                                                      aTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAc           >  2:73684/1‑149 (MQ=255)
                                                                                                                                      aTTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAc           >  2:73736/1‑149 (MQ=255)
                                                                                                                                       tgatagTTCTCCAGTCACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCAcc          >  1:255465/6‑149 (MQ=255)
                                                                                                                                           agTTCTCCAGTCACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAAc                               <  2:274282/123‑1 (MQ=255)
                                                                                                                                           agTTCTCCAGTCACAGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAAc                               >  1:274282/2‑124 (MQ=255)
                                                                                                                                             ttCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGata                                                     >  4:211474/1‑100 (MQ=255)
                                                                                                                                             ttCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGata                                                     <  3:211474/100‑1 (MQ=255)
                                                                                                                                             ttCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCgg    <  4:480243/149‑1 (MQ=255)
                                                                                                                                              tctcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCggg   >  4:247352/1‑149 (MQ=255)
                                                                                                                                              tctcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCggg   <  1:73736/149‑1 (MQ=255)
                                                                                                                                               ctcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGTGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGt  <  4:426102/149‑1 (MQ=255)
                                                                                                                                               ctcCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGTGCGCTCCACACATGCACCAGTCGGGt  <  1:294572/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAATTATCGCCGCGAACGCACAATTTGATATCGACGAGTAAAGTACTCAAACGGCGCGCTCCACACATGCACCAGTCGGGT  >  NC_000913/1535723‑1536014

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑