Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,835,586 C→T 7.7% Y24Y (TAC→TAT) ‡ ynjA → carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,835,5860CT 7.7% 161.9 / 2.0 51Y24Y (TAC→TAT) ‡ynjAcarboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base T (2/2);  total (26/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

TTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAA  >  NC_000913/1835439‑1835731
                                                                                                                                                   |                                                                                                                                                 
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                                                                                                                                              <  3:701173/106‑1 (MQ=255)
                                               tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                              >  2:930665/1‑106 (MQ=255)
                                               tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                              <  1:930665/106‑1 (MQ=255)
                                               tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                              >  4:701173/1‑106 (MQ=255)
                                               tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCAtt                                                                                                                              >  1:1063068/1‑122 (MQ=255)
                                               tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCAtt                                                                                                                              <  2:1063068/122‑1 (MQ=255)
                                                   ttGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttt                                                                                                <  3:303756/148‑1 (MQ=255)
                                                           cAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGc                                                                                                                 <  1:657065/123‑1 (MQ=255)
                                                           cAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGc                                                                                                                 >  2:657065/1‑123 (MQ=255)
                                                                   gAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttat                                                                                                    <  2:346643/128‑1 (MQ=255)
                                                                   gAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttat                                                                                                    >  1:346643/1‑128 (MQ=255)
                                                                             tGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtggctg                                                                      >  3:576022/1‑148 (MQ=255)
                                                                                  ttACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCgg                                                                                                               >  4:533701/1‑102 (MQ=255)
                                                                                  ttACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCgg                                                                                                               <  3:533701/102‑1 (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)
                                                                                                                                                 tACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCaaa  >  2:1094022/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAA  >  NC_000913/1835439‑1835731

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