Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 1,529,203 C→A pseudogene (1302/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,2030CA96.0% 54.3 / NA 25pseudogene (1302/2037 nt)rhsEpseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (12/12);  minor base G (0/1);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGG  >  NC_000913_3_bme_pgi/1529129‑1529348
                                                                          |                                                                                                                                                 
taaCTGGCGAGGGCGGTGGGGCCGGAATACACACCGGCGCGAAAAGTTCATTTTGATCACGGCGACCACCGGGGACTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAgtt                                                                         <  2:456404/148‑3 (MQ=11)
                                 ccGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGcc                                        <  1:715163/149‑1 (MQ=11)
                                      gtgAAAAGCTCAGCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTAcc                                   <  2:115745/147‑1 (MQ=11)
                                         aaaaGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTc                                <  1:45006/149‑1 (MQ=11)
                                              ctCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGccc                           >  1:17048/2‑148 (MQ=11)
                                              ctCATCTTTATCACTGCGACCACCGGGGACTACCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTACTGTGGAGAACCCGCATCCCCCGCACCAGCCGTACCGTCCGcct                           >  2:1284303/2‑148 (MQ=11)
                                                   ctttATCACTGCGACCACCGGGTACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgca                      <  1:702906/148‑1 (MQ=255)
                                                    tttATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcag                     >  1:993242/1‑149 (MQ=255)
                                                         cACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATAACGAAAAGGGCAACCAGCGTCATGAGGAGAACCCGCATCACCTCCACCAGCCGTACCGTCTGCCGGGgcagcagca                >  2:810889/1‑149 (MQ=14)
                                                          aCTGCGACCACCGGGGACTGCCGCTGGGGCTTATCCGCGAAGACGGCAATAGGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAGTGAGGAGAACCCGCATCACCTGCACCAGCCGTCTCGTCTGTCAGGgccgcggcag               >  2:902885/1‑148 (MQ=14)
                                                          aCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGGAGACGGCAATACGGCGTGGCGCGGGGAGTATGCTGAATGGGGCAACCAGCTCGATGAGGAGAACCCGCATCCCCTGCACCAGACGTCCAGTCTGGCGGGgcgccagcaa               >  2:811402/1‑148 (MQ=14)
                                                          aCTCCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGCCGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCat               <  1:241435/149‑1 (MQ=255)
                                                           cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCGATACGGCGTGGCGCGGGGAAGATGATGAATGGGGCAATCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatg              >  1:1197894/1‑149 (MQ=14)
                                                           cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTc                                >  1:162910/1‑131 (MQ=255)
                                                           cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTc                                <  2:162910/131‑1 (MQ=255)
                                                           cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatg              >  1:1264999/1‑149 (MQ=255)
                                                             gCGACCACCGGGGACTGCCGCTGGGGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat            <  2:1143412/149‑1 (MQ=255)
                                                              cGACCACGGGGGGCGGCCGGTGGCGCGTAGCAGCGTAGACGGTAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGGGGAGAACCCGCATCACCTGCAGCAGCCGTACCGTCGGTCAGGGCAGCAGCatgatg           <  2:1048880/149‑1 (MQ=14)
                                                               gACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGCATGGGGCAACCAGCTTAATGAGGAGAACCGGCATCACCTGCGCCAGCCGTACCGTCTGCAAGGGCTGCAGCAcgatga          >  1:504423/1‑149 (MQ=14)
                                                               gACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCAtgatga          <  1:811047/149‑1 (MQ=255)
                                                                  caccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgagga       >  2:1035159/1‑149 (MQ=255)
                                                                    ccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGt     <  2:1203644/149‑1 (MQ=255)
                                                                     cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                 <  2:1100791/88‑1 (MQ=255)
                                                                     cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagg                                                                 >  1:1100791/1‑88 (MQ=255)
                                                                       gggACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgg  <  1:1035159/149‑1 (MQ=255)
                                                                          |                                                                                                                                                 
CAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGG  >  NC_000913_3_bme_pgi/1529129‑1529348

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