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.95.7% 180.2 / 2.6 46coding (947/978 nt)ybeQSel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (21/23);  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.21e-01

TGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGC  >  NC_000913/675467‑675730
                                                                                                                                      |                                                                                                                                 
tGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGACTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTTCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTctgactg                                                                                                                              >  1:882763/1‑139 (MQ=255)
                 gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                                             >  2:254439/1‑139 (MQ=255)
                 gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                                             >  1:297457/1‑139 (MQ=255)
                 gCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGa                                                                                                             >  1:300273/1‑139 (MQ=255)
                         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                     <  2:297457/139‑1 (MQ=255)
                         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                     >  2:734953/1‑139 (MQ=255)
                         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                     <  2:772084/139‑1 (MQ=255)
                         gAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                     <  2:753655/139‑1 (MQ=255)
                                     tCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACa                                                                                         >  2:29114/1‑139 (MQ=255)
                                     tCTGATTTCGCACAGCAAGAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGGTACa                                                                                         >  1:612857/1‑139 (MQ=255)
                                         aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACt                                                                                                   <  2:804183/125‑1 (MQ=255)
                                         aTTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACt                                                                                                   >  1:804183/1‑125 (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                                                                                                     >  2:815962/1‑93 (MQ=255)
                                                                        gAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGaa                                                                                                     <  1:815962/93‑1 (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                                                <  1:492722/125‑1 (MQ=255)
                                                                                            tGTTTTTATCTTTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                >  2:492722/1‑125 (MQ=255)
                                                                                               ttttATCTTTACAATAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGc                               <  1:63145/139‑1 (MQ=255)
                                                                                                    tctttACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcctt                                                     >  1:578371/5‑112 (MQ=255)
                                                                                                    tcttTACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttcctt                                                     <  2:578371/112‑1 (MQ=255)
                                                                                                        tACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg                           <  2:821721/134‑1 (MQ=255)
                                                                                                        tACAGTAAATCTTCAGTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAgg                           >  1:821721/1‑134 (MQ=255)
                                                                                                                       gTCTCACGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTc       >  1:912965/1‑139 (MQ=255)
                                                                                                                            aCGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                <  2:351512/93‑1 (MQ=255)
                                                                                                                            aCGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCg                                                >  1:351512/1‑93 (MQ=255)
                                                                                                                            aCGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGc  >  2:751621/1‑139 (MQ=255)
                                                                                                                             cGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                        <  1:530252/84‑1 (MQ=255)
                                                                                                                             cGAATTTCCGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTcttc                                                        >  2:530252/1‑84 (MQ=255)
                                                                                                                                      |                                                                                                                                 
TGACCAGTTCCAGAGACGCCCGTTTGAGCCGTAATAGTCTGATTTCGCACAGCAATAAAGAGAATATGGGCCGAATGAAATTCCCCGGCCCTTGTTTTTATCTCTACAGTAAATCTTCAGTCTCACGAATTTCCTGACTGACATCTCTGGTTCCGGATGACTGAACTTTACGATACAGGTCGAGAGCCTGTCGGAGATCGACAGGTCTTCCTTCTCCGTAGTGATACATAATGGCCAGGTTATACCAGGCAGTGGCGTCATTGC  >  NC_000913/675467‑675730

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: