Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,334,336 C→T pseudogene (67/4605 nt)
*208* (TAG→TAA
yfaS ←
yfaT ←
pseudogene, bacterial alpha2‑macroglobulin YfaS variant family, putative membrane protein
DUF1175 family protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,334,3360CT98.2% 172.4 / NA 55pseudogene (67/4605 nt)
*208* (TAG→TAA
yfaS
yfaT
pseudogene, bacterial alpha2‑macroglobulin YfaS variant family, putative membrane protein
DUF1175 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/28);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGG  >  NC_000913/2334210‑2334469
                                                                                                                              |                                                                                                                                     
tgctAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGa                                                                                                                           <  1:1101132/139‑1 (MQ=255)
  ctAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGaaa                                                                                                                         >  1:1220461/1‑139 (MQ=255)
            tgTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGa                                                                                                               >  1:45486/1‑139 (MQ=255)
               cAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATaa                                                                                                            <  2:45486/139‑1 (MQ=255)
                       aaaaaGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAAt                                                                                                    <  1:1339187/139‑1 (MQ=255)
                               tGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTgg                                                                                            >  2:559911/1‑139 (MQ=255)
                                    ccccGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGa                                                                                                                           >  2:192569/1‑103 (MQ=255)
                                    ccccGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGa                                                                                                                           <  1:192569/103‑1 (MQ=255)
                                         gCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAAtt                                                                                              >  2:1508908/1‑127 (MQ=255)
                                         gCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAAtt                                                                                              <  1:1508908/127‑1 (MQ=255)
                                         gCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATc                                                                                  >  2:485143/1‑139 (MQ=255)
                                             gcggcgCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTg                                                                              >  1:770375/1‑139 (MQ=255)
                                             gcggcgCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTg                                                                              >  1:734078/1‑139 (MQ=255)
                                              cggcgCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTgg                                                                             <  1:589484/139‑1 (MQ=255)
                                               ggcgCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTggg                                                                            <  2:734078/139‑1 (MQ=255)
                                                 cgcATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTa                                                                          <  2:1220461/139‑1 (MQ=255)
                                                       gTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                  <  2:277371/93‑1 (MQ=255)
                                                       gTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                  >  1:277371/1‑93 (MQ=255)
                                                       gTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATc                                                                    <  1:559911/139‑1 (MQ=255)
                                                         tGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCgg                                                                  <  2:1301757/139‑1 (MQ=255)
                                                           cTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                  <  1:356141/89‑1 (MQ=255)
                                                           cTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                  >  2:356141/1‑89 (MQ=255)
                                                                aggaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTg                                                                                       >  2:21099/1‑111 (MQ=255)
                                                                aggaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTg                                                                                       <  1:21099/111‑1 (MQ=255)
                                                                  gaagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTc                                                         >  2:523398/1‑139 (MQ=255)
                                                                   aagCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTcc                                                        <  2:770375/139‑1 (MQ=255)
                                                                     gCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTc                                                         <  2:436385/136‑1 (MQ=255)
                                                                     gCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTc                                                         >  1:436385/1‑136 (MQ=255)
                                                                     gCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCAt                                                      <  1:485143/139‑1 (MQ=255)
                                                                      cGAATCATCATCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCa                                                                                                      <  2:911991/90‑1 (MQ=255)
                                                                      cGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCa                                                                                                      >  1:911991/1‑90 (MQ=255)
                                                                        aaTCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAAt                                                                                                          <  1:1275130/84‑1 (MQ=255)
                                                                        aaTCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAAt                                                                                                          >  2:1275130/1‑84 (MQ=255)
                                                                                 cATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAgttgt                                          <  1:523398/139‑1 (MQ=255)
                                                                                   ttAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGttgttg                                        >  1:1503010/1‑139 (MQ=255)
                                                                                     aGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                              >  2:1271025/1‑131 (MQ=255)
                                                                                     aGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAg                                              <  1:1271025/131‑1 (MQ=255)
                                                                                           cccTGTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGTCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAgactga                                >  2:746719/1‑139 (MQ=255)
                                                                                               gTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                  >  1:1193001/1‑53 (MQ=255)
                                                                                               gTTCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTaaa                                                                                                                  <  2:1193001/53‑1 (MQ=255)
                                                                                                ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGc                                                                                                        <  2:1464755/62‑1 (MQ=255)
                                                                                                ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGc                                                                                                        >  1:1464755/1‑62 (MQ=255)
                                                                                                ttCCCACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTgcg                           >  1:1237173/1‑139 (MQ=255)
                                                                                                   ccACAAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTgcgcgc                        <  2:1450775/139‑1 (MQ=255)
                                                                                                     acaAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAAt                                                                                                    >  2:508745/1‑61 (MQ=255)
                                                                                                     acaAGAGACAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAAt                                                                                                    <  1:508745/61‑1 (MQ=255)
                                                                                                          agaCAGGGACAAGCACGCTGTTATCGCGGGAGAAAATTTAAAGGATAAATGCGAATGAAATTGGGATTGGAATCGTTGGGTATACATCGGGGGTTCTTCCAGGTAATAAGTTTGTGCAGACTGACTGTGAGCAATCTGt                 >  2:355054/1‑139 (MQ=255)
                                                                                                            aCAGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTg               >  1:508531/1‑139 (MQ=255)
                                                                                                              aGGGACAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTc             <  2:1237173/139‑1 (MQ=255)
                                                                                                                   cAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCg                                                                                 >  2:1439776/1‑66 (MQ=255)
                                                                                                                   cAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCg                                                                                 <  1:983765/66‑1 (MQ=255)
                                                                                                                   cAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCg                                                                                 <  1:1439776/66‑1 (MQ=255)
                                                                                                                   cAAGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCg                                                                                 >  2:983765/1‑66 (MQ=255)
                                                                                                                     aGCACGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTc      <  2:1503010/139‑1 (MQ=255)
                                                                                                                         cGCTGTTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCTCGCATTCCGTTGTCAGTTTTCGTgg  <  2:682495/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                     
TGCTAAAACTGCTGTCAGCAAGCAAAAAGAATGTTCCCCCGGCAGGCGGCGCATAGTTGCTGGAAGGAAGCGAATCATCAGCATTAGCAAGCCCTGTTCCCACAAGAGACAGGGACAAGCACGCTGCTATCGCGCCAGAAAATTTAAACGATAAATGCCAATGAAATTGGGATTGGAATCGTTGGGTATCCATCGGGTGTCCTTCCATGTCATAAGTTGTTGCAGACTGACTGCGCGCATTCCGTTGTCAGTTTTCGTGG  >  NC_000913/2334210‑2334469

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: