Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,425,980 G→A pseudogene (351/453 nt) ydaY → pseudogene, Rac prophage,Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,425,9800GA100.0% 87.5 / NA 29pseudogene (351/453 nt)ydaYpseudogene, Rac prophage,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/13);  total (16/13)

CAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATT  >  NC_000913/1425859‑1426115
                                                                                                                         |                                                                                                                                       
cAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTa                                                                                                                        >  2:178552/1‑139 (MQ=255)
cAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTa                                                                                                                        >  2:111746/1‑139 (MQ=255)
   aaCGAGGGATCGCGCTGCTTAAGCAGAATGCAGAGTTGTACTGGCAGATCCTTGCAAAAGGAGCCAGCATCGGCCAGGGTAAAGAGTAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCg                                                                                                                     <  1:220375/139‑1 (MQ=255)
                      tAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTAt                                                                                                                       <  1:797144/118‑1 (MQ=255)
                      tAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTAt                                                                                                                       >  2:797144/1‑118 (MQ=255)
                                 cAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCa                                                                                       >  1:31138/1‑139 (MQ=255)
                                              gCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggt                                                                                                                             <  1:31837/88‑1 (MQ=255)
                                              gCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggt                                                                                                                             >  2:31837/1‑88 (MQ=255)
                                                     cTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACt                                                                   <  2:31138/139‑1 (MQ=255)
                                                                    agcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAgg                                                               <  1:145937/128‑1 (MQ=255)
                                                                    agcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAgg                                                               >  2:145937/1‑128 (MQ=255)
                                                                                     agagCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                                    >  2:473759/1‑42 (MQ=255)
                                                                                     agagCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                                    <  1:473759/42‑1 (MQ=255)
                                                                                         cagcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGAccc                               >  2:17721/1‑139 (MQ=255)
                                                                                          agcagAAGCAGGATACGATAAAAAAGACATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACTCg                              >  2:586633/1‑139 (MQ=255)
                                                                                            cagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGtt                                                                                                                <  1:115035/55‑1 (MQ=255)
                                                                                            cagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGtt                                                                                                                >  2:115035/1‑55 (MQ=255)
                                                                                                      aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGa                                                                                                                   <  2:181321/42‑1 (MQ=255)
                                                                                                      aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGa                                                                                                                   >  1:181321/1‑42 (MQ=255)
                                                                                                      aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg                  <  1:111746/139‑1 (MQ=255)
                                                                                                      aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg                  >  2:406671/1‑139 (MQ=255)
                                                                                                               aaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTAc         >  1:16412/1‑139 (MQ=255)
                                                                                                                 aaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGGAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGc       >  2:517363/1‑139 (MQ=255)
                                                                                                                   gCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCt     >  1:645419/1‑139 (MQ=255)
                                                                                                                    ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTa    <  1:744324/139‑1 (MQ=255)
                                                                                                                    ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTa    >  2:744324/1‑139 (MQ=255)
                                                                                                                      aTAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTAtt  <  1:639566/139‑1 (MQ=255)
                                                                                                                      aTAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTAtt  <  1:45269/139‑1 (MQ=255)
                                                                                                                      aTAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTAtt  <  1:178552/139‑1 (MQ=255)
                                                                                                                         |                                                                                                                                       
CAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATT  >  NC_000913/1425859‑1426115

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