Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,210,944 G→A pseudogene (1845/2001 nt) yehQ → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,210,9440GA100.0% 115.0 / NA 36pseudogene (1845/2001 nt)yehQpseudogene
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/16);  total (20/16)

ATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCC  >  NC_000913/2210853‑2211078
                                                                                           |                                                                                                                                      
atcaGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         >  1:903060/1‑139 (MQ=255)
atcaGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         >  2:634375/1‑139 (MQ=255)
     ttGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATgtgt                                                                                    <  1:482088/139‑1 (MQ=255)
                   ccAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGt                                                                      <  1:32522/139‑1 (MQ=255)
                                  tCAGCTGCTCCATAACGGCTAGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCCTAATACAGAAAGCGAg                                                       >  1:880609/1‑139 (MQ=255)
                                    aaCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                               <  2:374709/113‑1 (MQ=255)
                                    aaCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGa                                                                               >  1:374709/1‑113 (MQ=255)
                                     aCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         <  1:709863/102‑1 (MQ=255)
                                     aCTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         >  2:709863/1‑102 (MQ=255)
                                      cTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACGATTTCGTAATACAGAAAGCGAGGCAc                                                   >  1:249846/1‑139 (MQ=255)
                                       tGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCgg                                                                                                               >  1:611820/1‑78 (MQ=255)
                                       tGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCgg                                                                                                               <  2:611820/78‑1 (MQ=255)
                                       tGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCg                                                                                                  >  1:229206/1‑91 (MQ=255)
                                       tGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCg                                                                                                  <  2:229206/91‑1 (MQ=255)
                                         cTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTaa                                                                    >  1:130022/1‑119 (MQ=255)
                                         cTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTaa                                                                    <  2:130022/119‑1 (MQ=255)
                                                   gCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGTATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAAt                                      >  1:692996/1‑139 (MQ=255)
                                                      ggCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGaa                                   >  2:114966/1‑139 (MQ=255)
                                                             tCAGGAACAGAGTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGgt                            >  2:627907/1‑139 (MQ=255)
                                                                       agggCTATTGGGCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTaa                                                                    >  2:748547/4‑89 (MQ=255)
                                                                       agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCgg                                                                                                 >  1:844059/1‑60 (MQ=255)
                                                                       agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCgg                                                                                                 <  2:844059/60‑1 (MQ=255)
                                                                       agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTaa                                                                    <  1:748547/89‑1 (MQ=255)
                                                                       agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                  >  1:762788/1‑139 (MQ=255)
                                                                       agTGCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                  >  2:27376/1‑139 (MQ=255)
                                                                       agTGCTATTGGTCAGGCATAATTGCTGGCAAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTtttgcttt                  >  2:950916/1‑139 (MQ=255)
                                                                          gCTATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCg               <  2:880609/139‑1 (MQ=255)
                                                                            tATTGGTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGaa             <  1:634375/139‑1 (MQ=255)
                                                                                 gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         <  1:682381/58‑1 (MQ=255)
                                                                                 gTCAGGCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGcaca                                                                                         >  2:682381/1‑58 (MQ=255)
                                                                                  tCAGGTATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAatta       <  2:692996/139‑1 (MQ=255)
                                                                                     ggCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTa                                                                                                     <  2:1014392/42‑1 (MQ=255)
                                                                                     ggCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTa                                                                                                     >  1:1014392/1‑42 (MQ=255)
                                                                                      gCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            <  2:1129717/66‑1 (MQ=255)
                                                                                      gCATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACaa                                                                            >  1:1129717/1‑66 (MQ=255)
                                                                                       cATAATTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTccc  <  1:627907/139‑1 (MQ=255)
                                                                                           |                                                                                                                                      
ATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCTGGCGCTATCAGGAACAGAGTGCTATTGGTCAGGCATAGTTGCTGGCGAATGACCTCACCGCGGTGGGTTTTTATCGGCTGGCACATGTGTTGGGACAATTTCGTAATACAGAAAGCGAGGCACGGGTAGAAGCAATGAATAACGGTGTTTTGCTTTGCGAACAATTATTCCC  >  NC_000913/2210853‑2211078

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