Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 370,186 C→T 7.6% Q526* (CAA→TAA)  mhpA → 3‑(3‑hydroxyphenyl)propionate hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913370,1860CT 7.6% 158.7 / 2.9 53Q526* (CAA→TAA) mhpA3‑(3‑hydroxyphenyl)propionate hydroxylase
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base T (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

ATTCATACCGCACAGGATAACCACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGAC  >  NC_000913/370045‑370330
                                                                                                                                             |                                                                                                                                                
aTTCATACCGCACAGGATAACCACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCt                                                                                                                                           <  4:99211/149‑1 (MQ=255)
       ccGCACAGGATAACCACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAACCCCTGGGc                                                                                                                                       >  5:297134/1‑146 (MQ=255)
            caGGATAACCACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCt                                                                                                                               <  3:379518/149‑1 (MQ=255)
                aTAACCACGACGGCGTACCACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAAt                                                                                                                           <  2:264670/149‑1 (MQ=255)
                  acccACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCAGCAAACCCTGGGCAAGACCCTGAATa                                                                                                                          >  1:103496/3‑148 (MQ=255)
                  aaCCACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAg                                                                                                                                    >  7:151942/1‑138 (MQ=255)
                  aaCCACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAg                                                                                                                                    <  8:151942/138‑1 (MQ=255)
                    ccACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAAc                                                                                                                       >  3:57419/1‑149 (MQ=255)
                     cACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACt                                                                                                                      >  1:4116/1‑149 (MQ=255)
                      acgacgGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTg                                                                                                                     <  3:209224/149‑1 (MQ=255)
                                    cgcgTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGAc                                                                                                       >  8:186225/1‑149 (MQ=255)
                                            cGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGAccc                                                                                               <  7:49932/149‑1 (MQ=255)
                                              aTACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGc                                                                                             <  2:103496/149‑1 (MQ=255)
                                                aCGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGTAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGccc                                                                                           >  7:58501/1‑149 (MQ=255)
                                                aCGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGTAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGccc                                                                                           >  7:51873/1‑149 (MQ=255)
                                                  gCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCt                                                                                          >  4:185473/1‑148 (MQ=255)
                                                  gCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAAACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTg                                                                                         <  6:330496/149‑1 (MQ=255)
                                                     aGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATg                                                                                      >  4:151234/1‑149 (MQ=255)
                                                      gcgaccttGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGc                                                                                     <  7:422552/142‑1 (MQ=255)
                                                           ccTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACg                                                                                <  5:110215/149‑1 (MQ=255)
                                                           ccTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACg                                                                                <  4:57419/149‑1 (MQ=255)
                                                                 tAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTg                                                                          >  4:202340/1‑149 (MQ=255)
                                                                   gCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGt                                                                         <  3:77580/148‑1 (MQ=255)
                                                                   gCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTc                                                                        <  3:202340/149‑1 (MQ=255)
                                                                    cTGGTTCTCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGGCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCg                                                                       <  4:22753/149‑1 (MQ=255)
                                                                          cgcgCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGTAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAg                                                                  <  8:58501/148‑1 (MQ=255)
                                                                          cgcgCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGTAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAg                                                                  <  8:51873/148‑1 (MQ=255)
                                                                           gcgcAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGGCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGt                                                                <  7:186225/149‑1 (MQ=255)
                                                                            cgcAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTc                                                                                                                >  7:177154/1‑100 (MQ=255)
                                                                            cgcAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTc                                                                                                                <  8:177154/100‑1 (MQ=255)
                                                                                  acaATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGCCGTTTCTGTCGAAAAGGTAGCCTGa                                                         >  2:274808/1‑149 (MQ=255)
                                                                                          tCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGata                                                       >  1:11858/1‑143 (MQ=255)
                                                                                          tCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGata                                                       <  2:11858/143‑1 (MQ=255)
                                                                                           cGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCAc                                                  <  7:302692/147‑1 (MQ=255)
                                                                                                gtggtgATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATc                                           >  1:4516/1‑149 (MQ=255)
                                                                                                   gtgATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTc                                        >  8:163396/1‑149 (MQ=255)
                                                                                                    tgATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAACCTGGCGGCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCa                                       <  6:297134/149‑1 (MQ=255)
                                                                                                    tgATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCa                                       <  2:4116/149‑1 (MQ=255)
                                                                                                      aTGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCAc                                      <  3:151234/148‑1 (MQ=255)
                                                                                                           cccGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCtt                                >  4:392272/1‑149 (MQ=255)
                                                                                                            ccGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCa                                       >  5:269709/1‑141 (MQ=255)
                                                                                                            ccGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCa                                       <  6:269709/141‑1 (MQ=255)
                                                                                                                   gCTGTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTcccccgc                        >  2:340026/1‑145 (MQ=255)
                                                                                                                       tgttgCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTcgcc                     <  2:4516/148‑1 (MQ=255)
                                                                                                                           gccgccACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATAACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCtgg                >  3:242219/1‑149 (MQ=255)
                                                                                                                            ccgcctccgcgaTTCCGAAAACGATGGGCATGACCCTGTATATACTGGTGTCGGTGATGACGCTGACCCTACGTTATGCCTATTTTTCTGTCGATTATGTTGCCTGATATGAACGCTTATGTTCTCTGTCTTTCCCACTCGCCGCtgg                <  4:373721/148‑1 (MQ=255)
                                                                                                                                     ccaTTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGt      <  3:392272/149‑1 (MQ=255)
                                                                                                                                       aTTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTc     >  7:205607/1‑148 (MQ=255)
                                                                                                                                        ttCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGa   >  8:401106/1‑149 (MQ=255)
                                                                                                                                        ttCCCCAAACCCTGGGCAAGACCCTCAATAAACTGGCGTCGGTGATGACGCTCACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCg    >  8:23737/1‑148 (MQ=255)
                                                                                                                                         tCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGAc  >  8:158076/1‑149 (MQ=255)
                                                                                                                                            gCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTggg            >  2:140313/1‑136 (MQ=255)
                                                                                                                                            gCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTggg            <  1:140313/136‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
ATTCATACCGCACAGGATAACCACGACGGCGTACTACGCGTGGGCGATACGCAAGGTCGCCTGCGTAGCTGGTTCGCGCAACACAATGCTTCGCTGGTGGTGATGCGCCCGGATCGCTTTGTTGCCGCCACCGCCATTCCGCAAACCCTGGGCAAGACCCTGAATAAACTGGCGTCGGTGATGACGCTGACCCGCCCTGATGCCGACGTTTCTGTCGAAAAGGTAGCCTGATATGCACGCTTATCTTCACTGTCTTTCCCACTCGCCGCTGGTGGGGTATGTCGAC  >  NC_000913/370045‑370330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑