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% 140.5 / NA 45pseudogene (351/453 nt)ydaYpseudogene, Rac prophage,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/23);  total (22/23)

TGGAATATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTG  >  NC_000913/1425847‑1426102
                                                                                                                                      |                                                                                                                          
tGGAATATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGcc                                                                                                                       <  2:412172/139‑1 (MQ=255)
   aaTATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAg                                                                                                                    <  2:717546/139‑1 (MQ=255)
            cAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTa                                                                                                           >  2:86349/1‑139 (MQ=255)
            cAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTa                                                                                                           >  1:360031/1‑139 (MQ=255)
            cAGAACGAGGGATCGCGCCGCTTAAGCAGACTCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAACAGCAGTAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTa                                                                                                           >  2:326423/1‑139 (MQ=255)
               aaCGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCg                                                                                                        >  1:1108471/1‑139 (MQ=255)
                 cGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGa                                                                                                      >  2:1210767/1‑139 (MQ=255)
                    gggATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGtt                                                                                                   <  1:1129112/139‑1 (MQ=255)
                      gATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCg                                                                                                 <  2:525897/139‑1 (MQ=255)
                               gctTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGaaa                                                                                        >  1:820224/1‑139 (MQ=255)
                                             cAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCa                                                                          <  1:570890/139‑1 (MQ=255)
                                                                    gCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAg                                                   >  2:1126799/1‑139 (MQ=255)
                                                                    gCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAg                                                   >  1:754852/1‑139 (MQ=255)
                                                                         agcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTg                                                                                                                         >  1:1073607/1‑64 (MQ=255)
                                                                         agcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTg                                                                                                                         <  2:1073607/64‑1 (MQ=255)
                                                                          gcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCgg                                                                                                <  2:1003424/88‑1 (MQ=255)
                                                                          gcagcCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCgg                                                                                                >  1:1003424/1‑88 (MQ=255)
                                                                             gccagcATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCAcc                                          <  2:1108471/139‑1 (MQ=255)
                                                                                  cATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggt                                                                                                                >  2:158347/1‑64 (MQ=255)
                                                                                  cATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCggt                                                                                                                <  1:158347/64‑1 (MQ=255)
                                                                                        ccAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGa                                                                                 >  2:122704/1‑89 (MQ=255)
                                                                                        ccAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGa                                                                                 <  1:122704/89‑1 (MQ=255)
                                                                                                 agagCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGa                                                                                 >  2:108244/1‑80 (MQ=255)
                                                                                                 agagCAGCAGAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGa                                                                                 <  1:108244/80‑1 (MQ=255)
                                                                                                 agagCAGCAGAAGCAGGATACGACAAAAAAGCCATAACTGCCCAGCGGAGGTTATCGGAGTTCCGGGGAGAAAGAGGTGAAAAGGCAAGATGGAATGGAGAAAAACAGAGGTTGCCACCAAAAAAGCAACCATAAATac                      >  2:855323/1‑138 (MQ=255)
                                                                                                     cagcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGAccc                  <  2:754852/139‑1 (MQ=255)
                                                                                                      agcagAAGACGGCATACGATATAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGATATAGAccc                  <  1:855323/139‑1 (MQ=255)
                                                                                                      agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCg                                                                                                 <  2:420942/59‑1 (MQ=255)
                                                                                                      agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCg                                                                                                 >  1:420942/1‑59 (MQ=255)
                                                                                                      agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCg                 <  1:1210767/139‑1 (MQ=255)
                                                                                                      agcagAAGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACAGGAACAAGAAATAGACCCg                 >  2:924249/1‑139 (MQ=255)
                                                                                                           aaGCAGGATACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCt            <  1:86349/139‑1 (MQ=255)
                                                                                                                   aTACGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg     <  2:820224/139‑1 (MQ=255)
                                                                                                                      cGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCg                                   >  1:392874/1‑106 (MQ=255)
                                                                                                                      cGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCg                                   <  2:392874/106‑1 (MQ=255)
                                                                                                                      cGATAAAAAAGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAgttg  >  2:920498/1‑139 (MQ=255)
                                                                                                                          aaaaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGtt                                                                                                   >  2:246135/1‑38 (MQ=255)
                                                                                                                          aaaaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGtt                                                                                                   <  1:246135/38‑1 (MQ=255)
                                                                                                                            aaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg     <  1:1225751/130‑1 (MQ=255)
                                                                                                                            aaaaGCCATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAg     >  2:1225751/1‑130 (MQ=255)
                                                                                                                                 ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGAccc                  >  1:707317/1‑112 (MQ=255)
                                                                                                                                 ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTa          <  2:372948/120‑1 (MQ=255)
                                                                                                                                 ccATAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTa          >  1:372948/1‑120 (MQ=255)
                                                                                                                                 ccATAACGGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGAccc                  <  2:707317/112‑1 (MQ=255)
                                                                                                                                   aTAACTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATa                       <  2:47713/105‑1 (MQ=255)
                                                                                                                                      |                                                                                                                          
TGGAATATACGGCAGAACGAGGGATCGCGCTGCTTAAGCAGAATCCAGAGTTGTACTGGCAGATCCTTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTG  >  NC_000913/1425847‑1426102

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: