Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,067 G→A 8.7% P73P (CCG→CCA)  narI → nitrate reductase 1, gamma (cytochrome b(NR)) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,0670GA 8.7% 141.5 / 3.5 46P73P (CCG→CCA) narInitrate reductase 1, gamma (cytochrome b(NR)) subunit
Reads supporting (aligned to +/- strand):  ref base G (23/19);  new base A (2/2);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACGCGTGCGTGCAACC  >  NC_000913/1285923‑1286211
                                                                                                                                                |                                                                                                                                                
gctgCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccgc                                                                                                                                                 <  7:130195/146‑1 (MQ=255)
          gACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCACACTGGATg                                                                                                                                         <  1:160302/144‑1 (MQ=255)
          gACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCACACTGGATg                                                                                                                                         <  7:356371/144‑1 (MQ=255)
          gACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCACACTGGATg                                                                                                                                         >  8:356371/1‑144 (MQ=255)
          gACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCACACTGGATg                                                                                                                                         >  2:160302/1‑144 (MQ=255)
           aCTACGGGCAGTACACGTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGTTGAACCTGGCGTCGAAGCTGTTCCATATCGGGTTTCTTGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGAt                                                                                                                                          <  3:313676/142‑1 (MQ=255)
            cTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCTAAGGGATGAACGTGGCGTCGAACATGTTCCATATCGGGATTCTGGGGAGTATTGTCGGTCCCTTCTTCGGTATGCTGACGCGGCACTGGATg                                                                                                                                         >  4:313676/1‑142 (MQ=38)
               cGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAg                                                                                                                                  <  7:48784/146‑1 (MQ=255)
                       acacCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTCTGCTGCCGCCGCACTGGATGTATGAAGCCTgtttgtg                                                                                                                        >  8:357447/1‑142 (MQ=255)
                       acacCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAg                                                                                                                                  <  8:274880/138‑1 (MQ=255)
                       acacCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAg                                                                                                                                  >  7:274880/1‑138 (MQ=255)
                       acacCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGcc                                                                                                                        >  2:266138/1‑148 (MQ=255)
                             ggCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGTTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTg                                                                                                                   >  5:51418/1‑147 (MQ=255)
                             ggCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTg                                                                                                                   >  7:131502/1‑147 (MQ=255)
                                          cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaaa                                                                                                     <  8:196283/148‑1 (MQ=255)
                                          cagccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGaa                                                                                                      <  8:131502/147‑1 (MQ=255)
                                             ccagATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGTTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATg                                                                                                  <  6:51418/148‑1 (MQ=255)
                                                gATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGc                                                                                                >  3:38676/1‑147 (MQ=255)
                                                gATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCa                                                                                               <  1:266138/148‑1 (MQ=255)
                                                 aTGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCaa                                                                                              >  1:142397/1‑148 (MQ=255)
                                                                           ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTg                                                                    >  8:37689/1‑148 (MQ=255)
                                                                                cGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCtgtg                                                               >  4:307347/1‑148 (MQ=255)
                                                                                    ccTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTg                                                           <  2:142397/148‑1 (MQ=255)
                                                                                                    aTTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGt                                           >  4:115699/1‑148 (MQ=255)
                                                                                                    aTTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGt                                           >  7:86328/1‑148 (MQ=255)
                                                                                                     ttCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCCGTGTCTGATTGGCGGCGTGCTGtt                                          <  7:37689/148‑1 (MQ=255)
                                                                                                      tCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCtgtg                                                               >  7:40663/1‑126 (MQ=255)
                                                                                                      tCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCtgtg                                                               <  8:40663/126‑1 (MQ=255)
                                                                                                       cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCt                                                                  <  5:231894/122‑1 (MQ=255)
                                                                                                       cTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCt                                                                  >  6:231894/1‑122 (MQ=255)
                                                                                                          gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCt                                       <  4:38676/146‑1 (MQ=255)
                                                                                                             aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAc                                  >  8:280660/1‑148 (MQ=255)
                                                                                                              tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcg                                 >  7:153/1‑148 (MQ=255)
                                                                                                                tttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcgtc                               >  5:318949/1‑148 (MQ=255)
                                                                                                                tttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcgtc                               >  3:331721/1‑148 (MQ=255)
                                                                                                                   gTCGGTCACTTCTTCGGTATGCTGACGCCGCCCTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACgtcgtc                            >  1:94806/1‑148 (MQ=255)
                                                                                                                   gTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtg                                                                                     <  4:153661/91‑1 (MQ=255)
                                                                                                                   gTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCtg                                                                                     >  3:153661/1‑91 (MQ=255)
                                                                                                                       gTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGtt                        <  8:86328/148‑1 (MQ=255)
                                                                                                                       gTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGtt                        <  3:307347/148‑1 (MQ=255)
                                                                                                                          aCTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAg                     <  3:115699/148‑1 (MQ=255)
                                                                                                                             tcttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcgtcg                              >  2:244198/1‑136 (MQ=255)
                                                                                                                             tcttcGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAAcgtcg                              <  1:244198/136‑1 (MQ=255)
                                                                                                                                      tGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACgcgtgc          <  2:94806/147‑1 (MQ=255)
                                                                                                                                       gCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGaa                                    <  2:348631/120‑1 (MQ=255)
                                                                                                                                       gCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGaa                                    >  1:348631/1‑120 (MQ=255)
                                                                                                                                       gCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACgcgtgcgt        >  3:69934/1‑148 (MQ=255)
                                                                                                                                             gccgcACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACGCGTGCGTGCAAcc  <  8:153/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCTGCGTTATGACTACGGGCAGTACACCTGGCGCGCGGCGTCCAGCCAGATGCTGGATCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAACAGAAAATGGCAATGTTTGCTGGTGGTGCCAGCGGCGTGCTGTGTCTGATTGGCGGCGTGCTGTTGCTGAAACGTCGTCTGTTCAGCCCACGCGTGCGTGCAACC  >  NC_000913/1285923‑1286211

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