Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,323,166 A→G pseudogene (116/839 nt) phnE ← phosphonate ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,323,1660AG100.0% 77.4 / NA 24pseudogene (116/839 nt)phnEphosphonate ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/13);  total (11/13)

GGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCG  >  NC_000913/4323030‑4323300
                                                                                                                                        |                                                                                                                                      
ggTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCaa                                                                                                                                     <  2:384947/140‑1 (MQ=255)
     ggTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGcc                                                                                                                                <  1:1136902/140‑1 (MQ=255)
     ggTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGcc                                                                                                                                <  2:269162/140‑1 (MQ=255)
      gTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGCGCAACGGccc                                                                                                                               >  2:777091/1‑140 (MQ=255)
       tAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCa                                                                                                                              >  2:402633/1‑140 (MQ=255)
               ccACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGaa                                                                                                                      >  1:1077069/1‑140 (MQ=255)
                   tGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCga                                                                                                                  >  1:53036/1‑140 (MQ=255)
                     gCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGcc                                                                                                                                <  1:302464/124‑1 (MQ=255)
                     gCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGcc                                                                                                                                >  2:302464/1‑124 (MQ=255)
                               ggcgggAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTgcg                                                                                                      >  1:581901/1‑140 (MQ=255)
                                gcgggAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTgcgc                                                                                                     <  2:1037071/140‑1 (MQ=255)
                                                   gAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCgatggtga                                                                                  <  2:1077069/140‑1 (MQ=255)
                                                                     gccgTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCt                                                                <  2:581901/140‑1 (MQ=255)
                                                                                  aGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCt                                                                <  2:192115/127‑1 (MQ=255)
                                                                                  aGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCt                                                                >  1:192115/1‑127 (MQ=255)
                                                                                    cGTAAGCGGGGCCATTTCCGCGACCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCg                                                 >  1:950656/1‑140 (MQ=255)
                                                                                       aaGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGAt                                              <  2:53036/140‑1 (MQ=255)
                                                                                         gCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATaa                                            <  1:402633/140‑1 (MQ=255)
                                                                                              gCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGc                                                   <  2:165300/128‑1 (MQ=255)
                                                                                              gCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGc                                                   >  1:766812/1‑128 (MQ=255)
                                                                                              gCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGc                                                   >  1:165300/1‑128 (MQ=255)
                                                                                              gCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGc                                                   <  2:766812/128‑1 (MQ=255)
                                                                                                             tGCCACGAGACGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCgg                        >  1:1147196/1‑140 (MQ=255)
                                                                                                                   gagaCGACCAACACGGCGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCaa                  <  2:950656/140‑1 (MQ=255)
                                                                                                                                   cGAGTGCAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCg  >  2:867442/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                      
GGTGAGGTAATCCTGCCACTGGCTGAAATCGGGCGGGAAGAAGTCGGCGGCGAACGTCGCCATGTTGCCGCCGTCTTTGATCAGCGTAAGCGGGGCCATTTCCGCGCCCTGCCACGAGACGACCAACACGGCGAGTACAACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCG  >  NC_000913/4323030‑4323300

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