Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,606,758 T→C V220A (GTA→GCA)  lsrF → 3‑hydroxy‑2,4‑pentadione 5‑phosphate thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,606,7580TC95.7% 148.6 / ‑4.6 46V220A (GTA→GCA) lsrF3‑hydroxy‑2,4‑pentadione 5‑phosphate thiolase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (20/24);  total (21/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

GCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGTACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGA  >  NC_000913/1606623‑1606892
                                                                                                                                       |                                                                                                                                      
gCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTAGTATGTCGAGAAAAGTTGTGAACGGATTGTTGCCGGATGTCCGGCAccc                                                                                                                                    >  2:221073/1‑140 (MQ=255)
     ggCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGt                                                                                                                               >  1:117748/1‑140 (MQ=255)
          aGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTg                                                                                                                          >  2:814176/1‑140 (MQ=255)
          aGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTg                                                                                                                          >  2:789153/1‑140 (MQ=255)
                ggTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCg                                                                                                                    >  2:700130/1‑140 (MQ=255)
                      cgATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTaaaa                                                                                                              >  1:490820/1‑140 (MQ=255)
                         tCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAAt                                                                                                           >  1:317192/1‑140 (MQ=255)
                             cGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGTACCCATTGTTTTTGCTTGCGGTAAAAAATTAcc                                                                                                       >  2:355846/1‑140 (MQ=255)
                                 aTTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAg                                                                                                   >  1:51221/1‑140 (MQ=255)
                                     ctcgctcgCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGcgcg                                                                                               <  2:51221/140‑1 (MQ=255)
                                       cgctcgcGACTCGAATCGCCGCTGAGATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGa                                                                                              >  2:559915/1‑139 (MQ=255)
                                       cgctcgcGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGa                                                                                              <  1:559915/139‑1 (MQ=255)
                                       cgctcgcGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAg                                                                                             >  2:91356/1‑140 (MQ=255)
                                           cgcgACTCGAATCGCCGTTGAATTGGGGGCGTAAATTATTAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGgcgc                                                                                         <  2:209678/140‑1 (MQ=255)
                                             cgACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTg                                                                                       >  2:739158/1‑140 (MQ=255)
                                               aCTCGAATCGCCGCTTAAATGGTGGAGCATATTATCATAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGa                                                                                     <  1:221073/140‑1 (MQ=255)
                                                    aaTCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAAtgt                                                                                >  2:485595/1‑140 (MQ=255)
                                                         ccgcTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTgg                                                                           <  1:814176/140‑1 (MQ=255)
                                                           gcTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCa                                                                         <  1:91356/140‑1 (MQ=255)
                                                           gcTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCa                                                                         <  1:700130/140‑1 (MQ=255)
                                                           gcTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCa                                                                         <  1:2008/140‑1 (MQ=255)
                                                             tGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAgg                                                                       <  1:789153/140‑1 (MQ=255)
                                                             tGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAgg                                                                       <  1:485595/140‑1 (MQ=255)
                                                               aaaTGGGGGCGCAAAATATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGTGCGAGGCGCTGGAAATGTGCTGGCAGGCt                                                                     <  1:355846/140‑1 (MQ=255)
                                                                    ggggCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTa                                                                    >  1:179551/1‑136 (MQ=255)
                                                                    ggggCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTa                                                                    <  2:179551/136‑1 (MQ=255)
                                                                     gggCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGTACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGAt                                                               <  1:596319/140‑1 (MQ=255)
                                                                         gcAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCa                                                                         <  1:65591/126‑1 (MQ=255)
                                                                         gcAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCa                                                                         >  2:65591/1‑126 (MQ=255)
                                                                                  cAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTg                                                  <  2:490820/140‑1 (MQ=255)
                                                                                   aaaaCCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTgg                                                 <  2:148103/140‑1 (MQ=255)
                                                                                     aaCCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGgtg                                               <  2:317192/140‑1 (MQ=255)
                                                                                     aaCCTATTATGTCGAAAAAGGTTTTGAACGGATAGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGgtg                                               <  2:348504/140‑1 (MQ=255)
                                                                                       ccTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGtgtg                                             >  1:143615/1‑140 (MQ=255)
                                                                                         tattatGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGa                                           >  1:35490/1‑140 (MQ=255)
                                                                                          attatGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGat                                          >  1:811764/1‑140 (MQ=255)
                                                                                           ttatGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGata                                         <  1:460731/140‑1 (MQ=255)
                                                                                                tCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATgggg                                    <  2:117748/140‑1 (MQ=255)
                                                                                                      aaGGTTTTGAACGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAat                              <  2:143615/140‑1 (MQ=255)
                                                                                                                aCGGATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGt                    <  1:739158/140‑1 (MQ=255)
                                                                                                                   gATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCg                                  >  1:649243/1‑123 (MQ=255)
                                                                                                                   gATTGTTGCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCg                                  <  2:649243/123‑1 (MQ=255)
                                                                                                                        ttgCCGGATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCAt            <  2:476101/140‑1 (MQ=255)
                                                                                                                              gATGTCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGGTATGGGGGGTAATAGTTTCCAGTCTGACCATCCGGTg      >  1:247061/1‑140 (MQ=255)
                                                                                                                               aTGTCCGGCACCCATTGATATTCCTGGCGGTAAAACATTCCCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTgg     <  2:233869/140‑1 (MQ=255)
                                                                                                                                  tCCGGCACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATGTTGCAGTCTGACCATCCGGTGGCga  >  1:233869/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GCGTGGGCAAAGATATGGTGCGCGATCAGCGTTATTTCTCGCTCGCGACTCGAATCGCCGCTGAAATGGGGGCGCAAATTATCAAAACCTATTATGTCGAAAAAGGTTTTGAACGGATTGTTGCCGGATGTCCGGTACCCATTGTTATTGCTGGCGGTAAAAAATTACCGGAGCGCGAGGCGCTGGAAATGTGCTGGCAGGCTATCGATCAGGGCGCTTCTGGTGTGGATATGGGGCGTAATATTTTCCAGTCTGACCATCCGGTGGCGA  >  NC_000913/1606623‑1606892

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