Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,835,584 T→A 7.7% Y24N (TAC→AAC) ‡ ynjA → carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,835,5840TA 7.7% 182.6 / 5.1 52Y24N (TAC→AAC) ‡ynjAcarboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base A (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 = 1.00e+00

TTTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAA  >  NC_000913/1835438‑1835731
                                                                                                                                                  |                                                                                                                                                   
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                                                                                                                                              <  3:701173/106‑1 (MQ=255)
                                                tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                              >  2:930665/1‑106 (MQ=255)
                                                tGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGa                                                                                                                                              <  1:930665/106‑1 (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                                                                                                    >  1:346643/1‑128 (MQ=255)
                                                                    gAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttat                                                                                                    <  2:346643/128‑1 (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                                                                                                                                       >  1:758859/1‑75 (MQ=255)
                                                                                      ccgccgCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTa                                                                                                                                       <  2:758859/75‑1 (MQ=255)
                                                                                             ttAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatc                                                                                                   >  2:864984/1‑104 (MQ=255)
                                                                                             ttAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatc                                                                                                   <  1:864984/104‑1 (MQ=255)
                                                                                                           ttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctt                                                                                                 <  1:850880/92‑1 (MQ=255)
                                                                                                           ttAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctt                                                                                                 >  2:850880/1‑92 (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                                                                                    <  3:972739/94‑1 (MQ=255)
                                                                                                                      tCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCg                                                                                    >  4:972739/1‑94 (MQ=255)
                                                                                                                          caGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctttta                                                                                              <  1:559547/80‑1 (MQ=255)
                                                                                                                          caGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatctttta                                                                                              >  2:559547/1‑80 (MQ=255)
                                                                                                                                tGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                            <  3:194233/92‑1 (MQ=255)
                                                                                                                                tGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                            >  4:194233/1‑92 (MQ=255)
                                                                                                                                tGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                            <  1:810686/92‑1 (MQ=255)
                                                                                                                                tGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                            >  2:810686/1‑92 (MQ=255)
                                                                                                                                       aTCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAgcgcgcg             >  1:985388/1‑148 (MQ=255)
                                                                                                                                         cgtcgCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCa           >  1:953327/1‑148 (MQ=255)
                                                                                                                                          ctcgCCGTAATGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgtt                                                                           <  2:572700/82‑1 (MQ=255)
                                                                                                                                          ctcgCCGTAATGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgtt                                                                           >  1:572700/2‑83 (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)
                                                                                                                                                  |                                                                                                                                                   
TTTATTCTGTTTGCCCTGACCGTGGTGATTTTTATGGCGAAGAAAATATGGCTTGAACGCCAGAAGAGGAATGCCTGATGGGTTTACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGCATTGCGCAAA  >  NC_000913/1835438‑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: ‑