Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,540,993 T→C 9.1% H533R (CAT→CGT)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,540,9930TC 9.1% 124.3 / 3.0 44H533R (CAT→CGT) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (22/18);  new base C (2/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGAT  >  NC_000913/1540856‑1541140
                                                                                                                                         |                                                                                                                                                   
cGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATgg                                                                                                                                          >  7:43131/1‑149 (MQ=255)
    aTAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCgg                                                                                                                                      >  8:223351/1‑149 (MQ=255)
        tgtgCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGtt                                                                                                                                                <  1:39093/135‑1 (MQ=255)
        tgtgCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGtt                                                                                                                                                >  2:39093/1‑135 (MQ=255)
                 ccgccACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAAtt                                                                                                                          >  8:63512/1‑148 (MQ=255)
                   gccACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCAc                                                                                                                       >  7:113831/1‑149 (MQ=255)
                   gccACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCAGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAAGCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCAc                                                                                                                       <  8:179076/149‑1 (MQ=255)
                     cACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGAGGTAGTTCATGTCCATGTGATAACAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACgg                                                                                                                     <  3:23458/149‑1 (MQ=255)
                           tttGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCACGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATg                                                                                                               >  5:83404/1‑149 (MQ=255)
                           tttGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCACGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATg                                                                                                               >  5:23730/1‑149 (MQ=255)
                                 cGACACAACCGCAGAAGATCAGCATATTGATCCTCCCGCGGTAGTTCATGTCCATGTGATACCAGAGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTAGCGGCAAATTCACGGGAGATGGTTTc                                                                                                          >  5:19962/1‑148 (MQ=255)
                                       aaCCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATAc                                                                                                   >  3:215486/1‑149 (MQ=255)
                                         ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGATCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTg                                                                                                 <  6:126840/149‑1 (MQ=255)
                                         ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTg                                                                                                 >  3:53083/1‑149 (MQ=255)
                                         ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTg                                                                                                 >  6:40631/1‑149 (MQ=255)
                                            cAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGgcg                                                                                              >  2:176655/1‑149 (MQ=255)
                                                       cATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCgg                                                                                   <  8:113831/149‑1 (MQ=255)
                                                           ttGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATCGTTTCAATATACTGGCGCGGCACGCCGGTAAt                                                                               >  7:144499/1‑149 (MQ=255)
                                                                    ccGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCAcc                                                                      >  7:220033/1‑149 (MQ=255)
                                                                        ggTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAg                                                                  >  5:161694/1‑149 (MQ=255)
                                                                                   tCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtga                                                        <  4:215486/148‑2 (MQ=255)
                                                                                   tCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGCCTGGGgtga                                                        <  8:144499/148‑2 (MQ=255)
                                                                                    ccATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTAc                                                      >  8:227690/1‑149 (MQ=255)
                                                                                     cATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACg                                                     >  4:42280/1‑149 (MQ=255)
                                                                                            tACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCACGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAt                                              <  6:83404/149‑1 (MQ=255)
                                                                                            tACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCACGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAt                                              <  6:23730/149‑1 (MQ=255)
                                                                                             aCCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATggg                                                                                                                                         >  2:170216/1‑57 (MQ=255)
                                                                                             aCCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATggg                                                                                                                                         <  1:170216/57‑1 (MQ=255)
                                                                                             aCCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAtt                                             >  1:149497/1‑149 (MQ=255)
                                                                                               cAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGGTATATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGAttt                                            >  2:111497/1‑148 (MQ=255)
                                                                                                aggtgTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGTCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCa                                          <  8:116465/145‑1 (MQ=255)
                                                                                                aGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCa                                          <  6:161694/149‑1 (MQ=255)
                                                                                                 gTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                         >  6:184315/1‑149 (MQ=255)
                                                                                                 gTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAg                                         >  8:221576/1‑149 (MQ=255)
                                                                                                       aaCACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAAt                                   >  1:106968/1‑149 (MQ=255)
                                                                                                            cTGCGCCGAGGATAATCCTCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTc                              <  1:146561/149‑1 (MQ=255)
                                                                                                                     ggATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGttt                     <  5:40631/149‑1 (MQ=255)
                                                                                                                       aTAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTc                   <  7:227690/149‑1 (MQ=255)
                                                                                                                             atcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCtt              <  2:149497/148‑1 (MQ=255)
                                                                                                                             atcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc             <  7:63512/149‑1 (MQ=255)
                                                                                                                             atcGAGCGCCCATGCGGTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGGTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc             <  6:19962/149‑1 (MQ=255)
                                                                                                                                gAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCa           <  4:53083/148‑1 (MQ=255)
                                                                                                                                     cccATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCcgc     <  7:223351/149‑1 (MQ=255)
                                                                                                                                        aTGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGAt  >  2:3650/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
CGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGAT  >  NC_000913/1540856‑1541140

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