Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,544,956 T→A 15.6% N284Y (AAC→TAC)  hyuC ← bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,544,9560TA15.6% 101.7 / 12.2 45N284Y (AAC→TAC) hyuCbifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base A (5/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GGCCGTGGGCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCGGCATG  >  NC_002947/4544810‑4545103
                                                                                                                                                   |                                                                                                                                                   
ggCCGTGGGCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTAc                                                                                                                                                   >  4:135504/1‑149 (MQ=255)
ggCCGTGGGCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTAc                                                                                                                                                   >  4:115156/1‑149 (MQ=255)
 gCCGTGGGCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTAcg                                                                                                                                                  <  7:204537/149‑1 (MQ=255)
        gCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGccc                                                                                                                                           <  3:135504/149‑1 (MQ=255)
         cACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGa                                                                                                                                               <  1:113176/144‑1 (MQ=255)
         cACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGa                                                                                                                                               >  2:113176/1‑144 (MQ=255)
          aCCAACTCATTGCTCAGGCCTTTCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCATGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCgg                                                                                                                                         <  6:188371/149‑1 (MQ=255)
          aCCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCgg                                                                                                                                         <  8:52700/149‑1 (MQ=255)
            cAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGt                                                                                                                                       <  5:169521/149‑1 (MQ=255)
                tCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGC‑TCGCGTAACGTGATCCCgggcgaag                                                                                                                                   >  8:44318/1‑144 (MQ=255)
                tCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCACCCGGGC‑TCGCGTAACGTGATCCCgggc                                                                                                                                       <  8:222003/145‑2 (MQ=255)
                tCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCACCCGGGC‑TCGCGTAACGTGATCCCgggc                                                                                                                                       >  7:222003/1‑144 (MQ=255)
                      cTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGTGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAg                                                                                                                             <  5:114030/149‑1 (MQ=255)
                       tCAGGCCATGCTTGGCACAGGACGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCCCGTTACGCGAGCCCGGGTACGCCTGAAgg                                                                                                                            >  8:131682/1‑149 (MQ=255)
                            ccATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAAcc                                                                                                                       <  6:250107/149‑1 (MQ=255)
                                    ggCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGc                                                                                                               <  8:171354/149‑1 (MQ=255)
                                              gCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTg                                                                                                     <  1:35373/149‑1 (MQ=255)
                                               ccTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTgg                                                                                                    <  1:253228/149‑1 (MQ=255)
                                                cTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTggg                                                                                                   >  2:267003/1‑149 (MQ=255)
                                                       aCCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCtggtg                                                                                            <  2:92768/149‑1 (MQ=255)
                                                       aCCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCtggtg                                                                                            <  1:267003/149‑1 (MQ=255)
                                                                     ggCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGc                                                                              >  3:129254/1‑149 (MQ=255)
                                                                         atcatcGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGtt                                                                          <  2:152628/149‑1 (MQ=255)
                                                                                 aTTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCct                                                                  >  1:249136/1‑149 (MQ=255)
                                                                                   tCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGC‑TCGCGTAACGTGATCCCgg                                                                                                                                         <  8:71002/76‑1 (MQ=37)
                                                                                   tCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGC‑TCGCGTAACGTGATCCCgg                                                                                                                                         >  7:71002/1‑76 (MQ=37)
                                                                                   tCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCctc                                                                 >  8:104245/1‑148 (MQ=255)
                                                                                     aaCTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTAc                                                              >  4:248842/1‑149 (MQ=255)
                                                                                     aaCTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTAc                                                              >  8:122119/1‑149 (MQ=255)
                                                                                      aCTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTAc                                                              >  7:106647/1‑148 (MQ=255)
                                                                                       cTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACca                                                            <  3:248842/149‑1 (MQ=255)
                                                                                          cTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGCTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAg                                                         >  7:43284/1‑149 (MQ=255)
                                                                                                               gTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGt                                    >  8:193794/1‑149 (MQ=255)
                                                                                                                tCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTc                                   <  1:30501/149‑1 (MQ=255)
                                                                                                                  gAGGGTCATGCGCACTTCGCCCGGGAGTCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTcc                                  >  1:2580/1‑149 (MQ=255)
                                                                                                                     ggTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTgc                              <  8:106647/149‑1 (MQ=255)
                                                                                                                            cgcACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGt                        <  3:115156/148‑1 (MQ=255)
                                                                                                                               aCTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGccccc                                              <  3:33748/123‑1 (MQ=255)
                                                                                                                               aCTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGccccc                                              >  4:33748/1‑123 (MQ=255)
                                                                                                                                 ttCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAg                          >  8:13262/1‑141 (MQ=255)
                                                                                                                                 ttCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAg                          <  7:13262/141‑1 (MQ=255)
                                                                                                                                 ttCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACAGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCg                  >  6:112081/1‑149 (MQ=255)
                                                                                                                                  gtacccGGGC‑TCGCGTAACGTGATCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCgg                 >  1:76/4‑149 (MQ=255)
                                                                                                                                    gCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCg               >  5:65151/1‑149 (MQ=255)
                                                                                                                                               gcGTAACGTGATCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCgg      >  2:8360/2‑148 (MQ=255)
                                                                                                                                                  ttACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCGGCATg  <  2:249136/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCCGTGGGCACCAACTCATGGCTCAGGCCATGCTTGGCACAGGTCGCCTCGATCACCCCGCGCACCTCGGCAATCATCGCATTCAACTGCTCGCCCTCCAGGTGGCGGAAGTCGAGGGTCATGCGCACTTCGCCCGGGA‑TCACGTTACGCGAGCCCGGGTACGCCTGAAGGCAACCCACCGTGCCACAGGCGTGGGGCTGGTGGGCAAGCGCCGTGCGGTTGACCGCCTCTACCACAGCGGCTGCCCCCACCAGGGCGTCCTTGCGCAGGTGCATCGGCGTAGGGCCGGCATG  >  NC_002947/4544810‑4545103

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