Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,281,264 C→T 7.5% S467S (AGC→AGT)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,281,2640CT 7.5% 164.9 / 2.9 53S467S (AGC→AGT) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base T (2/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GGTAGCCAGGATGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGATTACCGGCGTTT  >  NC_000913/1281127‑1281406
                                                                                                                                         |                                                                                                                                              
ggTAGCCAGGATGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAACCGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTgg                                                                                                                                      <  6:314556/148‑1 (MQ=255)
ggTAGCCAGGATGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTg                                                                                                                                       <  2:127349/147‑1 (MQ=255)
  tAGCCAGGATGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTg                                                                                                                                    <  4:167466/148‑1 (MQ=255)
       aGGATGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAAAGCCTGCAACTGGCTGATGGCAGCACAGCCCTGGTGaccac                                                                                                                               >  2:1692/1‑148 (MQ=255)
        ggATGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACT‑GAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTg                                                                                                                                       >  8:49155/1‑138 (MQ=255)
        ggATGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACT‑GAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTg                                                                                                                                       <  7:49155/138‑1 (MQ=255)
          aTGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGt                                                                                                                            <  2:236688/148‑1 (MQ=255)
             agaTCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTTCCGGTGAAACGCCTGCAACTGGCTCATGGCAGCACCGCCCTGGTGACCACCGTTTa                                                                                                                         >  2:66675/1‑148 (MQ=255)
                  gCTGAGGTATGCTTCCCGTACTTTGGTTGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGTCTGATGGCAGCACCGCCCTGTTGACCACCGTTTATGAt                                                                                                                     <  7:328936/147‑1 (MQ=255)
                      aGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTg                                                                                                                  <  3:195012/146‑1 (MQ=255)
                            gCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTgg                                                                                                           >  5:38609/1‑147 (MQ=255)
                                          ggtggCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTg                                                                                             <  8:62005/147‑1 (MQ=255)
                                              gCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAAc                                                                                         <  3:262977/147‑1 (MQ=255)
                                                     cACGGAACACTTCAACTAAGTGGAACTGGAAAACGTGCTGGTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTg                                                                                 <  2:209578/148‑1 (MQ=255)
                                                         gAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGACCTGAac                                                                              >  3:149817/1‑147 (MQ=255)
                                                          aaCACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAacga                                                                            >  5:377942/1‑148 (MQ=255)
                                                            cacTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGACCTGAacgacg                                                                          <  4:149817/148‑1 (MQ=255)
                                                                tCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCt                                                                                                                                        <  4:231939/82‑1 (MQ=255)
                                                                tCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCt                                                                                                                                        >  3:231939/1‑82 (MQ=255)
                                                                      aaGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAacga                                                                            <  3:167576/136‑1 (MQ=255)
                                                                      aaGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAacga                                                                            >  4:167576/1‑136 (MQ=255)
                                                                       aGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTAc                                                                                                   >  2:37494/1‑112 (MQ=255)
                                                                       aGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTAc                                                                                                   <  1:37494/112‑1 (MQ=255)
                                                                              cTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACt                                                                                                     >  1:160769/1‑103 (MQ=255)
                                                                              cTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACt                                                                                                     <  2:160769/103‑1 (MQ=255)
                                                                                ggAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGTACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTAt                                                      >  2:244126/1‑148 (MQ=255)
                                                                                ggAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGTACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTAt                                                      >  6:24981/1‑148 (MQ=255)
                                                                                  aaaaCGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGa                                                    >  8:92927/1‑148 (MQ=255)
                                                                                     aCGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGa                                                 >  7:327311/1‑148 (MQ=255)
                                                                                         gctgctgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGAtgt                                              >  8:101877/1‑147 (MQ=255)
                                                                                            gctgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTg                                                                                             <  8:5502/97‑1 (MQ=255)
                                                                                            gctgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTg                                                                                             >  7:5502/1‑97 (MQ=255)
                                                                                             ctgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGaa                                           <  8:327311/146‑1 (MQ=255)
                                                                                                cacaAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGc                                        <  6:38609/146‑1 (MQ=255)
                                                                                                cacaAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCtt                                      <  3:216476/148‑1 (MQ=255)
                                                                                                cacaAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCtt                                      <  6:377942/148‑1 (MQ=255)
                                                                                                 acaAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCt                                                        <  7:302904/129‑1 (MQ=255)
                                                                                                 acaAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCt                                                        >  8:302904/1‑129 (MQ=255)
                                                                                                  caAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTtat                                    <  7:101877/148‑1 (MQ=255)
                                                                                                    aaCTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATAc                                  >  1:287515/1‑148 (MQ=255)
                                                                                                           ggTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGcc                           <  3:119547/148‑1 (MQ=255)
                                                                                                            gTGAAACGCCTGCAACTGGCTGATGGCAGTACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCt                          <  1:244126/148‑1 (MQ=255)
                                                                                                            gTGAAACGCCTGCAACTGGCTGATGGCAGTACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCt                          <  5:24981/148‑1 (MQ=255)
                                                                                                               aaaCGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGTAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGcc                           <  3:29994/144‑1 (MQ=255)
                                                                                                               aaaCGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACAACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGcc                           >  4:29994/1‑144 (MQ=255)
                                                                                                                aaCGCCTGCAACGTGCTGATGTCAGCACCGCCCTGGGGACCCCCGTTTATTATCTGACGCTGGCAAACTACGGTCTGGAACGGGGCCTGACCGACGTTAACTGTGCCACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTggg                       <  5:15118/147‑1 (MQ=255)
                                                                                                                       gCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGaa                                                                                          >  7:62549/1‑73 (MQ=255)
                                                                                                                       gCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGaa                                                                                          <  8:62549/73‑1 (MQ=255)
                                                                                                                       gCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAg                  >  1:129006/1‑145 (MQ=255)
                                                                                                                       gCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTCTGCAACCAGCTATGACGATGTCAAAGCTTATACCCCGGCCTGGGCCGAg                  >  5:266823/1‑145 (MQ=255)
                                                                                                                        cAACTGGATGATGGCAGCACCGCCCTGGTGACCACCGTTTTTGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCATGGGCCGAGCAg               <  1:66675/147‑1 (MQ=255)
                                                                                                                                 gATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGATTACCGGCg     >  8:252475/1‑148 (MQ=255)
                                                                                                                                    ggCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGATTACCGGCGttt  >  1:69322/1‑148 (MQ=255)
                                                                                                                                         |                                                                                                                                              
GGTAGCCAGGATGAGATCGCTGAGGTAGGCTTCCCGTACTTTGGTGGCGACGGCACGGAACACTTCAACAAAGTGGAACTGGAAAACGTGCTGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGATTACCGGCGTTT  >  NC_000913/1281127‑1281406

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: ‑