Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,594,132 (C)8→10 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.C100.0% 77.7 / NA 22pseudogene (2/1461 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (11/11);  total (11/11)
*NC_0009131,594,1242.C90.9% 70.1 / 1.2 22pseudogene (2/1461 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (9/11);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGA  >  NC_000913/1593989‑1594234
                                                                                                                                        ||                                                                                                               
gTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑‑ccc                                                                                                             <  1:230367/139‑1 (MQ=255)
 tcatcaGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑‑cccc                                                                                                            <  2:455897/139‑1 (MQ=255)
      aGGGCATTATGTCGGGTAAAGGTGACGGTAATGGTCTTTACATGCCCCGGTTTCCCCCTTTCGCCAGTCATTTGATCCGTGCCGCCATTTTCCAGGTTGCCGCCGATAGTGGTTTTAGCGAAACGCGGAG‑‑‑ccccacccc                                                                                                       <  1:481017/139‑2 (MQ=255)
                  ccGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑‑CCCCCCCCAGGCATAACGtt                                                                                            >  2:543779/1‑139 (MQ=255)
                       ataGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCCAGGCATAACGTTCat                                                                                         <  2:213812/139‑1 (MQ=255)
                           ttGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCAc                                                                                    >  1:343246/1‑139 (MQ=255)
                              ccGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                                 <  2:60580/139‑1 (MQ=255)
                                gTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCa                                                                               >  2:482315/1‑139 (MQ=255)
                                           gTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGaactt                                                                    >  1:154493/1‑135 (MQ=255)
                                                 tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGC‑‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAc                                                             >  1:62678/1‑139 (MQ=255)
                                                 tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑‑ccccccccccaggcata                                                                                               <  2:66568/104‑10 (MQ=255)
                                                 tCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑‑ccccccccccaggcata                                                                                               >  1:66568/1‑95 (MQ=255)
                                                    ccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACAt                                                           >  1:592611/1‑139 (MQ=255)
                                                           ccccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                    >  2:54709/1‑139 (MQ=255)
                                                                cTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                               <  2:422231/139‑1 (MQ=255)
                                                                  tCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTAt                                             <  2:343246/139‑1 (MQ=255)
                                                                   cGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATg                                            <  2:592611/139‑1 (MQ=255)
                                                                                     tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                               <  2:41984/118‑1 (MQ=255)
                                                                                     tGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGtt                                               >  1:41984/1‑118 (MQ=255)
                                                                                                    ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc           <  1:482315/139‑1 (MQ=255)
                                                                                                       tGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                                                 >  2:330917/1‑98 (MQ=255)
                                                                                                       tGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                                                 <  1:330917/98‑1 (MQ=255)
                                                                                                          cgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACa                                                            >  2:13783/1‑84 (MQ=255)
                                                                                                          cgccgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACa                                                            <  1:13783/84‑1 (MQ=255)
                                                                                                             cgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGa  <  1:54709/139‑1 (MQ=255)
                                                                                                                          aGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                                 <  2:164767/47‑1 (MQ=255)
                                                                                                                          aGCGACACGCAGTGCCCCCCCCCCAGGCATAACGTTCATATCACCTg                                                                                 >  1:164767/1‑47 (MQ=255)
                                                                                                                                        ||                                                                                                               
GTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGA  >  NC_000913/1593989‑1594234

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