Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 675,570 C→T *326* (TAG→TAA 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,5700CT95.5% 125.6 / ‑1.4 44*326* (TAG→TAAybeQSel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (18/24);  total (19/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 = 3.18e-01

GAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGC  >  NC_000913/675435‑675701
                                                                                                                                       |                                                                                                                                   
gAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTAc                                                                                                                                  <  2:178033/139‑1 (MQ=255)
      cTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTaaa                                                                                                                            >  1:6663/1‑139 (MQ=255)
                   cGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCAc                                                                                                               <  1:840493/139‑1 (MQ=255)
                    gAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGCCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACg                                                                                                              <  2:632826/139‑1 (MQ=255)
                         gtgGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAAttt                                                                                                         <  1:826796/139‑1 (MQ=255)
                                tGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGACTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTTCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTctgactg                                                                                                 >  1:882763/1‑139 (MQ=255)
                                                 gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                >  1:297457/1‑139 (MQ=255)
                                                 gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                >  1:300273/1‑139 (MQ=255)
                                                 gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                >  2:254439/1‑139 (MQ=255)
                                                         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                        <  2:772084/139‑1 (MQ=255)
                                                         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                        >  2:734953/1‑139 (MQ=255)
                                                         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                        <  2:753655/139‑1 (MQ=255)
                                                         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                        <  2:297457/139‑1 (MQ=255)
                                                                     tCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGtct                                                                                                                  <  1:375261/86‑1 (MQ=255)
                                                                     tCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGtct                                                                                                                  >  2:375261/1‑86 (MQ=255)
                                                                     tCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACa                                                            >  2:29114/1‑139 (MQ=255)
                                                                     tCTGATTTCGCACAGCAAGAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGGTACa                                                            >  1:612857/1‑139 (MQ=255)
                                                                         aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACt                                                                      >  1:804183/1‑125 (MQ=255)
                                                                         aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACt                                                                      <  2:804183/125‑1 (MQ=255)
                                                                         aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTc                                                        <  1:235828/139‑1 (MQ=255)
                                                                              gCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGagag                                                   >  1:18120/1‑139 (MQ=255)
                                                                               cacaGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGc                                                  <  1:29114/139‑1 (MQ=255)
                                                                                   gCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTGCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                        >  2:548725/1‑113 (MQ=255)
                                                                                   gCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                        <  1:548725/113‑1 (MQ=255)
                                                                                   gCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGt                                              <  2:466697/139‑1 (MQ=255)
                                                                                            agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTg                                                                          <  1:510872/102‑1 (MQ=255)
                                                                                            agaATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTg                                                                          >  2:510872/1‑102 (MQ=255)
                                                                                              aaTATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGAc                                   >  1:806761/1‑139 (MQ=255)
                                                                                               atatGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACa                                  >  2:295250/1‑139 (MQ=255)
                                                                                                 atGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAgg                                <  2:6663/139‑1 (MQ=255)
                                                                                                     gCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTctt                            >  1:225928/1‑139 (MQ=255)
                                                                                                      ccGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                           <  2:18120/139‑1 (MQ=255)
                                                                                                       cGAATGAAATTCCCCGGCCCTTGGTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACTATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcc                          <  2:612857/139‑1 (MQ=255)
                                                                                                        gAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGaa                                                                        <  1:815962/93‑1 (MQ=255)
                                                                                                        gAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGaa                                                                        >  2:815962/1‑93 (MQ=255)
                                                                                                         aaTGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcctt                        <  2:225928/139‑1 (MQ=255)
                                                                                                              aaTTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                   <  1:818337/139‑1 (MQ=255)
                                                                                                              aaTTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATATCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACTGATCATCCTTCTCCg                   <  2:622709/139‑1 (MQ=255)
                                                                                                                 tCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAg                <  2:300273/139‑1 (MQ=255)
                                                                                                                        cccTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACa         >  1:516773/1‑139 (MQ=255)
                                                                                                                          cTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATa       <  1:734953/139‑1 (MQ=255)
                                                                                                                            tGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                   >  2:492722/1‑125 (MQ=255)
                                                                                                                            tGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                   <  1:492722/125‑1 (MQ=255)
                                                                                                                               ttttATCTTTACAATAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc  <  1:63145/139‑1 (MQ=255)
                                                                                                                                    tcttTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcctt                        <  2:578371/112‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GAAGTGCTGACGCTATGCCCGAAATGTGGTCATGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGC  >  NC_000913/675435‑675701

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