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% 198.8 / NA 61pseudogene (475/942 nt)yjgXpseudogene fragment
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/29);  total (32/29)

ATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTCTTGA  >  NC_000913/4499926‑4500194
                                                                                                                                      |                                                                                                                                      
aTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTacac                                                                                                                                    >  1:600762/1‑139 (MQ=255)
      tGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGAttt                                                                                                                              >  1:1435552/1‑139 (MQ=255)
                            tCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCa                                                                                                        >  2:5608/1‑139 (MQ=255)
                               aGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGa                                                                                                     >  1:199000/1‑139 (MQ=255)
                                gCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGAtt                                                                                                                               <  2:131492/112‑1 (MQ=255)
                                gCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGAtt                                                                                                                               >  1:131492/1‑112 (MQ=255)
                                 cTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGaaa                                                                                                   >  2:372903/1‑139 (MQ=255)
                                   ggCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAAt                                                                                                 <  1:1081788/139‑1 (MQ=255)
                                    gCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAAtt                                                                                                >  1:1343924/1‑139 (MQ=255)
                                                 tCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                                   <  1:5608/139‑1 (MQ=255)
                                                  ccAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGa                                                                                  >  1:739923/1‑139 (MQ=255)
                                                      tttACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTa                                                                              >  1:1269385/1‑139 (MQ=255)
                                                             aTGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                                       <  1:1411491/91‑1 (MQ=255)
                                                             aTGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                                       >  2:1411491/1‑91 (MQ=255)
                                                              tGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAg                                                                      <  2:600762/139‑1 (MQ=255)
                                                                 tGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                   <  2:1435552/139‑1 (MQ=255)
                                                                  ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGAt                                                                                 <  1:686172/124‑1 (MQ=255)
                                                                  ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGAt                                                                                 >  2:686172/1‑124 (MQ=255)
                                                                    cGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACg                                                                <  2:739923/139‑1 (MQ=255)
                                                                     gAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGt                                                               <  2:199000/139‑1 (MQ=255)
                                                                      aGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGtt                                                              <  1:561937/139‑1 (MQ=255)
                                                                                   tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCa                                                                                                        >  2:1161498/1‑84 (MQ=255)
                                                                                   tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCa                                                                                                        <  1:1161498/84‑1 (MQ=255)
                                                                                     tttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATg                                                                                                                         <  1:646968/65‑1 (MQ=255)
                                                                                     tttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATg                                                                                                                         >  2:646968/1‑65 (MQ=255)
                                                                                      ttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTaacaa                                              >  1:725374/1‑139 (MQ=255)
                                                                                        cGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGcc                                                                        >  2:1358617/1‑111 (MQ=255)
                                                                                        cGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGcc                                                                        <  1:1358617/111‑1 (MQ=255)
                                                                                           ggATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAAc                                                                                                                    >  2:1479938/1‑64 (MQ=255)
                                                                                           ggATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAAc                                                                                                                    <  1:1479938/64‑1 (MQ=255)
                                                                                               aCTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgaga                                                                                                                      >  2:968905/1‑58 (MQ=255)
                                                                                               aCTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgaga                                                                                                                      <  1:968905/58‑1 (MQ=255)
                                                                                                 tGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCAt                                   <  1:372903/139‑1 (MQ=255)
                                                                                                     ccAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCAc                                                                 <  1:1447493/105‑1 (MQ=255)
                                                                                                     ccAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCAc                                                                 >  2:1447493/1‑105 (MQ=255)
                                                                                                        ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCcatc                                                                                                               <  2:380973/56‑1 (MQ=255)
                                                                                                        ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCcatc                                                                                                               >  1:380973/1‑56 (MQ=255)
                                                                                                        ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                   <  1:296470/100‑1 (MQ=255)
                                                                                                        ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                                   >  2:296470/1‑100 (MQ=255)
                                                                                                        ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGtttt                            >  2:1251209/1‑139 (MQ=255)
                                                                                                        ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGtttt                            >  1:1024770/1‑139 (MQ=255)
                                                                                                             ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCAtgt                                                                                                           <  1:136479/55‑1 (MQ=255)
                                                                                                             ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCAtgt                                                                                                           >  2:136479/1‑55 (MQ=255)
                                                                                                             ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTaac                                                >  2:1259614/1‑114 (MQ=255)
                                                                                                             ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTaac                                                <  1:1259614/114‑1 (MQ=255)
                                                                                                             ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACggg                       >  1:1155612/1‑139 (MQ=255)
                                                                                                                ttCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                                   >  1:449058/1‑76 (MQ=255)
                                                                                                                ttCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                                   <  2:449058/76‑1 (MQ=255)
                                                                                                                ttCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaag                    <  2:725374/139‑1 (MQ=255)
                                                                                                                  cGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGaa                                                                                                                     <  2:8281/40‑1 (MQ=255)
                                                                                                                  cGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGaa                                                                                                                     >  1:8281/1‑40 (MQ=255)
                                                                                                                   ggttAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTAc                          >  2:758474/4‑130 (MQ=255)
                                                                                                                   ggttAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTAc                          <  1:758474/127‑1 (MQ=255)
                                                                                                                     ttAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                                       >  1:611484/2‑35 (MQ=255)
                                                                                                                     ttAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                                       <  2:611484/34‑1 (MQ=255)
                                                                                                                              aaGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGc           >  2:867991/1‑134 (MQ=255)
                                                                                                                              aaGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGc           <  1:867991/134‑1 (MQ=255)
                                                                                                                              aaGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTc      >  2:311530/1‑139 (MQ=255)
                                                                                                                               aGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCg                                                                                         >  1:242659/1‑55 (MQ=255)
                                                                                                                               aGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCg                                                                                         <  2:242659/55‑1 (MQ=255)
                                                                                                                                  cTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGACAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTCTTGa  <  2:1343924/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
ATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTCTTGA  >  NC_000913/4499926‑4500194

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑