Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 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_0009131,529,2030CA88.9% 41.4 / NA 27pseudogene (1302/2037 nt)rhsEpseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (9/15);  minor base G (0/3);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAG  >  NC_000913/1529042‑1529468
                                                                                                                                                                 |                                                                                                                                                                                                                                                                         
ctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCg                                                                                                                                                          <  1:291886/275‑1 (MQ=12)
                  gaaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCat                                                                                                                                       >  1:359134/1‑276 (MQ=12)
                             aggccGCGGGAGGAGTGAAAGGCGAGGGGGGCTGGCGCAGGGGGGGGTGACGGTGGGGCAACTGGCCAGACATGTGGAGCCGGAATAAAAACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGt                                                                                                                             <  2:92109/272‑1 (MQ=11)
                             agacgggttgggCAGTGAAAGCGGGGCGTGGCGTGCGCAGGGGGGGGTGAAGGGGGAGCAAATGGCCAGAAAGGTGGAGCCGGAATACACACCGGGGCGAAAAGCACAGCTTTAGCACGGCGACCAACGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGt                                                                                                                             <  2:364728/264‑1 (MQ=11)
                                 cgcgtgggtgtggAAAGCCGGGGGTGGCTGGAGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCAGCTTTATCACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                        <  2:396415/264‑1 (MQ=11)
                                  gggTGGGCAGTGAAAGCCGGGGGTGGCGTGCGCAGGGCGGGCTGACGGGGGAGAAAATGGCCAGACAGGTGGGGCCGGAATAAACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                       <  2:377437/274‑1 (MQ=11)
                                    gTGAGCAGTGACAGCCGGTCGTGGCTTGCGCAGTGCGGGCGCACGGTGGAGCACCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTg                                                                                                                     <  1:146061/276‑1 (MQ=11)
                                         caggggaaGGCGGGCGGGGCTGGCGCAGTGCGGGCAGACGTGTGAGCAACTGGCCATACAGGTGGAGCCGGATGCCACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGTGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACt                                                                                                                 <  1:293775/269‑1 (MQ=11)
                                         cAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGCACCCGCATCACCTGCACCACCCGTACCGTCTGCCAGGGCAGCAGCATGATGCGGAGTCGGGGGGGTACt                                                                                                                 >  2:348356/1‑275 (MQ=11)
                                                    gggCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCAc                                                                                                      <  1:78989/275‑1 (MQ=255)
                                                       gtgggCTTGCGCAGTGCGGGCGGACGGGGGAGCAACTGGCCTGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCgg                                                                                                   <  2:404859/272‑1 (MQ=255)
                                                       cggggatggcgcAGTGCGGACTGACTGTGAGGAAAAGGGCAAGACAGGTGGACCAGGAATACACACCGGCGCGAAACGCTCATATTTATCACGGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCgg                                                                                                   <  2:37299/267‑1 (MQ=11)
                                                             ttGCGCAGTGCGGGCTGACGGTGGAGCAACGGGCCGGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCCCGGCGACCACCGGGGACTGCCGCTGGCGCTGATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGGCACCGGCACTACg                                                                                            <  1:20099/276‑1 (MQ=11)
                                                                           cTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                              <  1:16818/276‑1 (MQ=255)
                                                                            ggACGGTGGATCAACGGGCCAGACAGGTGGAGCCGGAATACACCCCGGCTCGAAAATCTCATCTCTATCACCGCGAGCACCGGGGACTGCCGCTGGCGTTTATCATCTAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                              <  1:150579/274‑1 (MQ=11)
                                                                               cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAATACGGCGTGGCGCGGGGCGTCTGATGAATGGGGCAACCAGCTTAATGAGGAGAGCCCGCATCACCTGAACCAGCAGTACATTATGCCAGGGCAGCCGCATGATGAGCAGTCGGGTCTGACCTATACCCGGCTCAGCCACTCGCATCCGTTGCAGggggggg                                                                          >  2:241407/1‑275 (MQ=11)
                                                                               cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCCGCAGACTGACGCGGAGTCGGGGCTGTACTATAACCGTACCCGGCACTACGACCCGTTGCCGGGGCgg                                                                           >  2:252601/1‑275 (MQ=11)
                                                                                           cggCCAGACAGGTGGAGCCGTAATCCACAGCGGCGAGAGAAGCTCATCTTAACCACTCCGACCGCCGAGGGCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATGCGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTAGCAGGGGCGGTATATCACCCCgg                                                              <  1:132351/275‑1 (MQ=14)
                                                                                                 caCATGTGGAGCCGAAATGACCGCCGGCTCTACAAGCGCATCTTTTGCACTGCCACCACCGGGGAGTGCCGCTGGCGCTGATCAGCGAAGACTGCAAGACGGCGTGGCGCGGGGAGTATGATGAATTGGTCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGa                                                        <  1:174309/275‑1 (MQ=14)
                                                                                                            ccGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCCCCGGCACTCCGATCCGTTGCAGGGCGGGTTTAACACCCCGTGCCCCGTAGGGTTgaga                                             >  2:404191/1‑276 (MQ=14)
                                                                                                                 aTACACGCCGTCGCGAAAAGCGCGTCTTTGACACTGTGACCGCCGGGGGATGCGGCGGGTGCGTATCAGGGAAGCCGGCAATACGGCGGGGGGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGGGGGGAGCCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtgg                                        <  2:332217/276‑1 (MQ=14)
                                                                                                                         cGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACAGTCACCGGCACTACGCTCAGTGGCAGGGGCGGTATATCCACCCGGACCCGATTGGGGTGGGCggtgggtgtgcgc                                >  2:330822/1‑271 (MQ=14)
                                                                                                                                        acctgtcgcagtGCGGCCACCGCGGGCGGCCGCGGGCGATGATCATCGAAGACGGCAATACGGCGGGGCGCGGGGAGTATGAGGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCg                  <  2:92010/264‑1 (MQ=17)
                                                                                                                                                    gcggacccgggggACTCCCGCTGGGGCCTAGCAGGGAAGACGGAGATGCGGCGTGGCGCGGGGAGTATGCTGGATGGGGCAACCATCTTAATGAGGAGAACCCGCATCCCCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATAc     <  2:375242/267‑1 (MQ=17)
                                                                                                                                                    gCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGGACTATAACCGTCACCGGCACTACGAGCCGTTGCAGGGGCGGTATATCACCCCGGACACGACTGGGGTGCGAGGTGGCTGGACTAGGTCTCAGTCTCCGTTGAATccccagc     >  2:328412/1‑272 (MQ=17)
                                                                                                                                                     cGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGCTTGGGTTGAGTGGTGGAAGGAATATGTAAACGTCTCCGTTGACTCCCCTACa    >  1:99090/1‑276 (MQ=17)
                                                                                                                                                       accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCCCCCCGGACCCGATTGGGTTGAGAGGTGGAGGGAATATGTATCAGTATCCGTGGATACCCATAAAAg  >  2:311360/1‑276 (MQ=255)
                                                                                                                                                                 |                                                                                                                                                                                                                                                                         
CTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAG  >  NC_000913/1529042‑1529468

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