Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 1,529,218 C→T pseudogene (1317/2037 nt) rhsE → pseudogene, Rhs family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,529,2180CT100.0% 46.3 / NA 22pseudogene (1317/2037 nt)rhsEpseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/9);  total (13/9)

CACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACC  >  NC_000913_3_bme_pgi/1529160‑1529363
                                                          |                                                                                                                                                 
cacCGGCGCGAAAAGCTCATTTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAg                                                         <  1:480416/149‑1 (MQ=255)
       gcgAAAAGCTCATCTTTGTCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTAcc                                                  <  1:209599/149‑1 (MQ=11)
         gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTTGCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGt                                                >  1:321817/1‑149 (MQ=11)
         gAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCCCCAGCCGTACCGt                                                >  1:218139/1‑149 (MQ=11)
                tCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATCCGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTACTGAGGAGAACCCGCATCACCTGCAACAGCCGTACCGTCTGCCAg                                         >  1:26893/1‑149 (MQ=11)
                  atctttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGACGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAggg                                       >  2:34433/4‑149 (MQ=11)
                    ctttATCACTGCGACCACCGGCGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGAAGAACCCGCATCACCTGCACCAGCCGTACCGTCCGCCAGGgca                                     >  2:527633/2‑149 (MQ=14)
                             tGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGCGTATGATGAAAGGGGCAACCAGCTTAATGAGAAGAACCCGCATCACCTGCACAAGCCGTACCGTCTGCCAGGGCAGTAGCatga                            >  2:478718/1‑149 (MQ=14)
                              gCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTACACCAGCCGTACCGTCTGCCAGGGCAGCAGCaagat                           >  1:50923/1‑149 (MQ=255)
                                 accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgag                        >  1:176221/1‑149 (MQ=255)
                                 accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgag                        <  2:433331/149‑1 (MQ=255)
                                     ccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggaga                    >  1:205138/1‑148 (MQ=255)
                                     ccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCCTACCGTCTGCCACGGCAGCAGCATACTGAGGCGt                    >  2:166216/1‑149 (MQ=14)
                                        gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGc                                            >  1:369429/1‑122 (MQ=255)
                                        gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGc                                            <  2:369429/122‑1 (MQ=255)
                                          gacTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg               <  1:34433/147‑1 (MQ=255)
                                           acTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAAGAGAACCCGCATCACCTGCACCAACCGTACCGTCTGCCAGGGCAGCAGCAGGAGGAGggggcggggc              >  2:321945/2‑149 (MQ=14)
                                            gggCCGTTGGGGCTTATGAGCGAAGAGGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCt             <  2:50923/147‑1 (MQ=17)
                                                cGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTAc         <  2:423245/149‑1 (MQ=255)
                                                       gcTTATCCGCGATGACGGCAATACGGCGTGGCGCGTGGAGTATGATGAATTGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAcc  <  1:228286/149‑1 (MQ=17)
                                                       gcTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                     <  2:313027/50‑1 (MQ=255)
                                                       gcTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAAt                                                                                                     >  1:313027/1‑50 (MQ=255)
                                                          |                                                                                                                                                 
CACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACC  >  NC_000913_3_bme_pgi/1529160‑1529363

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