Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 675,601 Δ1 bp coding (947/978 nt) ybeQ ← Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913675,6010T.97.9% 188.9 / ‑3.0 47coding (947/978 nt)ybeQSel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (21/25);  total (21/26)
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

TTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGCAT  >  NC_000913/675474‑675732
                                                                                                                               |                                                                                                                                   
ttCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACAtctc                                                                                                                         <  2:717719/139‑1 (MQ=255)
ttCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACAtctc                                                                                                                         >  1:496533/1‑139 (MQ=255)
          gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCgg                                                                                                                <  2:638487/138‑1 (MQ=255)
          gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCgg                                                                                                                >  1:638487/1‑138 (MQ=255)
          gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                                               >  1:347904/1‑139 (MQ=255)
          gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                                               >  2:427803/1‑139 (MQ=255)
          gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                                               >  1:445979/1‑139 (MQ=255)
          gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                                               >  1:654645/1‑139 (MQ=255)
                      cGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACttt                                                                                                   <  2:227133/139‑1 (MQ=255)
                                tGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAgg                                                                                         >  2:343564/1‑139 (MQ=255)
                                  aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTc                                                                                       >  1:306572/1‑139 (MQ=255)
                                  aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTc                                                                                       >  2:216422/1‑139 (MQ=255)
                                        cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                 <  2:92343/139‑1 (MQ=255)
                                        cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                                                 <  2:723571/139‑1 (MQ=255)
                                             cAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTc                                                                            >  2:574145/1‑139 (MQ=255)
                                                     agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCg                                                                    >  2:244729/1‑139 (MQ=255)
                                                      gaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTCATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGa                                                                   <  2:298716/139‑1 (MQ=255)
                                                       aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATc                                                                     >  1:381061/1‑136 (MQ=255)
                                                       aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATc                                                                     <  2:381061/136‑1 (MQ=255)
                                                       aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                                  >  2:470565/1‑139 (MQ=255)
                                                       aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                                                  >  2:131697/1‑139 (MQ=255)
                                                        atatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATTCAGGTCGAGAGCCTGTCGGAGATCGACa                                                                 <  2:496533/139‑1 (MQ=255)
                                                          atGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAgg                                                               <  2:261576/139‑1 (MQ=255)
                                                              gCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTctt                                                           >  1:46400/1‑139 (MQ=255)
                                                               ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                          >  2:176183/1‑139 (MQ=255)
                                                                  aaTGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcctt                                                       <  1:176183/139‑1 (MQ=255)
                                                                    tGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTtct                                                     <  2:341617/139‑1 (MQ=255)
                                                                       aaTTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                  <  1:343564/139‑1 (MQ=255)
                                                                       aaTTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                  <  2:576376/139‑1 (MQ=255)
                                                                       aaTTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                  <  2:306572/139‑1 (MQ=255)
                                                                       aaTTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                  <  1:593071/139‑1 (MQ=255)
                                                                          tCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAg                                               <  1:131697/139‑1 (MQ=255)
                                                                              cGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGAt                                           <  2:162253/139‑1 (MQ=255)
                                                                              cGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGAt                                           <  2:46400/139‑1 (MQ=255)
                                                                                    ttGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATaa                                     <  2:445979/139‑1 (MQ=255)
                                                                                        ttttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                                 <  2:347904/139‑1 (MQ=255)
                                                                                          ttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                                 <  2:83053/137‑1 (MQ=255)
                                                                                          ttATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                                 >  1:83053/1‑137 (MQ=255)
                                                                                           tATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                  >  2:523223/1‑119 (MQ=255)
                                                                                           tATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                  <  1:523223/119‑1 (MQ=255)
                                                                                            aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg                             <  1:505420/139‑1 (MQ=255)
                                                                                            aTCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg                             <  2:709489/139‑1 (MQ=255)
                                                                                                 tACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGa                                            <  2:187219/119‑1 (MQ=255)
                                                                                                 tACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGa                                            >  1:187219/1‑119 (MQ=255)
                                                                                                              cAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCg           >  2:218750/1‑139 (MQ=255)
                                                                                                                     aCGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGc    >  2:383794/1‑139 (MQ=255)
                                                                                                                       gAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGCAt  <  2:117412/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                   
TTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGCAT  >  NC_000913/675474‑675732

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