Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 734,833 A→C pseudogene (614/1434 nt) ybfO → conserved protein, rhs‑like YbfO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913734,8330AC97.7% 122.0 / ‑4.2 46pseudogene (614/1434 nt)ybfOconserved protein, rhs‑like YbfO
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base C (22/21);  minor base T (0/1);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGA  >  NC_000913/734723‑734968
                                                                                                              |                                                                                                                                       
gtTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGgtgt                                                                                                            >  1:1219992/1‑140 (MQ=14)
    ccgccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTacac                                                                                                        >  1:1152011/1‑140 (MQ=17)
      gccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGc                                                                                                      >  1:188587/1‑140 (MQ=17)
      gccgGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGc                                                                                                      >  1:506698/1‑140 (MQ=17)
                   cAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAAt                                                                                         <  1:960864/140‑1 (MQ=17)
                   cAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAAt                                                                                         <  2:39807/140‑1 (MQ=17)
                   cAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAAt                                                                                         <  2:344770/140‑1 (MQ=17)
                     gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATcc                                                                                       <  2:1219992/140‑1 (MQ=17)
                     gATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATcc                                                                                       <  1:963561/140‑1 (MQ=32)
                       tGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCAc                                                                                     >  2:485613/1‑140 (MQ=17)
                       tGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCAc                                                                                     <  2:451912/140‑1 (MQ=17)
                       tGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCAc                                                                                     <  2:1152011/140‑1 (MQ=17)
                          tCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTg                                                                                  >  1:1127977/1‑140 (MQ=17)
                          tCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTg                                                                                  >  2:424118/1‑140 (MQ=17)
                           cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGt                                                                                 >  2:824784/1‑140 (MQ=17)
                           cGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGGGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGt                                                                                 >  2:190708/1‑140 (MQ=17)
                                           gAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATcg                                                                 <  2:188587/140‑1 (MQ=17)
                                               tCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGc                                                             >  1:513646/1‑140 (MQ=17)
                                                    gCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGcc                                                        >  1:117252/1‑140 (MQ=17)
                                                        aCCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTg                                                    >  1:429421/1‑140 (MQ=17)
                                                           gggTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATcg                                                                 >  2:793647/1‑124 (MQ=17)
                                                           gggTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATcg                                                                 <  1:793647/124‑1 (MQ=17)
                                                           gggTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcg                                                 >  2:160453/1‑140 (MQ=17)
                                                                     ggAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCa                                       <  1:824784/140‑1 (MQ=255)
                                                                      gAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCAc                                      >  1:1018060/1‑140 (MQ=37)
                                                                       aaaGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACg                                     >  1:5687/1‑140 (MQ=37)
                                                                         aGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGa                                   <  1:72051/140‑1 (MQ=255)
                                                                          gccgccGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGaa                                  >  1:1224608/1‑140 (MQ=37)
                                                                           ccgccgCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAg                                 >  1:326589/1‑140 (MQ=37)
                                                                            cgccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgg                                <  2:1127977/140‑1 (MQ=255)
                                                                            cgccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgg                                <  1:160453/140‑1 (MQ=255)
                                                                              ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCGCCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                              <  2:1224608/140‑1 (MQ=255)
                                                                              ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                              <  2:513646/140‑1 (MQ=255)
                                                                              ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                              <  1:190708/140‑1 (MQ=38)
                                                                              ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                              <  1:202418/140‑1 (MQ=255)
                                                                              ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                              <  2:117252/140‑1 (MQ=37)
                                                                              ccgcTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAgggg                              <  2:1018060/140‑1 (MQ=255)
                                                                                 ctggctggCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCaa                           <  2:1238596/140‑1 (MQ=255)
                                                                                         cATCGTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATgg                   >  1:423971/1‑140 (MQ=255)
                                                                                             gTGCGGCCTGTCGGTGGGGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg               <  2:532155/140‑1 (MQ=255)
                                                                                             gTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg               <  2:5687/140‑1 (MQ=255)
                                                                                             gTGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcg               <  2:429421/140‑1 (MQ=255)
                                                                                              tGCGGCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTgcgc              >  1:925389/1‑140 (MQ=255)
                                                                                                 ggCCTGACGGTGGTGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGa           <  1:90957/140‑1 (MQ=255)
                                                                                                  gCCTGACGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGaa          >  1:1057404/1‑140 (MQ=255)
                                                                                                       aCGGTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGa     >  1:757851/1‑140 (MQ=255)
                                                                                                          gTGGCGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGa  <  1:587590/140‑1 (MQ=255)
                                                                                                              |                                                                                                                                       
GTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCGCTTGTCAGCACGGAAGGGGCAACAGAATGGTGCGCAGAATACGATGA  >  NC_000913/734723‑734968

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