Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 716,326 C→T 8.0% pseudogene (272/651 nt) ybfG ← pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913716,3260CT 8.0% 152.3 / 2.5 50pseudogene (272/651 nt)ybfGpseudogene
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (2/2);  total (25/25)
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

CAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAA  >  NC_000913/716184‑716469
                                                                                                                                              |                                                                                                                                               
cAGACTTTCCCAATTCAGCGTCTCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTtg                                                                                                                                             <  5:329984/147‑1 (MQ=255)
      ttCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCaa                                                                                                                                      >  8:222042/1‑148 (MQ=255)
       tCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAAc                                                                                                                                     >  1:379203/1‑148 (MQ=255)
                aGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCTATGGGTAAAGAACCCACGCCTACTAATTCAGGGTCTTGTGTCCAAATTAATGGa                                                                                                                             >  2:94789/1‑147 (MQ=255)
                aGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGa                                                                                                                             <  6:33908/147‑1 (MQ=255)
                aGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGa                                                                                                                             <  8:7541/147‑1 (MQ=255)
                aGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTACCCAAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCATCCGAATTCAGGGTGTTGTGGCAAACTTACTGGa                                                                                                                             >  4:312047/1‑147 (MQ=255)
                       cTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCAGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTgg                                                                                                                     >  8:217207/1‑148 (MQ=255)
                         ggCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGc                                                                                                                    >  3:339894/1‑147 (MQ=255)
                          gcgcAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGTGTTGTGTCCAACTTACTGGATCTCGTgg                                                                                                                     <  1:94789/145‑1 (MQ=255)
                          gcgcAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGTGTTGTGTCCAACTTACTGGATCTCGTgg                                                                                                                     <  3:312047/145‑1 (MQ=255)
                          gcgcAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGtt                                                                                                                                              >  1:375510/1‑120 (MQ=255)
                          gcgcAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGtt                                                                                                                                              <  2:375510/120‑1 (MQ=255)
                           cgcAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAAc                                                                                                                 <  2:379203/148‑1 (MQ=255)
                             cAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCg                                                                                                               <  8:316563/148‑1 (MQ=255)
                               aaaTAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCggg                                                                                                             <  7:222042/148‑1 (MQ=255)
                                          gCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGAt                                                                                                  <  4:339894/148‑1 (MQ=255)
                                                     cTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACtt                                                                                                                                   <  1:132541/104‑1 (MQ=255)
                                                     cTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACtt                                                                                                                                   >  2:132541/1‑104 (MQ=255)
                                                        ccccAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCt                                                                                            <  3:130107/140‑1 (MQ=255)
                                                        ccccAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCt                                                                                            >  4:130107/1‑140 (MQ=255)
                                                             aGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATAGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAg                                                                                                    <  3:137882/127‑1 (MQ=255)
                                                             aGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATAGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAg                                                                                                    >  4:137882/1‑127 (MQ=255)
                                                              gATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGAt                                                                                                  >  8:36904/1‑128 (MQ=255)
                                                              gATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGAt                                                                                                  <  7:36904/128‑1 (MQ=255)
                                                                   ccaccaTTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTat                                                                         >  4:99326/1‑148 (MQ=255)
                                                                      ccaTTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCACCTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTatcaa                                                                      <  3:295160/148‑1 (MQ=255)
                                                                       caTTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGCTGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTatcaat                                                                     >  3:285494/1‑148 (MQ=255)
                                                                          tAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGta                                                                                     <  2:223482/129‑1 (MQ=255)
                                                                          tAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACATACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGta                                                                                     >  1:223482/1‑129 (MQ=255)
                                                                                              gCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGa                                              >  5:220406/1‑148 (MQ=255)
                                                                                                 aCATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGAc                                             >  7:228220/1‑146 (MQ=255)
                                                                                                   aTCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgc                                         <  8:228220/148‑1 (MQ=255)
                                                                                                     cGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCCGACCTCCGGAGCCATAGACAAgcg                                        >  5:120970/1‑147 (MQ=255)
                                                                                                     cGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcgc                                       >  4:188307/1‑148 (MQ=255)
                                                                                                     cGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                        >  7:233679/1‑147 (MQ=255)
                                                                                                     cGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                        >  8:98048/1‑147 (MQ=255)
                                                                                                     cGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAgcg                                        >  3:302971/1‑147 (MQ=255)
                                                                                                                ggTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGt                             >  4:115661/1‑147 (MQ=255)
                                                                                                                 gTAAAGCCCACACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGa                           >  7:217364/1‑148 (MQ=255)
                                                                                                                    aaGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGAtt                         <  3:188307/147‑1 (MQ=255)
                                                                                                                      gCCCGCACGCCTGCGAATTCAGTGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttg                      <  7:365511/148‑1 (MQ=255)
                                                                                                                       cccGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGTGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgt                     <  7:217207/148‑1 (MQ=255)
                                                                                                                       cccGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgt                     <  3:99326/148‑1 (MQ=255)
                                                                                                                         cGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTttgttg                   <  4:302971/148‑1 (MQ=255)
                                                                                                                           cACGCCTGCTAATTCAGGGCGTTGTGTCCTACTTAGTTGATCTCGTGGCAACCGGGGGTGACGAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCTGAGCTCCGGAGCCATAGACAAGCGCACCCATCGGTGTTTTTGTTGTc                 <  6:220406/148‑1 (MQ=255)
                                                                                                                              gCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCAcc              >  4:234664/1‑148 (MQ=255)
                                                                                                                                       ttCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGa     <  4:285494/148‑1 (MQ=255)
                                                                                                                                        tCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGa     >  2:339840/1‑147 (MQ=255)
                                                                                                                                          aGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATaa  <  1:339840/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CAGACTTTCCCAATTCAGCGTATCTGGCGCAAAATAATAAACGCAGCCAGTATCTTCCCCAAGATCGCCACCATTAAGAGCAAAAATTACCCCTGCGACATCGTCAGCGATGGGTAAAGCCCGCACGCCTGCGAATTCAGGGCGTTGTGTCCAACTTACTGGATCTCGTGGCAACCGGGGGTGACCAGATCCTGCTATGCGTAACTAACCATTATCAATCAGCAGACCTCCGGAGCCATAGACAAGCGCACCCATCGGTGATTTTGTTGTCACCTGTAATTGATAA  >  NC_000913/716184‑716469

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