Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,500,060 C→T pseudogene (475/942 nt) yjgX ← pseudogene fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,500,0600CT100.0% 55.5 / NA 19pseudogene (475/942 nt)yjgXpseudogene fragment
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/10);  total (9/10)

AAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTT  >  NC_000913/4499924‑4500168
                                                                                                                                        |                                                                                                            
aaaTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtac                                                                                                            <  2:387118/139‑1 (MQ=255)
  aTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTacac                                                                                                          <  2:294114/139‑1 (MQ=255)
               cGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                             <  2:98520/139‑1 (MQ=255)
                                 aGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGa                                                                                                 >  2:671576/1‑117 (MQ=255)
                                 aGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGa                                                                                                 <  1:671576/117‑1 (MQ=255)
                                 aGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGa                                                                           >  2:355955/1‑139 (MQ=255)
                                   cTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGaaa                                                                         >  1:590211/1‑139 (MQ=255)
                                   cTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGaaa                                                                         >  1:742360/1‑139 (MQ=255)
                                    tGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGaaaa                                                                        <  1:593607/139‑1 (MQ=255)
                                                                  cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                           <  1:332274/74‑1 (MQ=255)
                                                                  cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                           >  2:332274/1‑74 (MQ=255)
                                                                     gCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                         <  2:400916/121‑1 (MQ=255)
                                                                     gCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                         >  1:400916/1‑121 (MQ=255)
                                                                                ggCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGTATGTATCCGCCAGCGCACGTTCTGTTTcc                            <  1:290939/139‑1 (MQ=255)
                                                                                    aTTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCcatc                                                                                     >  1:713678/1‑78 (MQ=255)
                                                                                    aTTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCcatc                                                                                     <  2:713678/78‑1 (MQ=255)
                                                                                                  cTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCa          <  1:83205/139‑1 (MQ=255)
                                                                                                          ggTAATTATTCGGGTAGGAGCTAAGGCTGTTACACGGATTTGGATGGCAACCCATCATGTGCAGGAAAATTATCTTCGGAGCGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGtttt  >  2:265477/1‑139 (MQ=255)
                                                                                                                         aGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCAtgtg                                                                                <  1:576630/46‑1 (MQ=255)
                                                                                                                         aGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCAtgtg                                                                                >  2:576630/1‑46 (MQ=255)
                                                                                                                                aaGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTcc                            <  1:800282/91‑1 (MQ=255)
                                                                                                                                aaGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTcc                            >  2:800282/1‑91 (MQ=255)
                                                                                                                                        |                                                                                                            
AAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTT  >  NC_000913/4499924‑4500168

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