Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,541,077 G→C 14.0% A505G (GCC→GGC)  narZ ← nitrate reductase Z subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,541,0770GC14.0% 103.9 / 12.5 48A505G (GCC→GGC) narZnitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base C (5/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

GTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTT  >  NC_000913/1540938‑1541218
                                                                                                                                           |                                                                                                                                             
gTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGCGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGCAGGGTTCAAGATGCTGGCGCGGCACGCCGGGAGATTGGCCGCCCCAgggcg                                                                                                                                            >  2:460127/1‑140 (MQ=25)
  ccATGTGACAGCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTggg                                                                                                                                          <  2:452064/143‑1 (MQ=255)
   cATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgggg                                                                                                                                         >  2:397505/1‑143 (MQ=255)
   cATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTgggg                                                                                                                                         >  1:435139/1‑143 (MQ=255)
     tgtgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGgtg                                                                                                                                       >  1:2025/1‑143 (MQ=255)
     tgtgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCCGGGgtg                                                                                                                                       >  2:574413/1‑143 (MQ=255)
          tACCAGTGGTTAACACCGGCGCCGAGGATAATCATCGAGCGCCCATGCGTCTTATGGGCGGTATCGGCAAAATCACGGGCGATGGTGTCGATATACTGGCGGGGCACGCCGGCCACTCGCCCACCCCGGGCTGGGGGGTACgg                                                                                                                                  >  2:646972/1‑143 (MQ=255)
            ccAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGtt                                                                                                                                <  1:397505/143‑1 (MQ=255)
               gccgTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGa                                                                                                                             <  1:471716/140‑1 (MQ=255)
               gTGGTTAACACCTGCGCCGAGGATAATCATCGTGCGCCTATTCGGTGTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGa                                                                                                                             <  1:533573/143‑1 (MQ=255)
                  ggtAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAttt                                                                                                                          <  1:260126/141‑1 (MQ=255)
                           tGCGCCGAGGATAATCATCGAGCGCCCATGCGCTTTATGGGCGGTATCGGCACATTCCCGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTCCGGTTTGATTTCAGCGTAAt                                                                                                                 >  1:240529/1‑143 (MQ=255)
                            gcgcCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                       >  1:303792/1‑120 (MQ=255)
                            gcgcCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                       <  2:299168/120‑1 (MQ=255)
                            gcgcCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                       <  2:303792/120‑1 (MQ=255)
                            gcgcCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtg                                                                                                                                       >  1:299168/1‑120 (MQ=255)
                                         tcatcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCACATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTCATCTTTCGCACTGttt                                                                                                   >  2:535663/1‑143 (MQ=255)
                                            tcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCa                                                                                                >  1:659660/1‑143 (MQ=255)
                                                 cacgcATGCGGTTTACGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc                                                                                           <  1:92773/139‑1 (MQ=255)
                                                  gcCCATGCGTTTTATGGGCGGTATCGTCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCCCCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTcc                                                                                          >  1:431899/1‑143 (MQ=255)
                                                  gcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATCTACTGGCGCGGCACGCCGGTAATTTGCTCGCCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTcc                                                                                          <  1:574413/143‑1 (MQ=255)
                                                      aTGCGTCTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGcc                                                                                      <  1:438435/143‑1 (MQ=255)
                                                          gTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGGTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTAGTCTTCCAGCCcgcg                                                                                  >  1:165162/1‑143 (MQ=255)
                                                            tttGTGGGCGGTATCGGCAAGTTCGCGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGAt                                                                                <  2:175176/143‑1 (MQ=255)
                                                              tctGGGCGGTACCGGCAAATTCACGGGCGGTGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCCCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCg                                                                              <  2:425861/141‑1 (MQ=255)
                                                                   gcgtgaTCGGCAAATTCACGGGAGCTGGTTTCAATCTACTGGCGCGGCACGGCGATCATTTGCTTAGCCCTGGCTGGTGTGTACGGTTTGATTTCAGCGTCATCTGTGCCACTGTTTTCATCTTCCAGCCCGCGATCGAGGcc                                                                         <  2:385780/138‑1 (MQ=255)
                                                                    cGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTATTCTTTCGCACGGGTCTCATCTTCCAGCCCGCGATCGAGGCCg                                                                        >  1:741491/1‑143 (MQ=255)
                                                                      gtctCGGCAAATGCGCGGGCGATGGTTTCAATCTACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTCTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTa                                                                      <  2:326067/140‑1 (MQ=255)
                                                                         tCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAAtt                                                                   >  1:630254/1‑143 (MQ=255)
                                                                           ggCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCg                                                                 >  1:484057/1‑143 (MQ=255)
                                                                               aaTTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAg                                                             <  1:241190/143‑1 (MQ=255)
                                                                                 ttCACGGGCGATGGTTTCAATATACTGGCGCGGCCCGCCGGTAATCCGCTCACCCCAGGCTGCGGTGCACGGTTTGATTCCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTa                                                           <  2:599422/143‑1 (MQ=255)
                                                                                          gttgctttCAATATACTCGCGCGGCACGCCGGTAATTTGCTCCCCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAACCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCAt                                                  <  1:466264/138‑1 (MQ=255)
                                                                                             ggTTTCAATAGACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATaaa                                               <  2:657002/143‑1 (MQ=255)
                                                                                              gTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAAc                                              >  1:677083/1‑143 (MQ=255)
                                                                                                ttCGATATACTTGCGCGGCTCGTCGGTACTTGGCTCACCTCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAAGCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGc                                            <  2:526481/143‑1 (MQ=255)
                                                                                                   attgtgcTGGCGCGGCACGCCGGTAATTTGCTTACCCCAGGCTGGGGTGTACGGTTTGATTGCAGCGTAATCTTTCGCACTGTTGTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGa                                         <  2:233657/137‑1 (MQ=255)
                                                                                                    atatACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGac                                        >  1:715655/1‑143 (MQ=255)
                                                                                                               cgGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGCACGAGACGTTATTTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAgg                             >  1:330444/1‑143 (MQ=255)
                                                                                                               cgGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACAGATTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGCCAgg                             >  1:331427/1‑143 (MQ=255)
                                                                                                                   aCGCCGGTAATTTGCTCCCCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGtt                         <  2:431899/143‑1 (MQ=255)
                                                                                                                      ccGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGCTTTCATCTTCCAGCCCGCGATCGAGGCCGTCATTCGCCAGTACCAAGTCATCAACGCCGACCaccggaccggggcggcc                      >  2:296913/1‑136 (MQ=255)
                                                                                                                            aaccgtacACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGAc                >  1:404319/8‑143 (MQ=255)
                                                                                                                                gCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAAc            <  2:240529/143‑1 (MQ=255)
                                                                                                                                  tCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTGGCCATCGAACAACGt          >  2:256134/1‑143 (MQ=255)
                                                                                                                                   cACCCCAGCCTGGGGTGTACGGTTTGATTTCAGCGTGATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACCGGTGTTGCCATCGACTAACGTc         >  1:206467/1‑143 (MQ=255)
                                                                                                                                     ccccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCCCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAg       <  2:609160/143‑1 (MQ=255)
                                                                                                                                      cccAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCAACGGACAGGTGTTGCCATCGACTACCGTCAgg      >  1:712236/1‑143 (MQ=255)
                                                                                                                                      cccAGGCTGGGGTGTACGGTTTGATTTCAGCGGAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGGCAgg      >  1:385777/1‑143 (MQ=255)
                                                                                                                                          ggCTGGGGTGTACGGTTTGATGTCAGCGGAATCGTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTCATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGtttt  <  2:177216/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTT  >  NC_000913/1540938‑1541218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: