Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,132 (C)8→9 pseudogene (3856/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,594,1241.C96.0% 80.6 / ‑3.2 25pseudogene (2/1461 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (12/12);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT  >  NC_000913/1593991‑1594248
                                                                                                                                      |                                                                                                                             
catcaGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑ccccc                                                                                                                         >  1:26807/1‑139 (MQ=255)
                           gCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                                               <  1:395059/139‑1 (MQ=255)
                            ccGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGACCCGGGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACAAGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGc                                                                                              <  2:26807/139‑1 (MQ=255)
                                          tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCAATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACAAGCAGTG‑‑CCCCCCCCAGGCATAACGTTCATATCAACAGCAATGGTGACCTTTcc                                                                               >  2:685867/1‑139 (MQ=255)
                                                    ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cccccccccag                                                                                                                   >  2:689299/1‑90 (MQ=255)
                                                    ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cccccccccag                                                                                                                   >  2:605380/1‑90 (MQ=255)
                                                    ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cccccccccag                                                                                                                   <  1:605380/93‑4 (MQ=255)
                                                    ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cccccccccag                                                                                                                   <  1:689299/93‑4 (MQ=255)
                                                        tCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                                  >  1:769608/1‑139 (MQ=255)
                                                                           ttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCAt                                                                                                                <  2:679148/74‑1 (MQ=255)
                                                                           ttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCAt                                                                                                                >  1:679148/1‑74 (MQ=255)
                                                                                ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCAt                                                                                                                >  1:32728/1‑69 (MQ=255)
                                                                                ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCAt                                                                                                                >  2:382635/1‑69 (MQ=255)
                                                                                ccGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCAt                                                                                                                <  2:32728/69‑1 (MQ=255)
                                                                                ccGTGCCGCCATTCGCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCAt                                                                                                                <  1:382635/69‑1 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATgg                                                         <  1:476112/122‑1 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATgg                                                         >  2:476112/1‑122 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cccccccc                                                                                                                      <  1:622304/59‑1 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cccccccc                                                                                                                      >  2:622304/1‑59 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cc                                                                                                                            >  1:45224/1‑53 (MQ=255)
                                                                                   tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑cc                                                                                                                            <  2:45224/53‑1 (MQ=255)
                                                                                     ccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACAt                                                                         >  2:384314/1‑103 (MQ=255)
                                                                                     ccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACAt                                                                         <  1:384314/103‑1 (MQ=255)
                                                                                             ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                                               <  2:426753/73‑1 (MQ=255)
                                                                                             ctcCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                                               >  1:426753/1‑73 (MQ=255)
                                                                                                      gccgccGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc                                         >  1:364839/1‑118 (MQ=255)
                                                                                                      gccgccGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc                                         <  2:364839/118‑1 (MQ=255)
                                                                                                       ccgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                                             <  1:790458/97‑1 (MQ=255)
                                                                                                       ccgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                                             >  2:790458/1‑97 (MQ=255)
                                                                                                                        aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                    <  2:746778/73‑1 (MQ=255)
                                                                                                                        aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                    >  1:746778/1‑73 (MQ=255)
                                                                                                                        aGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAttt  >  2:587957/1‑139 (MQ=255)
                                                                                                                             cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                                    >  1:811107/1‑100 (MQ=255)
                                                                                                                             cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                                    <  2:811107/100‑1 (MQ=255)
                                                                                                                                      |                                                                                                                             
CATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT  >  NC_000913/1593991‑1594248

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