Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,832,517 G→A 12.9% R15H (CGT→CAT)  norV → anaerobic nitric oxide reductase flavorubredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,832,5170GA12.9% 75.7 / 2.5 31R15H (CGT→CAT) norVanaerobic nitric oxide reductase flavorubredoxin
Reads supporting (aligned to +/- strand):  ref base G (11/16);  new base A (2/2);  total (13/18)
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.18e-01

AATTTTATTTATAGAGTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCGAAGAAAAAAACGTGCTGATCGACACCGTCGACCATAAATTCAGCCGCGAATTTGTGCAGAACCTG  >  NC_000913/2832373‑2832662
                                                                                                                                                |                                                                                                                                                 
aaTTTTATTTATAGAGTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATACCATTCATTGGGTTGGTCAACATGAc                                                                                                                                               >  1:453581/1‑149 (MQ=255)
aaTTTTATTTATAGAGTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGAc                                                                                                                                               <  2:299858/149‑1 (MQ=255)
aaTTTTATTTATAGAGTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACATGAc                                                                                                                                               >  1:84602/1‑149 (MQ=255)
   tttatttaTAGAGTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTCATATGTTTTTTGTGGTGAATAATAACATTCATTGGGTTGGTCAACGTGACTgg                                                                                                                                            >  1:644368/1‑149 (MQ=255)
    ttatttatAGAGTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACAATCATTGGGTTGGTCAACGTGACTggg                                                                                                                                           >  1:260841/1‑149 (MQ=255)
          ataGAGTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGc                                                                                                                                     >  1:779750/1‑149 (MQ=255)
            ggagTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCgg                                                                                                                                   >  2:595592/2‑148 (MQ=255)
              agTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGa                                                                                                                                 >  2:679804/1‑149 (MQ=255)
                           aGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTcacggcacg                                                                                                                    >  1:237052/1‑149 (MQ=255)
                           aGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTcacggcacg                                                                                                                    <  1:281718/149‑1 (MQ=255)
                              tAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGaa                                                                                                                 <  2:720917/149‑1 (MQ=255)
                                 aaaaCTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGGGATTTTCACGGCACGGAatat                                                                                                              >  2:200139/1‑149 (MQ=255)
                                  aaaCTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAAtata                                                                                                             >  1:428070/1‑149 (MQ=255)
                                     cTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATaaac                                                                                                          <  2:208793/149‑2 (MQ=255)
                                          aCGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTg                                                                                                     <  2:56479/149‑1 (MQ=255)
                                          aCGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATGCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTg                                                                                                     <  1:522671/149‑1 (MQ=255)
                                                     aaTTAGCAAGATATGTTTTTAGAACACGCGAAATTAATTGAGGTTGGTATGTCTTTTGTTGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGcagcag                                                                                          <  1:722631/149‑1 (MQ=255)
                                                     aaTTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATCCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGcagcag                                                                                          <  1:692321/149‑1 (MQ=255)
                                                      aTTAGCAAGACATCTTTTTAGAACACGCTGAATCAATTGATTTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCagcagc                                                                                         >  1:188677/1‑149 (MQ=255)
                                                      aTTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCggcagc                                                                                         >  2:156126/1‑149 (MQ=255)
                                                      aTTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCagcagc                                                                                         >  2:380351/1‑149 (MQ=255)
                                                                  tCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTag                                                                             <  2:398632/149‑2 (MQ=255)
                                                                  tCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTAc                                                                             <  1:592245/149‑1 (MQ=255)
                                                                  gcTTTTGAGAACACGCGGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAAGAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTAc                                                                             <  1:397477/148‑1 (MQ=255)
                                                                             cacGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTTACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCga                                                                  <  1:199144/149‑1 (MQ=255)
                                                                                    aaTACATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCGAAGaaaaa                                                           <  1:57854/149‑1 (MQ=255)
                                                                                    aaTAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGGGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCGAAGaaaaa                                                           <  1:419707/149‑1 (MQ=255)
                                                                                                                                    gggTTGGTCAACATTACTGGGAAGGGGGTCATTTTCACGGGACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCGAAAAAAAAAACGTGCTGATCGACACCGTCGACCATAAATTCAGCCGCGAATTtgtg           <  2:84602/149‑1 (MQ=255)
                                                                                                                                    gggTTGGTCAACATGACTGGGAAGGGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCGAAGAAAAAAACGTGCTGATCGACACCGTCGACCATAAATTCAGCCGCGAATTtgtg           <  2:453581/149‑1 (MQ=255)
                                                                                                                                       ttGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCGAAGAAAAAAACGTGCTGATCGACACCGTCGACCATAAATTCAGCCGCGAATTTGTGCAg        <  1:200139/149‑1 (MQ=255)
                                                                                                                                             aaCGTGATTGGGAAGTGCGTTATTTTCACGGCACGGAATATAAAACGCCGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCGAAGAAAAAAACGTGCTGATCGACACCGTCGACCATAAATTCAGCCGCGAATTTGTGCAGAACCTg  <  1:14805/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AATTTTATTTATAGAGTAAAAACAATCAGATAAAAAACTGGCACGCAATCTGCAATTAGCAAGACATCTTTTTAGAACACGCTGAATAAATTGAGGTTGCTATGTCTATTGTGGTGAAAAATAACATTCATTGGGTTGGTCAACGTGACTGGGAAGTGCGTGATTTTCACGGCACGGAATATAAAACGCTGCGCGGCAGCAGCTACAATAGCTACCTCATCCGCGAAGAAAAAAACGTGCTGATCGACACCGTCGACCATAAATTCAGCCGCGAATTTGTGCAGAACCTG  >  NC_000913/2832373‑2832662

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