Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 375,307 C→T intergenic (+426/‑152) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,3070CT100.0% 94.5 / NA 32intergenic (+426/‑152)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/15);  total (17/15)

AATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTC  >  NC_000913/375178‑375441
                                                                                                                                 |                                                                                                                                      
aaTTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCa                                                                                                                               >  2:303442/1‑139 (MQ=25)
  tttGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAgg                                                                                                                             <  2:997376/139‑1 (MQ=37)
            tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGcaca                                                                                                                   <  1:303442/139‑1 (MQ=25)
                     gcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCg                                                                                                          >  1:35638/1‑139 (MQ=37)
                        ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCg                                                                                                                                    <  2:514714/110‑1 (MQ=17)
                        ttAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCg                                                                                                                                    >  1:514714/1‑110 (MQ=18)
                                         aTCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                                                                                      <  2:35638/139‑1 (MQ=255)
                                           cTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCa                                                                                                                                  <  2:787284/93‑1 (MQ=18)
                                           cTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCa                                                                                                                                  >  1:787284/1‑93 (MQ=18)
                                                 aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                                                                 >  1:153894/1‑120 (MQ=255)
                                                 aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                                                                                                 <  2:153894/120‑1 (MQ=255)
                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCg                                                                                                    >  1:3260/1‑110 (MQ=38)
                                                        ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCg                                                                                                    <  2:3260/110‑1 (MQ=25)
                                                                ggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAgtgt                                                                                                                                      <  1:1302368/68‑1 (MQ=25)
                                                                ggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAgtgt                                                                                                                                      >  2:1302368/1‑68 (MQ=25)
                                                                  cgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAg                                                             <  2:1191135/139‑1 (MQ=255)
                                                                                             ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               >  2:291904/1‑110 (MQ=255)
                                                                                             ggCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  1:291904/110‑1 (MQ=255)
                                                                                                              gCCTGATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTcgc                 >  2:945205/1‑139 (MQ=255)
                                                                                                                  gATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               >  1:636851/1‑89 (MQ=255)
                                                                                                                  gATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGc                                                               <  2:636851/89‑1 (MQ=255)
                                                                                                                  gATAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGAtt             >  2:355611/1‑139 (MQ=255)
                                                                                                                   aTAAGACGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTc            >  1:1242001/1‑139 (MQ=255)
                                                                                                                        aCGCGTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTg       >  1:876198/1‑139 (MQ=255)
                                                                                                                           cgTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcc                                                                                 >  1:803886/1‑62 (MQ=255)
                                                                                                                           cgTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcc                                                                                 <  2:803886/62‑1 (MQ=255)
                                                                                                                            gTTAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACt   <  1:262365/139‑1 (MQ=255)
                                                                                                                             ttAGTGTCGCATCAGGCATCTGCACACGACTGTCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCaa                             <  1:829195/112‑1 (MQ=255)
                                                                                                                             ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                                                          <  2:436801/83‑1 (MQ=255)
                                                                                                                             ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCa                                                          >  1:436801/1‑83 (MQ=255)
                                                                                                                             ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCaa                             >  2:829195/1‑112 (MQ=255)
                                                                                                                             ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGt                    <  2:1139513/121‑1 (MQ=255)
                                                                                                                             ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGt                    >  1:1139513/1‑121 (MQ=255)
                                                                                                                             ttAGTGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTc  >  2:543515/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                      
AATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGCCCGTAGGCAGTCATTAAATAGTTCTGATTACGCGGCTGGCAATGTATCAGTCGCGATTCACCTGCACTC  >  NC_000913/375178‑375441

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