Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 297,699 A→G pseudogene (21/339 nt) yagP ← pseudogene, LysR family, fragment;putative regulator; Not classified; putative transcriptional regulator LYSR‑type

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913297,6990AG97.8% 154.0 / ‑4.0 45pseudogene (21/339 nt)yagPpseudogene, LysR family, fragment;putative regulator; Not classified; putative transcriptional regulator LYSR‑type
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (19/25);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

GGCAATTGTTATTCCTCCGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGACGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCAGCACGGAGAGT  >  NC_000913/297555‑297845
                                                                                                                                                |                                                                                                                                                  
ggCAATTGTTATTCCTCCGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGcc                                                                                                                                                >  1:501820/1‑149 (MQ=255)
 gCAATTGTTATTCCTCCGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGc                                                                                                                                             <  1:138949/151‑1 (MQ=255)
  cAATTGTTATTCCTCCGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGGAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGcc                                                                                                                                            <  2:292660/151‑1 (MQ=255)
    aTTGTTATTCCTCCGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGa                                                                                                                                          <  2:93753/151‑1 (MQ=255)
        ttATTCCTCCGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATac                                                                                                                                       >  2:404414/1‑150 (MQ=255)
             cctccGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTAt                                                                                                                                 <  2:104338/151‑1 (MQ=255)
             cctccGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTAt                                                                                                                                 <  1:95909/151‑1 (MQ=255)
              ctccGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATg                                                                                                                                >  1:93081/1‑151 (MQ=255)
                ccgccgGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATggg                                                                                                                              >  1:84993/1‑151 (MQ=255)
                ccgccgGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATg                                                                                                                                <  1:569466/149‑1 (MQ=255)
                  gccgGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGt                                                                                                                             >  1:340636/1‑150 (MQ=255)
                      gCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCTTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTa                                                                                                                            <  1:656266/147‑1 (MQ=255)
                        aaGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTccc                                                                                                                      <  1:323029/151‑1 (MQ=255)
                            gccagccTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCgg                                                                                                                  <  1:436574/151‑1 (MQ=255)
                             ccagccTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCgg                                                                                                                  >  2:370371/1‑150 (MQ=255)
                              cagccTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAg                                                                                                                <  2:211037/151‑1 (MQ=255)
                                             aCGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCActc                                                                                                  >  1:73736/1‑150 (MQ=255)
                                                 aGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCAc                                                                                             >  2:351996/1‑151 (MQ=255)
                                                                     cGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCa                                                                                                     <  1:626511/123‑1 (MQ=255)
                                                                     cGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCa                                                                                                     >  2:626521/1‑123 (MQ=255)
                                                                                      cGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCAt                                                        >  2:432866/1‑151 (MQ=255)
                                                                                      cGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCAt                                                        <  2:183391/151‑1 (MQ=255)
                                                                                       gAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATg                                                       >  2:645627/1‑151 (MQ=255)
                                                                                        aacgTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATg                                                       >  1:235564/4‑150 (MQ=255)
                                                                                        aaaGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATg                                                       <  1:25845/150‑1 (MQ=255)
                                                                                             tctcCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCTCGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAAtcc                                                  <  1:301121/150‑1 (MQ=255)
                                                                                                      ttcttcAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGAATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGAt                                        >  2:160156/1‑151 (MQ=255)
                                                                                                          tcAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGACAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATgg                                      >  1:114813/1‑149 (MQ=255)
                                                                                                          tcAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATgg                                      <  1:433068/149‑1 (MQ=255)
                                                                                                          tcAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGACGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATgg                                      >  1:594993/1‑149 (MQ=255)
                                                                                                            aaaTGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGgag                                    <  2:131617/149‑1 (MQ=255)
                                                                                                             aaTGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACg                                 <  2:298646/151‑1 (MQ=255)
                                                                                                              aTGGCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGt                                <  2:148024/151‑1 (MQ=255)
                                                                                                                ggCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGggcag                                                            >  2:186222/1‑121 (MQ=255)
                                                                                                                ggCCAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGggcag                                                            <  1:186221/121‑1 (MQ=255)
                                                                                                                  ccAGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACag                            <  2:544478/151‑1 (MQ=255)
                                                                                                                    aGCGATAAGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACag                            >  2:377507/1‑149 (MQ=255)
                                                                                                                          aaGGCGCAACCTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACt                    <  2:382328/151‑1 (MQ=255)
                                                                                                                                 aaCCGCCGGGGACTGGCGCCAGCCGATACACTGATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCa             <  2:114816/151‑1 (MQ=255)
                                                                                                                                    cTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGa                                         >  2:390334/1‑120 (MQ=255)
                                                                                                                                    cTCCGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGa                                         <  1:390330/120‑1 (MQ=255)
                                                                                                                                       cGGGGACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCAGCACgg       <  2:138513/151‑1 (MQ=255)
                                                                                                                                          ggACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCAGCACgg       >  2:646200/1‑148 (MQ=255)
                                                                                                                                           gACTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCAGCACGGagag   <  1:142087/151‑1 (MQ=255)
                                                                                                                                            aCTGGCGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCAGCACGGAGAGt  >  1:163514/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GGCAATTGTTATTCCTCCGCCGGCAAGAGCCAGCCTCAACATCAGACGCAGATCATTAGTCGTAATCTGCGGTTCAATCGCAAGGTCGAAAGTTCTCCCGTTTTCTTCAAATGGCCAGCGATAAGGCGCAACCTCCGGGGACTGACGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCAGCACGGAGAGT  >  NC_000913/297555‑297845

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑