Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 239,378 G→A pseudogene (273/273 nt) yafF → pseudogene, H repeat‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913239,3780GA100.0% 89.9 / NA 29pseudogene (273/273 nt)yafFpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/17);  total (12/17)

AGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGC  >  NC_000913/239243‑239512
                                                                                                                                       |                                                                                                                                      
aGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGtaataa                                                                                                                                     >  1:490923/1‑139 (MQ=25)
  ggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAataata                                                                                                                                   >  1:937567/1‑139 (MQ=37)
    aTACTGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGTAGCGGGCTTTCGTAATAATAcc                                                                                                                                 <  2:225193/139‑1 (MQ=255)
        ggCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGa                                                                                                                             <  2:937567/139‑1 (MQ=255)
               tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATAcc                                                                                                                                 <  1:977720/128‑1 (MQ=255)
               tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATAcc                                                                                                                                 >  2:977720/1‑128 (MQ=38)
                 cTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                                    <  1:841782/139‑1 (MQ=255)
                           tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAaca                                                                                                          >  2:858114/1‑139 (MQ=255)
                           tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAaca                                                                                                          >  2:7285/1‑139 (MQ=255)
                                   aTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAc                                                                                                      <  1:740374/135‑1 (MQ=255)
                                   aTGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAc                                                                                                      >  2:740374/1‑135 (MQ=255)
                                          ggTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacc                                                                                           <  1:759093/139‑1 (MQ=255)
                                             aTTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACcaca                                                                                        >  2:156819/1‑139 (MQ=255)
                                                   ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAccc                                                                                                    >  1:757445/1‑121 (MQ=255)
                                                   ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAccc                                                                                                    <  2:757445/121‑1 (MQ=255)
                                                                        gCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATg                                                             >  2:185371/1‑139 (MQ=255)
                                                                                 agcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACt                                                    <  2:490923/139‑1 (MQ=255)
                                                                                    cATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACtttt                                                 <  1:185371/139‑1 (MQ=255)
                                                                                       ggACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTtctc                                              <  1:156819/139‑1 (MQ=255)
                                                                                             aaaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacc                                                                                           <  1:283512/88‑1 (MQ=255)
                                                                                             aaaCTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacc                                                                                           >  2:283512/1‑88 (MQ=255)
                                                                                                   ccTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCa                                  >  2:211401/1‑139 (MQ=255)
                                                                                                      ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg                               <  1:858114/139‑1 (MQ=255)
                                                                                                      ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAAtg                               <  1:537315/139‑1 (MQ=255)
                                                                                                       gCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAc                                                                                                      <  2:353558/67‑1 (MQ=255)
                                                                                                       gCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAc                                                                                                      >  1:353558/1‑67 (MQ=255)
                                                                                                       gCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATc                                                                           <  2:288042/94‑1 (MQ=255)
                                                                                                       gCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATc                                                                           >  1:288042/1‑94 (MQ=255)
                                                                                                                                   cGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCAt                                                                            >  1:970370/1‑65 (MQ=255)
                                                                                                                                   cGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCAt                                                                            <  2:970370/65‑1 (MQ=255)
                                                                                                                                   cGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATgcgc  <  1:211772/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
AGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGACAGCTCCGCATCAATGCGC  >  NC_000913/239243‑239512

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