Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 372,113 T→C intergenic (+8/‑2) mhpC → / → mhpD 2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase/2‑hydroxypentadienoate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913372,1130TC97.7% 148.6 / ‑4.4 44intergenic (+8/‑2)mhpC/mhpD2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase/2‑hydroxypentadienoate hydratase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/20);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCATTCAGCACATAAAT  >  NC_000913/371978‑372249
                                                                                                                                       |                                                                                                                                        
gtctgcgtctgcTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATg                                                                                                                                      >  2:370448/1‑140 (MQ=255)
   tgcgtctgcTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACg                                                                                                                                   >  2:596236/1‑140 (MQ=255)
   tgcgtctgcTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACg                                                                                                                                   <  1:734832/140‑1 (MQ=255)
        ctgctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCa                                                                                                                              <  1:622361/140‑1 (MQ=255)
          gctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGAc                                                                                                                                    <  1:790539/132‑1 (MQ=255)
          gctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGAc                                                                                                                                    >  2:790539/1‑132 (MQ=255)
          gctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATa                                                                                                                            <  1:787706/140‑1 (MQ=255)
                cGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGaa                                                                                                                                 <  1:73066/129‑1 (MQ=255)
                cGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGaa                                                                                                                                 >  2:73066/1‑129 (MQ=255)
                        ccGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTg                                                                                                              >  2:690786/1‑140 (MQ=255)
                        ccGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTg                                                                                                              >  2:756447/1‑140 (MQ=255)
                                      cATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTAcg                                                                                                >  1:53414/1‑140 (MQ=255)
                                                gcgACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCa                                                                                      >  2:638172/1‑140 (MQ=255)
                                                gcgACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCa                                                                                      >  2:309510/1‑140 (MQ=255)
                                                       tgGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAgg                                                                               >  2:505897/1‑140 (MQ=255)
                                                                ggCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGa                                                                      <  2:169803/140‑1 (MQ=255)
                                                                   gcAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTg                                                                   >  1:42131/1‑140 (MQ=255)
                                                                    cAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGc                                                                  >  2:211936/1‑140 (MQ=255)
                                                                        gggAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCg                                                              <  2:352192/140‑1 (MQ=255)
                                                                               tGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGcgcg                                                       >  2:685359/1‑140 (MQ=255)
                                                                                 ccGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGAt                                                     <  2:469721/140‑1 (MQ=255)
                                                                                    aCGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTg                                                  >  2:252426/1‑140 (MQ=255)
                                                                                     cGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGa                                                 >  1:650208/1‑140 (MQ=255)
                                                                                     cGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGa                                                 <  1:505897/140‑1 (MQ=255)
                                                                                     cGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGa                                                 <  1:309510/140‑1 (MQ=255)
                                                                                     cGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGa                                                 <  1:211936/140‑1 (MQ=255)
                                                                                     cGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGa                                                 <  2:650208/140‑1 (MQ=255)
                                                                                          tcaatcaaCTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGt                                            <  1:199782/140‑1 (MQ=255)
                                                                                          tcaatcaaCTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGt                                            <  2:87808/140‑1 (MQ=255)
                                                                                               caaCTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGa                                       <  1:596236/140‑1 (MQ=255)
                                                                                               caaCTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGa                                       <  2:8296/140‑1 (MQ=255)
                                                                                                  cTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATaa                                    >  2:129533/1‑140 (MQ=255)
                                                                                                        cTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGa                                       >  2:252205/1‑131 (MQ=255)
                                                                                                        cTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGa                                       <  1:252205/131‑1 (MQ=255)
                                                                                                        cTGAATTTCCTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGa                              >  1:145838/1‑140 (MQ=255)
                                                                                                                ccTCGCACGCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCtt                      >  1:485657/1‑140 (MQ=255)
                                                                                                                        gCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCAtt              <  2:145838/140‑1 (MQ=255)
                                                                                                                        gCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCAtt              <  1:685359/140‑1 (MQ=255)
                                                                                                                        gCCCTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCAtt              <  1:129533/140‑1 (MQ=255)
                                                                                                                          ccTTAAGGAATGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCATTCa            <  2:42131/140‑1 (MQ=255)
                                                                                                                                  aaTGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCATTCAGCACATaa    >  1:151058/1‑140 (MQ=255)
                                                                                                                                    tGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCATTCAGCACATAAAt  >  2:672502/1‑140 (MQ=255)
                                                                                                                                    tGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCATTCAGCACATAAAt  >  1:78244/1‑140 (MQ=255)
                                                                                                                                    tGGCCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCATTCAGCACATAAAt  >  1:659847/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                        
GTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCGGATTTACGCCGCGCCGCAGAGCAGGGCGAAGCGATTGCACCGCTGCGCGATCTGATTGGTATCGATAACGCTGAAGCGGCTTACGCCATTCAGCACATAAAT  >  NC_000913/371978‑372249

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