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% 74.6 / NA 24pseudogene (273/273 nt)yafFpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/13);  total (11/13)

AGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGAC  >  NC_000913/239243‑239494
                                                                                                                                       |                                                                                                                    
aGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGtaataa                                                                                                                   >  2:612936/1‑139 (MQ=18)
                           tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAaca                                                                                        >  1:358130/1‑139 (MQ=255)
                           tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAac                                                                                         >  1:50732/1‑138 (MQ=255)
                           tttGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAac                                                                                         <  2:50732/138‑1 (MQ=255)
                                 gAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAccc                                                                                  <  1:760414/139‑1 (MQ=255)
                                     gagaAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCAc                                                                              <  1:612936/139‑1 (MQ=255)
                                                aaaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGGCAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAAGCCCCACTCACCACAAcc                                                                   >  1:515909/2‑139 (MQ=255)
                                                      aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACt                                                             >  1:94859/1‑139 (MQ=255)
                                                                  tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                  <  2:72266/90‑1 (MQ=255)
                                                                  tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                  >  1:72266/1‑90 (MQ=255)
                                                                   aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                  <  2:790878/89‑1 (MQ=255)
                                                                   aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                  >  1:790878/1‑89 (MQ=255)
                                                                     gATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGcc                                              <  2:515909/139‑1 (MQ=255)
                                                                                                cTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTccc                   <  2:94859/139‑1 (MQ=255)
                                                                                                     tGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGAct                                                                                                         <  2:3966/48‑1 (MQ=255)
                                                                                                     tGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGAct                                                                                                         >  1:3966/1‑48 (MQ=255)
                                                                                                        cGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGaaa                                    <  2:405338/114‑1 (MQ=255)
                                                                                                        cGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGaaa                                    >  1:405338/1‑114 (MQ=255)
                                                                                                                 ttGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGAc  <  2:590923/139‑1 (MQ=255)
                                                                                                                  tGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                  >  2:218463/1‑42 (MQ=255)
                                                                                                                  tGCGGGGAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCg                                                                                                  <  1:218463/42‑1 (MQ=255)
                                                                                                                        gAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACt                                                                             <  2:240985/57‑1 (MQ=255)
                                                                                                                        gAGCGGGCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACt                                                                             >  1:240985/1‑57 (MQ=255)
                                                                                                                              gCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCt                                                                  <  1:611239/62‑1 (MQ=255)
                                                                                                                              gCTTTCGTAATAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCt                                                                  >  2:611239/1‑62 (MQ=255)
                                                                                                                                       |                                                                                                                    
AGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACAACCTAAACTCATCCGCATCCTGCCATGCCGGAAACTTTTCTCTATATTCCCGCAATGCTGCCATCGAC  >  NC_000913/239243‑239494

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: