Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,835,582 A→G 7.7% H23R (CAT→CGT)  ynjA → carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,835,5820AG 7.7% 170.4 / 3.9 52H23R (CAT→CGT) ynjAcarboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base G (2/2);  total (27/25)
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.87e-01

TTTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCA  >  NC_000913/1835438‑1835722
                                                                                                                                                |                                                                                                                                            
tttATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTa                                                                                                                                           >  2:698504/1‑148 (MQ=255)
 ttATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTAc                                                                                                                                          >  1:568703/1‑148 (MQ=255)
        gTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAgg                                                                                                                                   <  2:568703/148‑1 (MQ=255)
           tGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGc                                                                                                                                <  1:698504/148‑1 (MQ=255)
                         tgATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGc                                                                                                                  >  4:303756/1‑148 (MQ=255)
                                                tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                     >  4:701173/1‑106 (MQ=255)
                                                tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                     >  2:930665/1‑106 (MQ=255)
                                                tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                     <  3:701173/106‑1 (MQ=255)
                                                tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                     <  1:930665/106‑1 (MQ=255)
                                                tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCAtt                                                                                                                     <  2:1063068/122‑1 (MQ=255)
                                                tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCAtt                                                                                                                     >  1:1063068/1‑122 (MQ=255)
                                                    ttGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttt                                                                                       <  3:303756/148‑1 (MQ=255)
                                                            cAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGc                                                                                                        >  2:657065/1‑123 (MQ=255)
                                                            cAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGc                                                                                                        <  1:657065/123‑1 (MQ=255)
                                                                    gAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttat                                                                                           >  1:346643/1‑128 (MQ=255)
                                                                    gAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttat                                                                                           <  2:346643/128‑1 (MQ=255)
                                                                              tGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtggctg                                                             >  3:576022/1‑148 (MQ=255)
                                                                                   ttACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCgg                                                                                                      <  3:533701/102‑1 (MQ=255)
                                                                                   ttACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCgg                                                                                                      >  4:533701/1‑102 (MQ=255)
                                                                                     aCCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAg                                                      <  4:576022/148‑1 (MQ=255)
                                                                                      ccgccgCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTa                                                                                                                              <  2:758859/75‑1 (MQ=255)
                                                                                      ccgccgCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTa                                                                                                                              >  1:758859/1‑75 (MQ=255)
                                                                                             ttAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatc                                                                                          <  1:864984/104‑1 (MQ=255)
                                                                                             ttAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatc                                                                                          >  2:864984/1‑104 (MQ=255)
                                                                                                           ttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctt                                                                                        >  2:850880/1‑92 (MQ=255)
                                                                                                           ttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctt                                                                                        <  1:850880/92‑1 (MQ=255)
                                                                                                              aTTTTACGTCCACCGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCt                                            >  1:845983/1‑133 (MQ=255)
                                                                                                              aTTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCt                                            <  2:845983/133‑1 (MQ=255)
                                                                                                                  tACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttttatc                                                                                   <  2:714016/90‑1 (MQ=255)
                                                                                                                  tACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttttatc                                                                                   >  1:714016/1‑90 (MQ=255)
                                                                                                                      tCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATcc                                                                                                  <  4:484352/71‑1 (MQ=255)
                                                                                                                      tCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATcc                                                                                                  >  3:484352/1‑71 (MQ=255)
                                                                                                                      tCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCg                                                                           >  4:972739/1‑94 (MQ=255)
                                                                                                                      tCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCg                                                                           <  3:972739/94‑1 (MQ=255)
                                                                                                                          caGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctttta                                                                                     <  1:559547/80‑1 (MQ=255)
                                                                                                                          caGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctttta                                                                                     >  2:559547/1‑80 (MQ=255)
                                                                                                                                tGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                   <  1:810686/92‑1 (MQ=255)
                                                                                                                                tGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                   >  4:194233/1‑92 (MQ=255)
                                                                                                                                tGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                   >  2:810686/1‑92 (MQ=255)
                                                                                                                                tGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                   <  3:194233/92‑1 (MQ=255)
                                                                                                                                       aTCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAgcgcgcg    >  1:985388/1‑148 (MQ=255)
                                                                                                                                         cgtcgCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCa  >  1:953327/1‑148 (MQ=255)
                                                                                                                                          ctcgCCGTAATGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgtt                                                                  >  1:572700/2‑83 (MQ=255)
                                                                                                                                          ctcgCCGTAATGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgtt                                                                  <  2:572700/82‑1 (MQ=255)
                                                                                                                                           tcgCCGTAATGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAg          >  2:38005/1‑138 (MQ=255)
                                                                                                                                           tcgCCGTAATGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAg          <  1:38005/138‑1 (MQ=255)
                                                                                                                                               cATTACGGCGTGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCg                  >  3:651665/1‑126 (MQ=255)
                                                                                                                                               cATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGctggctgg                                                        >  1:848806/1‑88 (MQ=255)
                                                                                                                                               cATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGctggctgg                                                        <  2:848806/88‑1 (MQ=255)
                                                                                                                                               cATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCg                  <  4:651665/126‑1 (MQ=255)
                                                                                                                                               cATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTc            >  3:785130/1‑132 (MQ=255)
                                                                                                                                               cATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTc            <  4:785130/132‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
TTTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCA  >  NC_000913/1835438‑1835722

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