Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 658,153 A→G pseudogene (123/195 nt) ybeM → pseudogene, C‑N hydrolase superfamily homology;putative enzyme; Not classified; putative amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913658,1530AG94.4% 167.2 / ‑0.7 54pseudogene (123/195 nt)ybeMpseudogene, C‑N hydrolase superfamily homology;putative enzyme; Not classified; putative amidase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (25/26);  total (26/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 = 8.99e-01

GTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAA  >  NC_000913/658058‑658248
                                                                                               |                                                                                               
gTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGc                                                                                             <  2:404719/100‑1 (MQ=255)
  tACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTGGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTg                                                                                           >  2:850469/1‑100 (MQ=255)
      tCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGgcgc                                                                                       >  1:756746/1‑100 (MQ=255)
        tgtgtgGGAAAAGAACGGTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGgcgcgc                                                                                     <  2:1229831/100‑1 (MQ=255)
         gtgtgGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGcgcgcg                                                                                    <  1:828582/100‑1 (MQ=255)
          tgtgGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCga                                                                                   <  2:305892/100‑1 (MQ=255)
              ggAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCgatgat                                                                               <  1:852178/100‑1 (MQ=255)
               gAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGAGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATc                                                                              >  1:314311/1‑100 (MQ=255)
                 aaaGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCa                                                                             >  1:702049/1‑99 (MQ=255)
                   aGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCatga                                                                          >  2:875037/1‑100 (MQ=255)
                    gAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCatgat                                                                         >  2:297922/1‑100 (MQ=255)
                       cGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCa                                                                      >  2:930258/1‑100 (MQ=255)
                        gCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAg                                                                     <  2:1124145/100‑1 (MQ=255)
                         cTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGa                                                                    <  2:202801/100‑1 (MQ=255)
                         cTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGa                                                                    >  1:742926/1‑100 (MQ=255)
                         cTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGa                                                                    >  1:1243252/1‑100 (MQ=255)
                           gagaTTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATc                                                                  >  1:356899/1‑100 (MQ=255)
                           gagaTTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATc                                                                  >  2:32187/1‑100 (MQ=255)
                             gaTTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTa                                                                >  2:725408/1‑100 (MQ=255)
                                 tgtgCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCgg                                                            <  1:326346/100‑1 (MQ=255)
                                 tgtgCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCgg                                                            <  2:194125/100‑1 (MQ=255)
                                   tgCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGtt                                                          <  1:1210824/100‑1 (MQ=255)
                                   tgCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGtt                                                          >  2:1049859/1‑100 (MQ=255)
                                    gCCTCGGTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGGTa                                                         >  2:986712/1‑100 (MQ=255)
                                      cTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTaa                                                        >  2:725261/1‑99 (MQ=255)
                                          ttGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAg                                                   <  1:1260276/100‑1 (MQ=255)
                                                gcgcAGGCGGCGGAAAACGACGCTTGGCTGTTTGCACTGCCGGAAGAGTGGCTGGCGCGCGATGATCGTTATTCAGATCTGTCGGTTAAATCAGCagag                                              >  2:1300751/1‑99 (MQ=255)
                                                      gcggcgGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTgg                                       <  1:32188/100‑1 (MQ=255)
                                                       cggcggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGa                                      <  2:158028/100‑1 (MQ=255)
                                                       cggcggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGa                                      <  2:658412/100‑1 (MQ=255)
                                                         gcggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAg                                    <  2:849232/100‑1 (MQ=255)
                                                         gcggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAg                                    <  2:1274209/100‑1 (MQ=255)
                                                          cggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAgg                                   >  1:167412/1‑100 (MQ=255)
                                                          cggAAAACGACGCATCGCTGTTTGCCCTGCCGAAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGaaag                                   >  2:164117/1‑98 (MQ=255)
                                                            gAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCg                                 >  1:909220/1‑100 (MQ=255)
                                                            gAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCg                                 <  1:987647/100‑1 (MQ=255)
                                                              aaaCGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGaa                               <  2:998489/100‑1 (MQ=255)
                                                                   acgCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCt                          >  1:1120372/1‑100 (MQ=255)
                                                                    cgCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTc                         >  2:387182/1‑100 (MQ=255)
                                                                    cgCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTc                         <  1:227381/100‑1 (MQ=255)
                                                                     gCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCg                        <  1:1097871/100‑1 (MQ=255)
                                                                     gCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCg                        <  2:715972/100‑1 (MQ=255)
                                                                      cATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCgg                       >  2:812033/1‑100 (MQ=255)
                                                                        tCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGAc                     <  1:328721/100‑1 (MQ=255)
                                                                         cGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACt                    <  1:461042/100‑1 (MQ=255)
                                                                          gCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACtt                   >  1:1226674/1‑100 (MQ=255)
                                                                             gTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTAc                <  1:812678/100‑1 (MQ=255)
                                                                                 gCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCg            >  2:62018/1‑100 (MQ=255)
                                                                                   ccTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCgag          >  1:286602/1‑100 (MQ=255)
                                                                                        cGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGt     <  1:329960/100‑1 (MQ=255)
                                                                                        cGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAg      <  2:744837/99‑1 (MQ=255)
                                                                                         ggAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTa    <  2:314308/100‑1 (MQ=255)
                                                                                           aaGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTaaa  >  1:252494/1‑100 (MQ=255)
                                                                                           aaGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTaaa  <  1:666197/100‑1 (MQ=255)
                                                                                               |                                                                                               
GTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAA  >  NC_000913/658058‑658248

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